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Abstract 
Several aspects of the clustering of earthquakes related to 
precursory seismicity patterns are examined. Using a simple 
definition for an aftershock, a self-similarity or scaling 
law in mainshock - aftershock clustering is proposed. The 
scaling law suggests that number of aftershocks does not 
depend on the magnitude of the mainshock if the aftershocks 
are counted in a fixed magnitude range below that of the 
mainshock. This scaling law reduces the number of 
observations necessary for statistical studies of aftershock 
Sequences and imposes a constraint on acceptable physical 
models of aftershocks. On the average, foreshocks seem to 
generate larger numbers of aftershocks than do mainshocks. 
The scaling law supports the hypothesis that the 
distribution of irregularities in material strengths ona 
fault zone shows no characteristic scale. 

The possibility of local clustering in earthquake 
sequences in Southern California with longer time and 
distance scales than the obvious "mainshock-aftershock" 
clustering is explored. To study the possibility of 
long-term clustering "foreshocks" and "aftershocks" are 
removed from the catalogue using a simple box car filter 
definition for the foreshocks and aftershocks. The lengths 
of the box car filters were dependent on the magnitude of 
the foreshock and mainshock. It is demonstrated that even if 
foreshocks and aftershocks are elimiminated through such a 


filtering process, the residual catalogue still shows weak 
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but significant clustering characteristics. This clustering 
is evident on time scales of the order of 3 to 6 years and 
distance scales of the order of 100 kilometers for residual 
events ("mainshocks") with magnitudes greater than or equal 
to 4, This observation supports suggestions that sequences 
of earthquakes are interdependent over longer times and 
distances than is widely accepted. The possibility of longer 
period clustering with periods greater than 6 years could 
not be explored due to the short length of the catalogue. 

Anomalous periods of clustering and anti-clustering in 
the sequence of residual events were juxtaposed with the 
occurrences of large events in the same region. Anomalous 
groups of residual events, Swarms or concentrated clusters 
of "mainshocks", preceded by quiescence were found to 
precede, most but not all, large earthquakes in the Southern 
California region. The quiescence > activation pattern 
observed here for mainshocks is Similar to precursory 
patterns observed by other investigators, although the time 
scale is in general longer. The mathematical technique of 
cluster analysis is shown to be useful in the study of 
anomalous seismicity patterns. 

Three proposed precursory seismic patterns termed 
"bursts of seismicity" are applied to retrospective 
prediction of large earthquakes which have occurred in the 
vicinity of the Cocos - North American - Caribbean triple 
junction. The three patterns, "bursts of aftershocks”, 


Swarms, and sigma, are based on abnormal clustering of 
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earthquakes in time, Space and energy. The patterns give 
respectable results, although the available data set is not 
of the highest quality. All occurrences of the patterns are 
close to the proposed location of the unstable triple 
junction. Success of the patterns may be due to the complex 


nature of the tectonics in the triple junction area. 
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1. Introduction 


1.1 Preliminary remarks 

Variations in seismicity patterns before large 
earthquakes have been studied by many investigators. These 
Studies serve dual purposes. Very little is known about the 
physical mechanism of earthquake failure and in particular 
about the nature of any inter-relationships between two or 
more earthquakes in a given region. Premonitory seismicity 
patterns place constraints on acceptable physical models 
which attempt to explain these phenomena. Also, seismicity 
patterns may be useful for earthquake prediction. 

In general the study of seismicity patterns can be 
approached differently depending on the purpose of the 
study. If the purpose is to develop an algorithm which is 
directly useful for earthquake prediction it is necessary to 
treat the seismicity patterns in as rigorous a manner as is 
Statistically possible. Completeness and homogeneity of 
earthquake catalogues, and uniformity of magnitude estimates 
are critical. Most studies have searched for anomalous 
seismicity before individual events without attempting to 
generalize the observed pattern. In this approach rigorous 
definitions of seismicity patterns and uniformity of the 
catalogue are less important. As a result most observed 


precursory seismicity patterns are only loosely defined. 
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Precursory seismicity patterns can be grouped into 
three basic categories, quiescence (including seismic gaps), 
activation (including swarms and foreshocks), and migration. 
Large earthquakes, at least those along simple plate 
boundaries, apparently repeat with intervals of a few tens 
to hundreds of years (for example Rikitake, 1976a). A number 
of investigators have tried to find periodities (recurrence 
intervals) and regularities in their occurrence. The main 
content of this chapter will be a review of observed 
precursory patterns and a review of studies of recurrence 
intervals for large earthquakes. 

An important feature of earthquake occurrence with 
great relevance to the failure process is the apparent self 
Similarity or scale invariance of earthquake catalogues 
(Kagan and Knopoff, 1980a). An aspect of this scale 
invariance is studied in chapter 2, a self similarity in the 
occurrence of aftershocks. Chapter 3 studies the clustering 
of mainshocks in Southern California. A precursory pattern, 
based on the quiescence and clustering of mainshocks, 
premonitory to large (M26.4) earthquakes in the region is 
proposed. Cluster analysis is applied to the study of the 
mainshock clustering. Chapter 4 applies the "bursts of 
seismicity" patterns proposed by Keilis-Borok et al. (1980b) 
to retrospective prediction of earthquakes in the region of 


thesCocos =" North American - Caribbean ®tuiple junction. 
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1.2 Implicit assumptions 

The direct use of premonitory patterns in earthquake 
prediction rests on several postulates. Foremost is the 
postulate that earthquakes in a catalogue do not follow a 
Poisson distribution. The occurrence of aftershocks is the 
most prominent deviation from Poisson behavior. With 
aftershocks removed there is some question as to whether the 
residual catalogue follows a Poisson distribution (Gardner 
and Knopoff, 1974). There is amibiguity in the exact 
definition of aftershocks. The question of longer period 
clustering with the obvious short term clustering removed is 
explored in Chapter 3. If the residual catalogue after short 
term clustering is removed is Poisson in character there is 
little hope of finding useful precursory seismicity 
patterns. 

Other postulates concern the nature of the failure 
process. The site of a large earthquake will not be 
physically isolated from surrounding regions. Presumably the 
state of stress, which occurs in the epicentral area of a 
major earthquake prior to failure, will have a significant 
effect on the state of stress in Surrounding regions, 
possibly to rather large distances. Since seismicity results 
from a release of accumulated strain; it seems reasonable 
that the nature of the seismicity in a region should depend 
on the amount of unreleased strain in that region. In 
particular, the seismicity when the average unreleased 


Strain is high might be expected to show some non-random 
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differences from the seismicity when the average unreleased 
strain is low. Presumably, the highest average unreleased 
strains in a region will be achieved just prior to the 
occurrence of a large earthquake in that region. 

In a given region large magnitude earthquakes are rare 
and the time spans of earthquake data sets are in general 
short compared to the return times for great earthquakes. 
This means that the statistical basis of any phenomenon 
associated with large magnitude earthquakes will at best be 
very weak for a given area. The direct use of premonitory 
patterns in earthquake prediction will likely fail unless 
there are systematic precursory patterns which do not vary 
Substantially in form from event to event and region to 


region. 


1.3 Some definitions 

The problem of searching for premonitory patterns can 
be summed up as a search for patterns in the background 
seismicity that precede the largest earthquakes (M2M,) ina 
region. The background seismicity is defined as all 
earthquakes with magnitudes less than M, and greater than 
some minimum threshold. The minimum threshold is usually 
determined by the detection capabilities of the seismic 
network being used. Aftershocks of the largest events are 


excluded from the background seismicity. 
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Some rules are needed to analyze the success of a given 
pattern. If a pattern is valid then the occurrences of the 
pattern should show a marked tendency to fall within some 
time window preceding the occurrences of the major events. 
The length of this time window may depend on several 
Parameters such as the regional rate of strain accumulation, 
the minimum magnitude of the major events, and the type of 
pattern being considered. Call the length of this precursory 
time window t.. A success for a pattern is then defined as 
the occurrence of the pattern which is followed by a major 
event within time t,.. A failure is defined as a major event 
which was not preceded by the occurrence of the pattern 
Withinetinest, 0A false alarm is defined as an occurrence of 
a pattern which is not followed by a major event within time 


Es 


1.4 Seismicity patterns and plate tectonics 

One of the most prominent features of world seismicity 
is that the vast majority of earthquakes occur along narrow, 
well defined, interconnecting belts. This feature was 
eaters by the global plate tectonic theory. The narrow 
seismic belts coincide with major plate boundaries. 
Relatively little deformation occurs in most plate 
interiors. 

The majouity of large (M275) ¥earthquakes occur at 


shallow depths along either convergent boundaries where 
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active subduction of one lithospheric plate is occurring or 
along strike-slip type boundaries. Large events account for 
most of the total seismic slip, and seismic energy release 


for these types of boundaries (Brune, 1968). 


1.4.1 Seismic Gaps 

Presumably strain energy is gradually accumulated at 
the plate boundaries due to the relative motion between 
plates. Accordingly, the largest earthquakes in the same 
region will usually be separated by many years. It has been 
shown that rupture zones for large earthquakes do not 
overlap appreciably (Mogi, 1968a; Sykes, 1971). Since a 
large portion of the relative plate motion along convergent 
and strike slip boundaries is accounted for by the slip 
during major earthquakes, it follows that, over a suitably 
long period of time, rupture zones of the largest 
earthquakes will tend to cover a seismic zone. Gaps in the 
Spatial distribution of these rupture zones are then likely 
locations for large earthquakes in the future. These gaps 
are termed seismic gaps or seismic gaps of the first kind 
(Mogi; 1979). The seismic gap hypothesis can be summed up as 
follows: segments of major plate boundaries that have not 
been the sites of large earthquakes for tens to hundreds of 
years are more likely to be the sites of future large 
earthquakes than segments that have experienced major 


rupture more recently. 
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Seismic gaps were first observed by Fedotov (1965) and 
successfully applied in the Kamchatka, Kurile Island and 
northeastern Japanese seismic regions. Mogi (1968a) 
independently applied the pattern to the Japanese and Alaska 
seismic zones, noting that gaps in Seismic activity in the 
Japanese seismic zone have been successively filled within 
several tens of years by great earthquakes without 
Significant overlap of their rupture zones. Subsequent 
papers (Sykes, 1971; Kelleher, 1970, 1972; Kelleher et al., 
1973: Ohtake et al., 1977: McCann et al., 1979) have applied 
the technique with great success to other seismic regions. 
McCann et al. (1979) gives a comprehensive review. 

The seismic gap technique provides estimates of the 
location and maximum likely size for some future large 
earthquakes, however, it does not constrain the time of 
occurrence to closer than decades (McCann et al., 1979). 
Another complication is the apparent existence of permanent 
seismic gaps (Kelleher and McCann, 1976,1977). Kelleher and 
McCann (1976,1977) found, contrary to what would be expected 
from plate tectonic considerations, sections of some island 
arcs have experienced very infrequent or no large 
earthquakes in the time spanned by present earthquake data. 
Most of these zones are at or near regions where groups of 
seamounts, aseismic ridges or other topographic features on 
the subducting plate intersect the subduction zone (Kelleher 
and McCann, 1976). They propose that in these regions the 


subduction process is modified due to the relative buoyant 
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material associated with the topographic features. The 
seismic gap method is not applicable in these regions 


(MeCannmet als ©1979). 


1.4.2 Frequency laws 

Most studies of seismic risk begin with an estimat2 of 
the frequency of occurrence of earthquakes with different 
magnitudes in the region of interest. These estimates are 
usually based on the historical records of seismicity. 
UnroOrtiunarely, nastorical records are too short to evaluate 
the frequency of occurrence with any confidence. Molnar 
(1979) suggests that intervals of the order of 10? years in 
most places and 10* years in some places are required. 

Gutenberg and Richter (1954) and others have shown that 
the number of earthquakes, N(M), with magnitudes greater 


than or equal to M is approximated well by 


log, .N(M) = a-bM . 


where a and b are constants for a given region. Due to the 
short lengths of time spanned by earthquake data sets, the 
constants a and b probably can not be accurately determined 
unless large regions are considered. If large regions are 
taken the analysis loses much of its usefulness. 

A useful parameter to characterize the size of an 


earthquake is the seismic moment: 


M, = yYU.oA 2 
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Wie Leia Gut nemshearemodulus, eA asethe fault area and u, 1s 
the average slip on the fault during the earthquake. The 
Summation of the seismic moments for individual earthquakes 
occurring on a fault can be used to estimate the average 
Slip on a fault (Brune, 1968). Let L be the length of the 
fault, and W the width. Brune (1968) showed that the average 
displacement on the fault U, after a sufficiently long 


period of time is given by. 


U = 2 (Mo/YLW) 3 

where the sum is over all earthquakes on the fault. The 
above assumes that no motion occurs by fault creep. Average 
rates of slip estimated tectonically or otherwise can be 
used to regionally constrain the frequency of occurrence of 
earthquakes with different magnitudes. Using this 
constraint, equation 3, and the Gutenberg Richter frequency 
Magnitude relation Molnar (1979), derives a frequency of 
OceuLLences Lawawath mores, Localeappluicab ilwuany. 


(1-8) ¥AUM, = 


N(M,) = ———-—__— 
: Mo(max yl 3) 


N(M,) is the number of events with seismic moments greater 
than or equal to M,, y is the average shear modulus on the 
faulfazone, —AsiSathesareapof ithey faulttezonesmvieisethe 
average slip rate excluding creep, £ is a constant which can 
be determined independently and M,(max) is the maximum 


possible seismic moment for an event in the region. Anderson 
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(1979) uses a Similar method with geological estimates of 
regional strain rates as the constraint. With only presently 
available data there are large uncertainies involved in 
estimating the values of the above parameters (8, y, Vv , 
etc.). Equation 4 probably results in only a modest 
improvement over eStimates which use historical data alone 
(Molnar, 1979). The Gutenberg Richter relationship may not 
be valid for largest magnitude events in some regions (Singh 
Stealer Toit) elterius suS athe caseuthe methods of Molnar 
(1979) and Anderson (1979) are not valid in these zones 


CSa ngs 6 tira lee 1 98.19)2 


1.4.3 Recurrence intervals 

Equations 1 and 4 are, at present, probably only useful 
in estimating mean seismicity patterns for large regions. 
Estimates of recurrence intervals (repeat times) for 
individual fault segments would be more useful. The 
recurrence of large earthquakes which rupture nearly the 
Same portion of a plate boundary have long been documented 
even before the introduction of plate tectonics. Early 
investigation of historical Japanese data (Imamura, 1928) 
indicated that large earthquakes in southwest Japan occurred 
repeatedly, at approximately the same location, with repeat 
times of 100 to 150 years. Similar regularities have been 
report in other regions by many authors (for example, Mogi, 
1968a; Fedotov et al., 1970; Kelleher, 1970, 1972; Sykes, 


1971+ Utsu, 1974: Rikitake, 1976a; Kelleher et al., 1973, 
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W972; Ando, 19797eSiehy 1978-"McNally and’ Minster, 1981: 
Singmeeteala, "1964 Sykes¥and Quittmeyer, heh iy he 

PMeEMSnOnemLImewspane Of* historical records for most 
regions make the statistical determination of repeat times 
at a location unreliable. For simple plate boundaries 


average repeat times ( T ) can be estimated by 


Tea / av 5 
where u is the average displacement which occurs 
seismically at the location, v is the relative plate 
Verocuivyera0OmOmI Sete elatio, of sensmicaslip to wotale slip 
(Sykes and Quittmeyer, 1981). The above assumes that 
relative plate velocities are uniform over long periods of 
time (Sykes and Quittmeyer, 1981), and that the majority of 
seismic slip in a region occurs during large earthquakes 
(Singh et al., 1981). Both assumptions seem to be reasonable 
(Sykes and Quittmeyer, 1981; Minster and Jordan, 1978: 
Brune, 1968). 

Equation 5 can be used to estimate return times from 
the history of previous events if a is constant locally over 
many earthquake cyles. This requires that some focal 
parameters of previous events are known from which u can be 
estimated. For example, u can be derived from the seismic 
moment, M,, uSing equation 2. Equation 5 then becomes 
Te=uM a vAVOlEOLMS Of equation > telatingeluslosothermrocal 
Parameters are given by Sykes and Quittmeyer (1981), Acharya 


ios weandmoincneet al. (198 1))/. 
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Time intervals between successive earthquakes at a 
given location are irregular enough that for most regions a 
knowledge of average recurrence interval does not greatly 
constrain future times of occurrence (Rikitake, 1976a; Ito, 
1980; Ellsworth et al., 1981). Ito (1980) suggests that much 
of the irregularity in recurrence intervals could be due to 
interaction between seismic events. 

The stress drop in an earthquake is proportional to the 
coseismic slip (for example, Chinnery, 1969). Using this 
Shimazaki and Nakata (1980) suggest a simple model for 
estimating recurrence intervals of individual events. Let T, 
be the shear stress at which the earthquake initiates, and 
let tT, be the final value after the earthquake. If the rate 
of accumulation of stress is constant in time, the cyclic 
nature of stress relief associated with a sequence of 
earthquakes will be related to the manner in which 1T, and T, 
change from event to event (Shimazaki and Nakata, 1980). A 
Strictly periodic sequence of events would result if T, and 
T, are both time independent. There will be no regularity if 
bothim, andvteavaryeinstime*: (itis evariesswithetameuwhw cma, 
remains constant a sequence of events results in which the 
time of occurrence of an event would depend on the amount of 
coseismic displacement of the preceding event. Shimazaki and 
Nakata (1980) call this the "tame=predictable™ models 272i 7, 
remains constant and tT, varies "Slip-predictable" recurrence 
occurs. The amount of slip of subsequent events will be 


DECOOnC1 Ona ImcOmthemtotal time welapsed tromethe: preceding 
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earthquake, the time of occurrence will be uncertain. The 
time-predictable and slip-predictable models are illustrated 
Mek Fe sWepbrat=: ahecie. 

Several studies find support for the time-predictable 
model. Bufe et sak, (CI) found that the time interval 
between small shocks along a segment of the Calaveras fault 
in California depended on the displacement of previous 
shocks. Similar findings are reported for three sequences of 
earthquakes and earthquake related phenomena in Japan 
(Shimazaki and Nakata, 1980). Sykes and Quittmeyer (1981) 
found good agreement between calculated and actual 
recurrence times for sequences of events along many simple 
boundaries. They also use the time-predictable and 
slip-predictable models to calculate a of equation 5. The 
time-predictable model gave realistic values whereas the 
Slip-predictable model did not. 

If the slip-predictable model is valid the occurrence 
times of future large events along simple plate boundaries 
can be greatly constrained provided the coseismic slips of 
several previous events at the location are known (Sykes and 


Quittmeyer, 1981). 


1.5 Quiescence 
Quiescence refers to a Spatial-temporal reduction (gap) 
of seismic activity in comparison to the normal activity. In 


the literature, quiescence is often confused with seismic 
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STRESS 


CUMULATIVE COSEISMIC 
DISPLACEMENT 





(a) TIME (b) 


Piguremiart Sean sillustration of therl"time-predictablemma) 
and the "slip-predictable” (b) models of recurrence (after 
Shimazakiwand "Nakata, }1980)2°inemodel "(aye thesstress at 
which the earthquake initiates Tt, is constant for all shocks 
Ducmonestinalestness 1, Valles. [hiSeresultsaiieampaLrtennmain 
which the time to the next event ( T )is determined by the 
SlipeGeUe) Of thee last event. In model (bjmthemstiescmat 
which the earthquake initiates varies while the final stress 
remains constant. This results in a sequence in which the 
Slip of the next event in Dproportionalsitosthemtamemirconmmene 
last event. V is the relative plate velocity, and a is the 
PeplOsOmeseloniCucilDELOuLOta lois 





gaps. The term seismic gap will be reserved for gaps in the 
Spatial A eenmbueionwet rupture zones of the largest 
earthquakes in a seismic belt. Mogi (1979) refers to 
Spatial-temporal gaps in the seismicity pattern as seismic 
gaps of the second kind. 

A reduction of smaller magnitude seismic activity ina 
focal region of a future large earthquake is commonly 
observed. This pattern was first noted by Inouye (1965) for 
several large shallow earthquakes in Japan. Kelleher and 
Savino (1975) demonstrated that seismic gaps are also often 
gaps for smaller magnitude activity and that such gaps 
commonly persist until the time of the future large event. 
Other studies have reported the pattern for earthquake with 
magnitudes as low as 5 in virtually all tectonic settings 
including intraplate seismic zones (for example, Mogi, 
iSocapemivova me bDOnOVIUKEet al. 2 19/0" McEVillyeetual., 967+ 
Hvis Ole Silapmelord One 1977 C*mSekl ya = 197 7s Engdalleand 
Kisslinger, 1977; Mizone et al., 1978; Kelleher and Savino, 
1975; Ohtake et al., 1977; Ishida and Kanamori, 1977; Wei et 
al., 1978; Kristy and Simpson, 1980; Yamashina and Inoue, 
1979: ‘Kanamori, 1981; Wyss and Habermann, 1979; and Khattri 
and Wyss, 1978). Kanamori (1981) gives a tabulation of 
earthquakes for which the pattern has been observed, noting 
that preseismic quiescence appears to to be the most 
commonly observed precursory pattern. 

The above studies suggest that seismic activity in the 


region of the eventual rupture zone of a large earthquake 








. - : sap Sinaise wt er See ag . 
? sepeve i a4 = i wt ic rel 4oo) 16s) Ae a 
» 34 ‘ pro F | a slaodetrs to supe? ae 

~»! ae 7 . 7 ~a, Lau. donc 
*) SGau. 










eirren) fi é : ae te 





sativa! t bmn 

7 oan Meh? 8. vkG? oc Ras 

See LOC SEEE, Pail Hh ceye 2 42) 8 am 
aay 


7 F ie ys oe ; : : 
on 7 






7 _ 


- —_ 
7 


decreases more or less abruptly before the major event. The 
nature of this decrease is poorly understood and seems to 
vary Significantly from event to event. Mogi (1979) notes 
that the reduced levels of seismicity may extend only 
several magnitude ranges below the future large event; 
smaller earthquakes may or may not show normal activity. 
Also, for an individual event the spatial extent and time 
interval of the quiescence may vary significantly depending 
on the magnitude interval used to define the seismic 
activity. Ishida and Kanamori (1978) and Kanamori (1981) 
note that there are many quiet zones and quiescent periods 
which are not followed by a large earthquake. Quiescence 
alone is not a reliable indicator of a future large 


earthquake. 


1.6 Activation 

Activation refers to a spatial-temporal increase in the 
seismic activity of a region in comparison to normal 
activity. Activation can be manifest through increases in 
seismicity rates, increases in the total seismic energy 
poneaser or both. Reported increases are regional or local 
in nature. Swarms and foreshocks are examples of local 


aGtuvatlon. 
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1.6.1 Foreshocks 

Foreshocks are the most obvious premonitory phenomenon 
preceding major earthquakes. Foreshocks are events or swarms 
which are close in space and time to the future mainshock. 
However, there iS no unambiguous, widely accepted definition 
of what is meant by "close in space and time". Some 
earthquakes are preceded by very clear increases in the 
number of smaller earthquakes in their epicentral area a few 
days or weeks before they occur. Such events can 
unambiguously be called foreshocks. Such a sequence was very 
instrumental in the successful prediction of the 1975 M=7.3 
Haicheng, China earthquake (Scholz, 1977). In other cases 
increased activity may be weaker or more spread in time, yet 
the events may be causally related to the mainshock and so 
should also be considered as foreshocks (Kanamori, 1981). 

Jones and Molnar (1976), defining foreshocks as events 
occurring within 100km and 40 days of the mainshock, found 
that 44% of all large M 2 7 shallow earthquakes in the world 
from 1950 to 1973 were proceeded by foreshocks large enough 
to be teleseismically reported, and that 21% of the large 
earthquakes that occurred in China between 1900 and 1949 had 
foreshocks large enough to be noted in local records. They 
note that the percent of large events with foreshocks is 
probably much higher than 44%. Many events with locally 
reported foreshocks (i.e. the 1975 Haicheng, China 


earthquake) had no foreshocks recorded teleseismically. 
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Many examples of possible broad sense foreshock 
activity have been reported. Mogi (1969a) and Kelleher and 
Savino (1975) show that for great earthquakes, seismic 
activity prior to failure tends to cluster around the future 
epicenter. Other examples of tight clustering of activity 
around future epicenters are given by Engdahl and Kisslinger 
(1977), Ishida and Kanamori (1978), and Fuis and Lindh 
(1979). A tabulation of many events for which foreshocks 
have been reported can be found in Rikitake (1975b, 1979) 
and Kanamori (1981). 

The unambiguous identification of foreshocks sequences, 
as in the case of the 1975 Haicheng earthquake, would 
provide an extremely useful short term predictor for many 
earthquakes. Unfortunately, the occurrence of a foreshock 
sequence is generally not obvious Gaeit after the mainshock 
has occurred. Laboratory studies and theoretical 
considerations suggest that the ratio of smaller fractures 
to larger fractures should decrease as effective 
differential stress in an area increases (Scholz et al., 
1973). This has lead to the postulate that foreshock 
sequences are characterized by smaller b values in the 
Gutenberg-Richter relation (equation 1) than are aftershock 
sequences or the background seismicity. Decreases in b have 
been reported for several foreshock sequences (for example, 
Suyehiro,1966; Papazachos et al., 1967; Suyehiro and Sekiya, 
Nod 2eeesCHhOLZ Giada 19/3) 3Papazacnos, 975° WySs and sLee, 


1973; HasSagawa et al., 1975). A large number of events are 
















ox AD 
(3462425) @anet * ‘ae jae 4 1S Teeaeee Ss 
ecetiat Sor | \ got naosane ased eve ce < 
ihe at Terre son ie ae cor’ caleaaaaE 


= a F * = 4 on +e _ 
a gogsu2 aig ft 7 Oe Cer ; | + 765 6g Bee. i?) 


bert tip rns, ) ceio 2682) Seta 





ha { itTei) 
— ; - * j . (ete) 7 
f¢ : - 3s 1995 @veng 


eet One 


nop 44m . 
agin 2 | a en: hhh eieg a 


@rac oF + 





a 1s + j en 2 | oupeteit6 
> Le. ie alusa 724 Ltr tere 






Se 6 G' 


- 


Eigener Sb Rens 4° sopdapps) 0un.1e66y tei eg 
va ae bruni ideo i> ae 


ait dope maar re ia 


api. — 7 - 










7 


needed for the reliable determination of b. Unfortunately, 
the number of recorded events in most foreshock sequences is 
small. The usefulness of monitoring local b values as a 
function of time to detect foreshock sequences is debatable. 
Monitoring this would require much more sensitive seismic 
networks than presently exist in most regions. Also, it is 
notaclearsthateequationdieissapplicable toslocaleseismicity 
pateenns: 

On the basis of a statistical study, Yamashina (1981) 
Suggests that a sequence of events occurring within a short 
time interval of each other is more likely to be a foreshock 
sequence if the magnitude difference between the largest and 
second largest event is small. In ehaptere2eotethis ithesis- 
it 1S Suggested that, on the average, foreshocks are 
characterized by a greater rate of production of 
aftershocks. The above effects, if present, are weak. Their 
usefulness in identifying foreshock sequences is doubtful. 

Ishida and Kanamori (1980) found that waveforms of 
foreshocks of the 1971 San Fernando and 1952 Kern County 
earthquakes are characterized by a concentration of energy 
at higher frequencies relative to earlier events. They 
suggest that such effects, although subtle, may be useful 


for monitoring the stress state of selected regions. 
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1.6.2 Swarms 

A swarm can be loosely defined as a sequence of 
earthquakes which are close to each other in space, time, 
and magnitude (i.e. a sequence of tightly clustered 
earthquakes similar to an aftershock sequence but lacking 
the mainshock). McNally (1977) found that clusters of small 
earthquakes occurred close to the future epicenters of 
several moderate size earthquakes in California two to ten 
years before the mainshock. Ohtake (1976) and Sekiya (1977) 
reported similar findings for some larger Japanese 
earthquakes. Evison (1977a,b,c) finds that for many 
earthquakes worldwide which are characterized by quiescence, 
the quiescent period was preceded by a characteristic swarm 
of minor activity. He suggests that identification of swarms 
which are followed by quiescence may be useful for 


earthquake prediction. 


1.6.3 Bursts of seismicity 

Three related seismicity patterns, precursory to the 
largest earthquakes in a region, are proposed by 
Keilis-Borok et al. (1980a,b,c). These three patterns, 
collectively termed "bursts of seismicity", consist of the 
abnormal clustering of earthquakes in time, energy and space 
before a major earthquake. 

The simplest of the three patterns is pattern B, 
"bursts of aftershocks". Pattern B consists of a medium 


magnitude mainshock with an anomalous number of aftershocks 
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concentrated at the beginning of the aftershock sequence. 
Pattern S consists of a "Swarm" of mainshocks where a 
"Swarm" is defined as a group of moderate events 
concentrated in space and time and occurring during a time 
interval when the overall seismicity is not below average. 
The final pattern, Sigma, consists roughly of an increase in 
the cummulative seismic energy released to the 2/3 power in 
a sliding time window. Sigma is roughly proportional to a 
summation of the fracture areas of the earthquake sources 
(Keilis-Borok and Malinovskaya, 1964). Pattern Sigma is 
identified as a peak in this summation. 

In algorithmic form the three patterns have been 
applied with some success to many regions. Pattern Sigma was 
first observed for 20 large earthquakes in the Central Asia, 
Eastern Mediteranean, Pamir and Assam regions (Keilis-Borok 
and Malinovskaya, 1964). Gasperini et al. (1978) found that 
an overwhelming majority of strong earthquakes in Italy were 
preceded by anomalous "Swarms" of weak events within 6 years 
of the mainshock. Keilis-Borok and Rotvain (1979) show that 
pattern S precedes strong (M27.2) earthquakes in Anatolia 
and surrounding regions. Sauber and Talwani (1980) apply 
pattern S to small earthquakes in the Lake Jocassee, South 
Carolina region. They show that the majority of M22.0 
earthquakes are preceded by the pattern; a precursory time 
of 10 days was used. Keilis-Borok et al. (1980b) apply 
patterns B, S, and Sigma to retrospective prediction of 


strong (M26.5) earthquakes in Southern California. The 
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patterns were found to precede 5 out of 6 strong earthquakes 
in the region. Patterns Sigma and B were found to precede 
strong earthquakes in New Zealand; pattern B was reported 
for Italy. A precursory interval of 3 years was used in the 
Southern California, and New Zealand studies. Keilis-Borok 
etmalmni980a)afoundethat™ 18)o0fm23astrongeearthquakess in 
five regions worldwide, were preceded by pattern B with a 
precursory time of 3 years. Keilis-Borok et al. (1980c) 
apply patterns Sigma and S to strong (M28) earthquakes in 
Tibet and the Himalayas. 

In the above studies many parameters used in the 
definitions of the bursts of seismicity patterns were 
determined a posteriori. Some examples of free parameters 
which may or may not have been chosen a posteriori include: 
minimum magnitude of a strong earthquake in a region; 
precursory time interval; time windows for defining 
anomalous activity; boundaries of the region; definition of 
Spatial size for the swarms of pattern S; minimum thresholds 
above which activity is considered anomalous; minimum 
magnitude considered for each pattern; and definition of 
aftershocks. The large number of free parameters involved 
constitutes the largest single difficulty in the evaluation 
of the significance of the above studies (Keilis-Borok et 
al., 1980a). When attempting to construct a hypothesis after 
the fact, the danger of self deception is considerable. By 
constructing a hypothesis with enough free parameters it may 


always be possible to fit the hypothesis to a known data set 
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even if the hypothesis is incorrect. 

Bursts of seismicity patterns are regional in nature, 
unlike foreshocks, quiescence, and the swarms of McNally and 
Evison. As such, they do not indicate the location of the 
future earthquake. In most cases the patterns take place at 
some distance (hundreds of kilometers) from the mainshocks 
which they precede. Accordingly, the occurrence of patterns 
B, S and Sigma may be interpreted as indicating a 
Significant increase in the probability of strong earthquake 
somewhere in the region within the next several years 


(Keil us-Bonokmetaaley, 1980b). 


1.7 Quiescence and Activation 

Many observed seismicity patterns can be interpreted as 
superpositions of quiescence and activation. The pattern of 
a precursory swarms followed by quiescence in the epicentral 
region suggested by Evison (1977a,1977b) is an example. 

Quiescence in the epicentral region accompanied by 
increased activity in surrounding regions has been termed 
the "doughnut pattern". Mogi (1969a) reported this pattern 
LOL ee ere large earthquakes in Japan. Yamashina and Inoue 
(1979) report a similar pattern before the magnitude 6.1 
1978 Shimane earthquake. 

An algorithm uSing a pattern Similar to the doughnut 
pattern is given by Keilis-Borok and Rotvain (1979). This 


premonitory pattern is termed A-Q (activation-quiescence). 
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Pattern A-Q is defined as a combination of activation and 
quiescence on opposite sides of an earthquake prone site 
which is in a state of quiescence. The pattern presumes that 
the scheme of lineaments and sites where strong earthquakes 
are possible are known for the region of interest. Examples 
of regional lineaments and the determination of sites where 
strong earthquakes are possible are given by Gelfand et al. 
(97 6 mpereGamkrfornia MGelfandset alew( 1974 ).for “several 
regions in Central Asia, and Southeastern Europe; and by 
GGuShkOvec teal mm lo.) efOresltaly, 

Keilis-Borok and Rotvain (1979) show that for Anatolia 
and adjacent regions seven of nine M27.2 earthquakes which 
occurred from 1900 to 1976 were preceded by pattern A-Q with 
an average precursory time of 13 years. There were, however, 
many false alarms. Joint use of this pattern and pattern § 
("Swarms") gave better results (fewer false alarms) with a 
shorter average precursory time of 6 years. 

Keilis-Borok and Rotvain (1978) suggest that pattern 
A-Q may be useful in the determination of soon-to-break 
seismic gaps. 

Talwani (1979) gives an empirial earthquake prediction 
model for small earthquakes ( 2.0 $ M S 2.6) in the Lake 
Jocassee area. He finds that premonitory seismicity is 
characterized by two distinct phases. The first phase, a 
phase, is characterized by a slow or no increase in overall 
seismicity and by quiescence in the epicentral area. The 


second phase, B phase, is characterized by an increase in 
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seismicity with a clustering of events around a "target" 
area defining the site of the future earthquake. The £ phase 
is similar to the Doughnut pattern defined above. The 
mainshock occurred after a period of quiescence in the B 
phase. A Similar pattern was observed before several large 
earthquakes in the Middle America region (Ohtake et al., 
1977). The similarity between the micro-seismicity before 
small events at Lake Jocassee and reported precursory 


patterns for large events is remarkable. 


1.8 Migration 

Many investigators have noted an almost linear 
progression in time of large earthquakes along various 
seismic zones. An apparent northeasterly migration of large 
earthquakes through Northeastern China at a rate of 110km/yr 
was instrumental in the successful prediction of the 1975 
M=7.3 Haicheng, China earthquake (Scholz, 1977). Other 
examples are numerous. Several authors have noted east to 
west migration at a rate of greater than 50km/yr along the 
North Anatolian fault of Northern Turkey (Mogi, 1968c; 
eniste ce 1970s Toksoz et) al. sio7 Q)melOkSOzmet salemmaio gD 
report a slower (<10km/yr) easterly migration for the same 
zone. Mogi (1968c) observes a southeasterly migration 
pattern towards the 1933 Sanriku, Japan earthquake at a rate 
of 150km/yr. After the occurrence of the Sanriku earthquake 


the migration pattern changed direction towards the 
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southwest and continued at about the same rate away from the 
mainshock. Mogi (1968b,c) suggests a worldwide pattern of 
migration for large earthquakes (M28.2) in three branches 
between 1900 and 1964. In the Aleutains the apparent 
migration of two sequences of large earthquakes (M27.7) at a 
rate of around i00km/yr are noted by Kelleher (1970). Wood 
and Allen (1973) observe an apparent migration of 
earthquakes above magnitude 5 in Central California at a 
rate of 3.3km/yr. Bufe et al. (1974) observe an apparent 
migration of micro-earthquakes to greater depths before 
moderate earthquakes in Central California. Engdahl and 
Kisslinger (1977) find that the foci of foreshocks 
precursory to a magnitude 5 earthquake in the Central 
Aleutians gradually migrated towards the future mainshock 
epicenter. Using spectral analysis techniques Delsemme and 
Smith (1979) observe an apparent migration of large 
earthquakes (M27.7) in South America at a rate of 95km/yr 
BEOnMmSoOUCHmLOsMiOG Lh 

The existence of a migration pattern is difficult to 
prove. Migration patterns are uSually based solely on visual 
observations and are therefore not amenable to statistical 
tests. An example of a quantitative statistical treatment of 
earthquake space time correlations iS given by Kagan and 
Knopoff (1976, 1978). Using the second-order moment ("the 
number of events in different magnitude ranges that occur in 
a certain time-distance interval from a given shock") they 


demonstrate that a weak but statistically significant 
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migration pattern exists for large shallow earthquakes with 


migration velocities in the range of 300 to 2000km/yr. 


1.9 Mainshock magnitude and. precursory time 

An apparent exponential dependency of precursor time on 
mainshock magnitude for some types of seismic precursors has 
been noted by several authors. Keilis-Borok and Malinovskaya 
(1964) find that the time difference between anomalous peaks 
in pattern Sigma and the subsequent mainshock increased with 
mainshock magnitude. A logarithmic dependence of precursor 
time on mainshock magnitude for several different types of 
long term precursory phenomena to moderate and strong 
earthquakes is reported by Scholz et al. (1973), Whitcomb et 
al. (1973), Evison.(1977a,b), and Kristy and Simpson (1980). 
Talwani (1979) shows a logarithmic dependence of precursory 
time and mainshock magnitude for small (2.0<M2>2.6) 
earthquakes in the Lake Jocassee region. 

Several authors have reported apparent relationships 
between the maximum magnitude observed in the precursory 
seismicity and the magnitude of the mainshock. Papazachos 
C1974, 01975 mandi Wu eteal .o(1976)esudgest thatethercmisea 
direct relationship between the magnitude of the mainshock 
and that of the largest foreshock. However, Jones and Molnar 
MITTS) HEindGnossuchirelationship-ebvasoneGiSe7arb) efindswa 
linear relationship between the maximum magnitude of a 


precursory Swarm and the magnitude of the mainshock. Kristy 
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and Simpson (1980) find a similar result for bursts of 
activity preceding large aerehquakes in Central Asia. 
Rikitake (1975b, 1979) shows that reported earthquake 
precursors can be grouped into separate categories on the 
basis of precursory time-mainshock magnitude dependence. 
Precursors ofthe first kind are characterize by a direct 
dependency of precursor time on the magnitude of future 
mainshock. This category includes most types of anomalous 
seismicity patterns with the exception of foreshocks. For 
Short term precursors, including obvious foreshocks, no 
Systematic correlation between precursor time and mainshock 
magnitude was found. A single exponential relationship 
between precursory time and mainshock magnitude seems to 
exustmt Onmal le precursors) of theltinstekind, although there 
is significant, scatter in the data (Rikitake, 1975b, 1979). 
Other factors such as the regional strain rate may be 


important. 


1.10 Criticism of previous works 

The previous discussion has given examples of reported 
SOUS patterns. Most of the patterns are based solely 
on visual observations of epicenter maps or space-time 
diagrams of often poorly located events, as such the 
patterns are not well defined. A common procedure followed 
in looking for the patterns is to simply look at the 


seismicity in the immediate space-time vicinity of a large 
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event and to note any possible irregularities which seem to 
precede the large event. Different investigators have 
reported different and sometimes contradictory patterns for 
the same event. Little statistical significance can be given 
to patterns developed from observations on single events. 
Examples of more systematic studies are given by 
Keilis=Borok et al. (1980a,b,c), Keilis-Borok and Rotvain 
(1979), Sauber and Talwani (1980), Talwani (1979) and 
Kanamori (1981). Such systematic studies assume that there 
are systematic precursory patterns which do not vary 
substantially from event to event and region to region. This 
assumption may not be true. 

What is needed in precurSory studies is a quantitative 
definition for the patterns. For example, in the precursory 
patterns swarms (for example, Evison 1977a) and quiescence 
(for example, Mogi 1979) there is no quantitative definition 
of what a constitutes a precursory Swarms or precursory 
quiescence. No studies are made to test for the occurrences 
of false alarms (i.e. the occurrences of "anomalous" swarms 
and periods of quiescence which are not followed by large 
events). Such studies need to be made before any statistical 
Significance can be assigned to many of the reported 
patterns. Also needed is some physical model which explains 
the observed patterns. 

In spite of the previous objections to the statistical 
Significance of some observed precursory patterns, it should 


be noted that precursory patterns (in particular precursory 
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quiescence and precursory swarms) have been studied 
extensively by many independent investigators. Occurrences 
of precursory patterns are well documented, any model of 
sequences of earthquakes which fails to explain the 


observations suffers a severe deficiency. 


1.11 On models of earthquake sequences 

All observed precursory seismicity patterns even 
quiescence can be thought of as manifestations of clustering 
of earthquake sequences. Little is known about the physical 


Nature, of this clustering. 


1.11.1 Stochastic models of clustering 

In the literature many stochastic models for earthquake 
occurrence have been proposed which result in clustering of 
earthquakes in time. Such models are usually used to explain 
the statisics of aftershock sequences. A review of several 
such models can be found in Utsu (1972) and Vere-Jones 
(1975). An example is a time-dependent Poisson process in 
which the mean number of events per unit time is a function 
of time. 

Another class of stochastic models with inherent 
clustering are the branching process models. Branching 
models start with one or a series of primary "events". These 
primary events generate secondary events with a probability 


that depends on the elapsed time from the primary event. 
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Secondary events generate subsequent event SaandesoOnOll, 
Branching models have been applied to models of crack 
expansion (Vere-Jones, 1976) as well as to statistical 
models of earthquake sequences (for example, Lomnitz, 1974; 
Kagan and Knopoff, 1981). The model of Kagan and Knopoff 
(1981) displays several interesting features and will be 
discussed below. 

In the Kagan-Knopoff model the probability that one 
evel aqivecmol Eth ietOmanothereinatnestimesinterval dt is 


given by 


P(t)dt=o, ae 


aS = We gare 


The singularity at time zero is eliminated by assuming that 
a parent earthquake can not produce other shocks in time t, 
where t, is the assumed rupture time of an elementary 
earthquake. They suggest that a value of B equal to .5 gives 
results statistically similar to shallow earthquake 
sequences. 

Using plausible assumptions for what would constitute 
an observed earthquake in such a model (ie the superposition 
of the effects of many elementary shocks which are spaced 
close in time), they show, through the statistical analysis 
of earthquake sequences generated by the model, that several 
time-magnitude properties similar to those observed in real 


Catalogues are reproduced. These include the 





eZ 


frequency-magnitude law and the Omori law of the rate of 
aftershock production. 

If such a model of earthquake occurrence is correct 
then present shocks can be thought of as a being the 
"children" of all previous earthquakes, or in a sense any 
earthquake is an "aftershock" of some earthquake occurring 
before it. All or most earthquakes will be interdependent. 
This presents a conceptual problem in dividing earthquakes 
into foreshock - mainshock - aftershock sequences if a 
purely statistical definition is assumed in which mainshocks 
are defined simply as those events which are statistically 
independent of other events and aftershocks and foreshocks 
are defined as statistically dependent events. This point 
will be discussed later. There is a related problem. We have 
only a short catalogue. In discussing statistical properties 
of the catalogues we are limited to time scales that are 
much less than the length of the catalogue, except for 


speculative purposes. 


1.11.1.1 Deficiencies of stochastic models 

A stochastic model of earthquake occurrence on the 
aroha ta setS available can not be demonstrated to be 
correct. All that can be shown is that the model agrees with 
some observed properties of the sequences which are also 
extracted from short data sets. It should be remembered that 
the statisical properties of earthquake sequences are all 


based on a short data set. The short data set does not allow 
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for the description of longer term properties of the 
Sequence. Omori's law postulates one such long term 
behavior. It should be remembered that Omori's law is based 
on a relatively short data set, there is no guarantee that 
it is valid for very long times. There is great ambiguity in 
what constitutes an aftershock; this ambiguity increases 
with the time from the mainshock. Any theory which rests on 
the vague notion of what is or is not a long term aftershock 
rests on a very weak foundation. 

Stochastic models of earthquake occurrence ignore one 
very important aspect of the earthquake problem, namely the 
mechanical properties of the system. Little physical input 
is used to constrain the free parameters inherent to such 
models. There is no physical input in such models which 
accounts for the causes of seismicity (the energy source) or 
which takes into account the physical properties of the 
fault zone. Nothing is said about the fundamental cause of 
seismicity, namely the release of strain which accumulates 
due to plate motions. Stochastic models of earthquake 
sequences do not, nor can they, explain the observed 
precursory patterns. What 1S needed iS some model which also 
incorporates the physical properties of the failure process 


and the mechanical properties of fault zones. 
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1.11.2 On the physical basis of clustering 

Several qualitative physical explanations for 
earthquake clustering have been proposed. These include: 
inhomogeneities in the material strengths along a fault zone 


and non-linear dynamic friction on the fault zone. 


1.11.2.1 Non-linear dynamic friction 

Using a simplified two-dimensional model of the 
expansion of a crack Barenblatt et al. (1981) show that, 
with an appropriate assumption for the functional form of 
Enegdynamucmirretvoneasvastunction Of=slip rate andthe 
modulus of cohesion as a function of the rate of expansion 
Olb@thescrack,~expansioneofethe crack is@realized iby 
alternating slow and fast phases. The number of alternations 
becomes large in the presence of local maxima in the shear 
stressvor *localiminima in the friction “along*the ‘fault. They 
suggest that the effect may be a possible explanation for 
the occurrence of premonitory clusters at large distances 
relative to subsequent strong earthquakes. 

The® functional form of the dynamic friction as a 
function of slip rate has the following characteristics: for 
ube lis rate the dynamic friction increases with 
increasing slip rate; for intermediate values of slip rate 
the dynamic friction decreases with increasing slip rate to 
some minimum; after the minimum is achieved dynamic friction 
increases steadily with increasing slip rate. The modulus of 


cohesion: 1S assumed to increase with rate of crack 
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expansion. There is little data supporting the assumed 
functional form of the dynamic friction as a function of 


slip rate. 


1.11.2.2 The barrier or asperity model 

Irregularities in the physical or mechanical properties 
along a fault zone provide another possible explanation for 
clustering. Regions of high strength along a fault zone are 
usually termed asperities Or equivalently barriers. The 
concept of asperities has been used to qualitatively explain 
many features of seismicity. These include, for example: 
complex events (Das and Aki, 1977; Aki, 1979: Kanamori and 
Stewart, 1978; Lay and Kanamori, 1980; and others); 
non-uniform seismicity along fault zones (Wesson and 
Ellsworth, 1973); the dependence of the frequency of 
foreshocks on the time to the mainshock (Jones and Molnar, 
1979); and seismic clustering (Ishida and Kanamori, 1978, 
EEO. 

Kanamori (1981) shows that a simple asperity model can 
explain many features of the various observed precursory 
Se1smicity patterns. The followingewill) give a discussiconyot 
this mode Since it has great relevance to this work. 

The model assumes that a fault is composed of a series 
of subfaults with varying strengths. The strength of a 
subfault is defined as the maximum loading stress it can 
acommodate before it fails. A region which is characterized 


by the occurrence of large earthquakes is a zone which is 
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characterized by an increase in the strength of subfaults. 
Accordingly, an asperity is introduced as a region of high 
strength of subfaults’relative to the "normal" distribution 
of strengths. As a simple example Kanamori (1981) considers 
the case with one asSperity. The strengths on the subfaults 
are assumed to follow a gaussian distribution with mean S 
and standard deviation a. The strengths of segments on the 
asperity are assumed to follow a similar distribution with 
larger mean (S,) and with different standard deviation (a,). 
The model further assumes that a section of the fault fails 
when the stress on that section exceeds the strength and 
once a section fails the loading stress it accommodated is 
held uniformly by the remaining unfailed sections. Once a 
subfault fails the stress there drops to zero and no healing 
occurs. In such a model the rate of stress increase on 
unbroken sections accelerates as the number of broken 
sections increases. 

Consider qualitatively what happens when the loading 
stress increases linearly in time in such a situation. When 
the tectonic loading stress is low compared to S only a 
small numbers of subfaults will fail as a scattered 
activity. As the loading stress increases approaching S the 
process accelerates until most subfaults outside the 
asperity are broken. The increaSed activity at this stage 
could result in swarm-like activity since the loading stress 


will exceed S rapidly. 























sPIANFEUS: FG HAD seonsegt! as gh Gee 


eet) 't ror“? 6 ' :; & 4 ve) ne ae 


Zrabienan ¢{ ERO ; : . Sdspeetze 36 
/ 7 
is F a4 rises. ot 1% @869 OAT 


t - 
& sens Pe o38 7 7 


‘ohvnsée Sag 
F _ 


4 veo Peeggays 


ee IEOweTe 
7 


» BUC onT . 


oh) non, 


nd 


¢ a] i? iio : 


sii pein. ia “4k seeuee 
Pe 

7 di 

. ih Eine vs . ‘bojoes HS 
7 7 
—_ Ose 7 ' ; : ¥* oeipe rz ' 





te ren ees HO2997501 €aec NiVoel aif e& aivids ft 


: y 





—% ' it sPiasne Eine n ipsa rean bh tity ovinseteree 
. a bere betas 


7 


ss 


3% 


After this stage quiescence could result since the 
sections of the fault which remain are relatively strong. 
The stress will be concentrated on the asperity, the area 
Surrounding the asperity will be essentially decoupled. This 
may also result in loading of adjacent fault zones which 
could produce a doughnut pattern in the seismicity. In the 
final stage as the loading stress approaches the strength of 
the asperity S, the asperity itself will begin to break in 
an accelerated process possibly resulting in a concentrated 
Swarm which terminates with the mainshock. Kanamori (1981) 
gives a simple numerical simulation which reproduces the 
above behavior. The time length of the periods of quiescence 
and of the precursory swarm which terminates with the large 
everitmalescOontnolled by the parameters S,/S fandsa,.- AS S,/S 
decreases the length of the quiescence decreases. AS a, 
increases the length of the precursory swarm which 
terminates with the large event increases. If a, is too 
small no such precursory Swarm appears. 

Such a process, although simplistic, gives a possible 
qualitative explanation of swarms, quiescence, and 
activation before large earthquakes. It does not explain the 
large distances and times which are sometimes attributed to 
precursors. 

Furthermore, the model suffers from too many free 
parameters. In such a qualitative picture it might be 
possible to "model" any sort of activity by varying the 


number and characteristics of the asSperities. Nothing is 
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Said about the nature of the interaction between different 
fault zones. No healing mechanism is present. Such a model 
can not explain the apparent recurrence of large earthquakes 
at nearly the same location without some sort of healing 
mechanism for the asperities. However, the model is perhaps 
useful as the simplest possible model which could explain 


the complexity of observed seismicity patterns. 


1.12 Concluding remarks 

Anomalous seismicity has been observed before a large 
number of earthquakes in nearly all magnitude ranges. Most 
patterns are based solely on visual observations of 
epicenter maps or space-time diagrams. Such visual patterns 
are only qualitatively defined. Some examples of attempts at 
Quantitative descriptions of seismicity patterns are given 
by Keilis-Borok et al. (1980a,b,c), Keilis-Borok and Rotvain 
(1979), Gasperini et al. (1978), Sauber and Talwani (1980), 
and Talwani (1979). Different investigators have reported 
different and sometimes contradictory patterns for the same 
event. Quiescence seems to be the most commonly observed 
petrerae 

In a global survey of seismicity patterns before major 
earthquakes along subduction zones Kanamori (1981) finds 
Significant regional variations in the pattern of which 


areas have foreshocks, quiescence, and premonitory Swarms. 
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Although there are significant variations in the 
premonitory patterns for individual events, on the average 
different magnitude ranges show remarkable similarities if 
appropriate time and space scaling is taken into account. 
This may be another manifestation of the self similarity of 
the earthquake process proposed by Kagan and Knopoff (1980). 
If so, the study of micro-seismicity preceding relatively 
numerous small events may be useful in determining the 
nature of seismicity before large earthquakes. 

Many patterns may not be sufficiently obvious to be 
observed before the occurrence of the mainshock. Also, it is 
not known how many false alarms would result from the use of 
the various patterns (for example, there have been no 
studies on the numbers of anomalous swarms and periods of 
quiescence which were not followed by large earthquakes). 
This along with variations in patterns from event to event 
limits the usefulness of seismicity patterns alone in 
earthquake prediction. However, with other constraints such 
as estimates of recurrence intervals and locations of 
possible strong earthquakes, seismicity patterns may be 
useful in forecasting the times, places, and perhaps 
magnitudes of future earthquakes at least on a long-term 


basis. 
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2. A scaling law for the occurrence of aftershocks in 


Southern California 


2.1 On self similarity 

Mhew concep or stochastic self-similarity is a familiar 
and important concept in the field of turbulent flow. The 
Scliacimiiatityeotmturbulent flows implies that the "picture" 
of turbulence shows no characteristic scale. The picture of 
flow will look statisically similar at any magnification 
between an outer scale, determined by the driving force, and 
an inner scale, where viscosity becomes important. Recent 
WOGKSmOteKagdaleandehnopotia(1976,.1977, 1978" 1980a. 1980b) 
Suggest that seismic process shows similar stochastic 
self-similarity, (i.e. there seems to be an absence of any 
particular scale connected with time - distance - magnitude 
patterns of earthquake occurrence; epicentre maps have the 
Same appearance regardless of scale). Kagan and Knopoff 
(1978, 1980b) report that this self-similarity holds over 
magnitude ranges extending from M=1.5 to the largest 
earthquakes. Insufficient data 1S available to test the 
Mepoes eee WK igen 

Andrews (1980) shows that the Gutenberg Richter 
frequency magnitude relation of earthquake occurrence 
results by assuming self-similarity alone. Furthermore, 
self-similarity manifests itself in the form of known power 


law distributions of many features of seismicity. These 
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include: the power law distribution of seismic energies and 
Omori's law for the rate of occurrence fi aftershocks and 
foreshocks of shallow events (Kagan and Knopoff, 1978); the 
power-law dependence of the energies of foreshocks and 
aftershocks (Kagan and Knopoff, 1978); and the inverse 
power-law dependence of the spatial moment on the separation 
between two foci (Kagan and Knopoff, 1980a). 

The apparent self-similarity of the failure process has 
great relevance to the identification of anomalous 
seismicity patterns. In a given region large events are 
rare. The time spans of earthquake data sets are usually 
short compared to the return times of the largest events in 
a region. This means that the statistical basis of any 
phenomenon associated with large magnitude earthquakes will 
at best be very weak. Self-similarity implies that the study 
of anomalous seismicity in the occurrence of small events 
before relatively numerous moderate events may be useful in 
the identification of anomalous seismicity in the sequence 
of the moderate events themselves, provided that appropriate 
time and space scaling is taken into account. 

This chapter explores an aspect of the self-similarity 
of the earthquake process, a Scaling law for the occurrence 
of aftershocks in Southern California. The research which 
forms the basis of this chapter was carried out at the 
Seismological Laboratory, California Institute of 
Technology. The chapter is based on the paper "A scaling law 


for the occurrence of aftershocks in Southern California" by 
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Reva moOueauxe sel uhelliS—-BOorok, and ke, H.. Hutton which, has 
been submitted to Physics of the Earth and Planetary 


Interiors. 


Zc ntroduction 

The occurrence of aftershocks has been intensively 
Studied in many branches of seismology. These include: 
Statistical models of the earthquake source (Utsu, 1972; 
Vere sJvOnes, mic Gsmnagan and Knopott-.11976)-- mechanical 
models of the earthquake source (Aki, 1979; Kanamori, 1981; 
Barenblatt et al., 1981): and estimation of seismic risk and 
earthquake prediction. Some of the premonitory seismicity 
patterns discussed in Chapter 1 such as seismic gaps (Sykes, 
1971) and "bursts of aftershocks" (Keilis-Borok et al., 
1980a, 1980b) are directly based on the observed pattern of 
aftershock occurrence. | 

The above mentioned studies require a knowledge of 
Statistical properties of aftershock sequences. Not many 
such properties are known. Utsu (1972) gives a review. One 
of the few comparatively well established regularities is 


the Omori Law 


n(T(e)) =C/[Tle) +B] 7 


where n(T(e)) is the average number of aftershocks during 
theseth day after main «shock, C) avand) Bare constants. The 


coefficient C obviously depends on the magnitude of the 
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mainshock and the magnitude range in which aftershocks are 
counted. It also depends on the location of the hypocenter 
and probably on the mechanism of the main shock. However, 
even when all the above factors are similar, n(T(e)) shows 
strong variations from event to event. Omori's law only has 
meaning aS a rough average for many events. As such, it is 
valid for different formal defintions of aftershocks 
(Gardner and Knopoff, 1974; Vere-Jones, 1976; Keilis-Borok 
et al., 1971). The commonly accepted values for a and Bf are: 
a=1 and B=0. The large number of factors which influence 
n(T(e)) presents a great difficulty in the estimation of C, 
a and B. There are probably too many parameters for the 
available data. 

This chapter will test the hypothesis that the number 
of aftershocks does not depend on the magnitude Mm of the 
mainshocks if aftershocks are counted in a magnitude 
ineervaletrom (Mm-A) to Mm, where A iS sa constant. The total 
number of aftershocks over a period of the order of a year 
is considered. In terms of Omori's law, the hypothesis means 
that for a fixed A, C does not depend on Mm. This 
hypothesis, if correct, provides a self-similarity which 
imposes an additional constraint on the study of the 
physical nature of aftershocks. Also, the elimination of one 
parameter (Mm) reduces the number of observations necessary 
for statistical studies of aftershock sequences. 

In the next section the above mentioned scaling law is 
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catalogue 1932 through 1980 (Hileman et al., 1973; Friedman 
et al., 1976; Subsequent Preliminary Epicenter Listings) 

provided the data for the study. Finally the scaling law is 
used for the analysis of a long-term premonitory seismicity 
pattern termed "bursts of aftershocks" (Keilis-Borok et al., 
1980a, 1980b). This pattern is also discussed in Chapters 1 
and 4, Keilis-Borok and Prozoroff (1981) give a similar 


study for the strongest earthquakes worldwide (M 2 7). 


2.3 The scaling law 

To test the hypothesis formulated in the introduction, 
the number of aftershocks for different mainshock magnitude 
Mm was counted for several values of A. The following simple 
definition for an aftershock was used: Consider two 
earthquakes with time Sequence numbers i and j, i>j. The 
second earthquake is an aftershock of the first if the 
following conditions are satisfied: The distance between 
their epicenters is less than R(Mi); the time difference 
a esl M 1) ed nOeM js <0Mi.e T(M)eander(M)manesempunicad 
functions. There is no physical basis for this definition. 
The aamarete home T(M) and R(M) can be viewed as simple box car 
windows in which obvious "mainshock-aftershock" clustering 
is evident. The thresholds R(M) and T(M), assumed after 
Gardner and Knopoff (1974) are given in Table 2.1. They seem 
to give an estimation from above (i.e. for most mainshocks 


R(M) and T(M) are larger than necessary). An exception is 
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Windows for the identification of aftershocks 
Wet ter eGardner tand Knopot fee 19741)" 


M R (Km) T(days) 
25 toa 6an0 
CaO 2 es Ess.) 
SEO 2650 2270 
AO) 10h, LO) 
ah S50 8350 
Sel!) aa (8 155:0 
5) AI. 230.0 
6.0 S40 ori) 48 
65 Ga, 7'S0.0 
vis) TORO Satan 16 
Fis Bre 0 960.0 
Sr 94.0 985.0 


Table 2.1 gives the empirical windows for the identification 
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the Kern County earthquake, July 1952, M=7.2. For this 
earthquake R(M) and T(M) are greater than those values 
indicated in table 2.1. A map of mainshocks with magnitudes 
greater than or equal to 4.0 for the years 1932-1980 using 
Bhesparamerenssor table 2:4 ais givenin figures 2.1. 

The total number of aftershocks (bi) were considered. 


It follows that 

T(M) 

bj=). n(T(e)) 8 

e=1 
where n(T(e)) is given by (1). The results of the 
Computaulonssare shown in tables) 2.29-a224. It ist presumed 
that the Southern California catalogue is reasonably 
Complete slOrsMucw 3. Sefromethe begdinning, (1932) .aThis allowed 
the computation of bi for A=1, Mm24.5 for mainshocks from 
1932, For the years 1970 - 1980 earthquakes with magnitudes 
Greater than or equal to 3.0 were considered and for 1977 - 
1980 events with M greater than or equal to 2.5 were 
considered. These divisions are approximate. 

The above places obvious constraints on the miminum 
mainshock magnitude which could be considered for the 
differént values of A. The number of smaller magnitude 
mainshocks occurring after 1977 was insufficient to test the 
hypothesis for A=2.0. To test the hypothesis for sa=2.0 the 
above constraint was relaxed; M22 enon 1970 and Mz22..08from 
1977 were included. The results for sa=2.0 could be biased 
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Figure 2.1 gives an epicentral map of mainshocks occurring 
in Southern California from 1932 through 1980 as determined 
by the aftershock criteria of table 2.1. The approximate 
CONLOCUDMOLA them vVoOlcanicuzone afttem(Allen teal.) 1965) 15) 
shaded. 
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Frequency distribution of numbers of aftershocks 
per mainshock, delta = 1.0 
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Taplow 2 yeni seattabulation oft thes: requencypon occUUumencesrOr 
the number of aftershocks per mainshock with A = 1. The 
mainshock magnitude range 4.0 to 4.4 uses data after 1970, 
other mainshock magnitude ranges use data from 1932. 
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Frequency distribution of numbers of aftershocks 
per mainshock, delta = 1.5 
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tables2 3 'smantabulation of the frequency of occurrence for 
the number of aftershocks per mainshock with 4 = 125 sethe 
Mainshock magnitude range 4.0 to 4.4 uses data cet Clemo 
4.5 to 4.9 data after 1970, other mainshock magnitude ranges 
USerddtamrrome 1937. 
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Frequency distribution of numbers of aftershocks 
DeGemainshochemde lta s=—2.0 
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Table 2.4 is a tabulation of the frequency of occurrence for 
the number of aftershocks per mainshock with a = 2. The 
mainshock magnitude range 4.0 to 4.4 uses data after 1977, 
4.5 to 5.4 data after 1970, other mainshock magnitude ranges 
LSemCaitasen tl. OMnesl, os 2. 
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Tablecmenc wee swand. 2.4 Show trequency distributions of 
bevVersusemainchnocksmagnatude for A=(.0) 1.5, and 2.0 
respectively. Figures 2.2 -2.4 show the corresponding 
normalized cumulative frequency distributions b(q), where q 
is the number of mainshocks with b or less aftershocks. The 
hypothesis predicts that the frequencies tabulated in tables 
2.2 - 2.4 should be independent of mainshock magnitude 
interval. To test this each table was divided into magnitude 
intervals such that each interval contained similar numbers 
of events. The frequencies in the magnitude intervals were 
then compared by the Chi-squared criteria. The results are 
shown in table 2.5. The hypothesis that the frequencies are 
the same cannot be rejected at the ..28 level of significance 
HOT hegworst. case | (A=1.0) . 

bavtesme 26>) 2, 5and figures 2.2 — 274 seem to confirm 
GUmsny pothesis sat least) for Asi. 5vand) Mn=4.0. The catalogue 
1s insufficiently complete to make conclusions for smaller 
Mm or larger A. For A=2.0 the results are encouraging, 
although for later years the minimum aftershock magnitude 
used in this case is probably not justified by catalogue 
completeness. The cumulative frequency distributions for 


different A are juxtaposed on figure 2.5. 
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Figure 2.2 gives the normalized cumulative frequency 
distribution b(q)" A = 1 for various magnitude ranges. 
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Figure 2.3 gives the normalized cumulative frequency 
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Figure 2.4 gives the normalized cumulative frequency 
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Chi-square comparisons of the frequency distributions 


Key delta X q p 

(A) 10 AOS 2 (he 
(B) ies5 5 real) 3 .04 
To) 220 PV TE 3 Ou 
(D) #25 20 3 .94 
lasy, iG (es 1 sis) 


(A) gives a comparison between the magnitude range 4.0-4.9 
for the years 1970-1980 and the magnitude range >4.9 for 
for the years 1932-1980. 

(B) compares the magnitude ranges 4.0-4.4 after 1977 plus 
4.5-4.9 after 1970 to the ranges >4.9 after 1932. 

(C) compares the magnitude range 4.0-4.4 after 19 7V/aons 
4.5-5.4 after 1970 to the combined distribution for 
all magnitudes greater than 5.4 occurring from 1932. 

(D) compares the combined distribution of the magnitude ranges 
4.0-454 (1977-1980), 4.5-4.9 (1970-1980) and >4.9 
(1932-1980) which occurred in the volcanic region 
outlined in figure 2.1 to the similar Gistimiow ti one ror 
the non-volcanic region. 

(E) compares the frequency distributions of ear thquakes 
found to be foreshocks to the expected distribution 
for mainshocks. The number of aftershocks was counted 
for two days. Earthquakes with magnitudes greater than 
4.4 were considered for al] years. Those with magnitudes 
4.0-4.4 were considered after 1970. 


a ee Se el ee eee 


q gives the number of degrees of freedom 

X gives the value of Chi-squared 

P is the Probability that a random variable which follows the 
Chi-square distribution with n degrees of freedom is less 
than or equal to x. 


Table 2.5 gives the Chi-square test of goodness of fit 
between various magnitude intervals from the frequency 
Gabe cml nae 
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Figure 2.5 shows a comparison of the normalized cumulative 
Prequency Cistrabutions for diiterentmvalucsmoleAs mine 
distributions are based on all mainshocks for which the 
catalogue was considered sufficiently complete (ie Mm 2 
Mate Newer emi nt) =>. Lisenioy Ome mOm romeo y0mtOlnLOd 7 
ance oma Leoni Pet hemor St oii onsmnormAac. OsandeA=205 
are based on small numbers of points (64 and 31 


respectively). 
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2.4 Regional variations 

Examination of the catalogue reveals that the 
elustering Ofsearthquakes 1s different in) different parts of 
Southern California. For example, the area of Quaternary 
volcanism is known as "aftershock prone" (Allen et al., 
1965). This area is characterized by more intensive 
GIUSEGrING;wine. clusters with larger numbers of events. 
Accordingly, b should be larger on the average in this 
region. 

Theme DpaOximaLemcOntOur ol themvyolcanic area latter 
Aleve teal elo )mLSeshownmin figures, 1. The cumulative 
PECGuenCVecE Ser MouiGronwOtaberonstnemvolcanic and 
non-volcanic areas with A=1.5 are given in figure 2.6. The 
differences between the two regions are evident from the 
figure. For a=1.5 the frequency distribution for the 
volcanic zone is characterized by fewer mainshocks with zero 
aftershocks (23% compared to 33% for the non-volcanic 
region) and more mainshocks with greater than 20 aftershocks 
(7% compared to 0% in the non-volcanic region). 

The frequency distributions for the two areas were 
compared by Chi-squared calculations. The results are shown 
idetables2.5.) Theshypothesis thatethnesdistributronsmare 


Similar can be rejected at the .06 level of significance. 
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STMBOLS 
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Figure 2.6 give a comparison of the normalized cumulative 
frequency distributions b(q) for the volcanic and the 
non-=volcanic, regions )(A=1/.5). 
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2.5 Premonitory pattern B 

According to Keilis-Borok et al. (1980a, 1980b) the 
Occurrence Of vam burst of aftershocks” (pattern B) indicates 
a Significant increase in the probability of a strong (M > 
M,) earthquake in the same region within the next t, years. 
"Pattern B" consists of a mainshock in the medium magnitude 
range with an anomalous number of aftershocks concentrated 
at the beginning of the aftershock sequence. Specifically 
the pattern is a mainshock with magnitude Mi such that M, < 
Mi < M, and bi(e) 2 B where bi(e) is the number of 
aftershocks in the first e days following the mainshock and 
B is a threshold above which bi(e) is considered to be 
anomalous. The following values were selected for Southern 
California by retrospective juxtaposition of strong 
earthquakes and *the values of bi: M,=6.5: M,=M,-:1; 
M,=M,-1.; t.=3 years; B =13; and e=2days (Keilis-Borok et 
al., 1980b). Aftershocks were counted in the magnitude range 
BEones tosMm. 

The conclusion of the previous section allows us to 
choose B independently of Mm and juxtaposition with strong 
earthquakes, although not unambigiously. We define B(q) as 
the q% quantile of the empirical distribution of the values 
of bi (i.e. B(q) is the value of b which is greater than the 
bi values of q% of the mainshocks). 

Completeness of the catalogue restricts the maximal 
value of A we can use. For Mm=5.5 we are restricted to s<2.0 


for early years. We considered the distribution of bi(2 
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days) for A=2.0, Mm>5.5, for the years 1932 - 1970, Mm>5.0 
for 1970 - 1977, and Mm>4.5 for 1977 - 1980 and obtained 
B(95) = 13 and B(90) = 11. 

Mainshocks with M>5.0 and bi(2 days)>11 are juxtaposed 
with strong earthquakes in table 2.6. As expected the 
results are very similar to those obtained by Keilis-Borok 
et al. (1980b). Most events found to be anomalous previously 
are found to be anomalous here. The exceptions are listed in 
table 2.7. Seven of nine strong events are preceded by 
anomalous bi. There are two false alarms, one in 1935 and 
one in 1958. The false alarm in 1935 is generated by long 
distance aftershocks of the 1934 M=6.5 and M=7.1 strong 
events. 

For smaller A the identification of anomalous events 
becomes more sensitive to the magnitude estimates, 
particularly to the mainshock magnitude. For example, 
Keilir s Bocokmeteal, \(l980b) Sidentit tedrthe Oct. 495 1978 
Mm=5.8 earthquake aS anomalous with bi(2 days)=35. For s=2.0 
bi(2 days)=9 which is not identified as anomalous at the 90% 
quantile level. However, if the mainshock magnitude were 
over estimated by .1 two more aftershocks should be included 
and bi=11 which would be identified as anomalous. An 
increase in A would make possible magnitude errors less 
Significant as well as increase the aftershock numbers. This 
would make the separation of anomalous events from 
non-anomalous events clearer. It is therefore suggested that 


as large a magnitude range as is possible should be used to 
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Mainshocks with anomalous bi (delta=2.0) 
and strong earthquakes 


Year day Ni lat ese, Veni) Ashi 
1933 BO SMO ome SO eD le if 7aeO jee GAs 613 
1934 BOOuay eel eo 5 eee 1 15 250 

1935 UUel  Biclt) = SSS Sapeiarer yi = 
1935 Abe TaN) Sey eS aed a, = 
1940 VA OMG (S27 3 ee 115. 50 

1944 era erilet meu eaters ala 18 
1942 2OCmeOn OM ISO Oye. 6o00 

1947 AM ey © fe Aten SVC PE ST ou oily. kay ZS) 
1948 SOI@SCE O igs. OShiea 16068 

1950 POST eot 4 ti 3S Obwianl 158 oven ete = 
1950 2 OPO OmsoS alee] 1 5e vee No iS 
1952 2OSMMe 2 i350 4 O0res8 192,02 

1954 SL OMRICl 3 Sie S0NGEi 60 0 t1'7 19 
1955 SMM 4am o sO0te lt Se 50a 14 = 
1956 AUR OF OMeO ARG Oma Oo! 6 

1958 SO AMIDE ES 2 Sra Se oe iM 
1969 BOmeor Cost 20m 114.20) .44 49 
1971 4 Oma 49004 4h E1840 

1979 ioc OMe OC erie owe 1 4 

1979 (Reo ee aioe wi oes A 

1979 PoOmsOmOUmNG © Olam i ops? 

1980 (Gis Tish) Sesion Oo ee 


bi was counted for a 2-day interval after the 
mainshocks. 
*bi indicates values of bi 
(1980b). 


TLOmaKel Sa bOLOK re tarad ae 


Table 2.6 lists mainshocks with anomalous bi for a = 2.0 
along with strong earthquakes M,26.4. Strong earthquakes are 


underlined. 
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Mainshocks with anomalous bi from Keilis-Borok 


et al. (1980b) not found to be anomalous for 
Oe ttam—aee 0 
Year day Mi lat Long Oh a 8y) 
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bi was counted for a 2-day interval after the 
mainshocks. 

*bi indicates values of bi from Keilis-Borok et al., 
(1980b). 
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identify anomalous bi. 


2.6 Aftershocks of foreshocks 

In the test of the scaling law we have eliminated 
events found to be foreshocks from our analysis. The 
definition of foreshock used here is similar to the 
previously mentioned definition for aftershocks with the 
exception that the condition Ma<Mm is reversed to Mf<Mm, 
where Mf is the magnitude of the foreshock. The same 
thresholds R(Mf) and T(Mf£) were assumed. 

The number of aftershocks bi(e) with A=1.0 was computed 
POrmeach fOresnocks fOr whichytm=tiee, A vaiue of 2 days was 
chosen for e. Figure 2.7 shows the corresponding normalized 
cumulative frequency distributions of bi(2 days) for 
foreshocks and mainshocks. The frequency distribution for 
foreshocks was compared to the expected frequency 
distribution generated by mainshocks using the Chi-squared 
criteria. The result is given in table 2.5. The hypothesis 
that the distribution for foreshocks is the same as that for 
mainshocks can be rejected at the .01 level of significance. 

racers 2.7 leaves the impression that on the average 
foreshocks are characterized by larger b. However, the total 
number of foreshocks, 27, is too small for definite 


conclusions. This difference deserves further attention. 
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Figure 2.7 shows the normalized cumulative frequency 
distributions b(q) for foreshocks and mainshocks with s=1.0. 
Aftershocks were counted for 2 days. 
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2.7 Summary 

The scaling law hypothesis seems to be Supported, with 
the qualification that apparently the distribution of 
aftershocks is regionally dependent. Regional dependence, if 
present, increases the amount of data which is necessary for 
Statistical studies of aftershock occurrence. The regional 
dependence may be related to changes in local tectonic 
Style. Areas characterized by predominately extensional 
tectonics may show systematic differences in aftershock 
occurrence from areas characterized by predominately by 
compressionial tectonics. 

It is suggested that, on the average, foreshocks 
generate larger numbers of aftershocks than do mainshocks. 
More data is needed before this observation can be tested. 

The scaling law presented here as well as 
self-similarity in general have several implications for 
premonitory pattern "bursts of aftershocks". Anomalous 
bursts of aftershocks should be identifiable even for 
smaller magnitude mainshocks, provided the catalogue used is 
complete enough. The regional dependence of the scaling law 
may complicate the identification of what constitutes 
anomalous activity at a particular location. Keilis-Borok et 
al, (1980b) suggest that in the entire Southern California 
region anomalous bursts of aftershocks in the mainshock 
Tange Of 1.5 to 6.4 aren precurnsony,toeM, 2) 6.5) earthquakes, 
With a precursory time of the order of 3 years. If this is 


correct the scaling law and self-similarity suggest a 
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further question: Are anomalous bursts of aftershocks with 
smaller mainshock magnitudes precursory to intermediate 
magnitude events on a shorter time-distance scale? The 
Southern California catalogue is probably insufficiently 


complete to test this hypothesis at present. 
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3. Swarms and clusters of mainshocks in Southern California 


"Figures often beguile me, particularly when I have the 
arranging of them myself; in which case the remark 
attributed to Disraeli would often apply with justice and 
force: 'There are three kinds of lies, lies, damn lies and 
Staristicsis 


Mark Twain 
Autobiography, 
Volume 1 page 246. 


Serene roducti on 

Clustering of earthquakes in space and time is one of 
the most prominent features of seismicity. Phenomenology and 
Statistical models of clustering are described in many 
Studies (i.e. Utsu,. 1972; Knopoff, 1971; Keilis-Borok et 
al., 1971; and’Dziewonsky and Prozorov,.1981). Many of the 
proposed precursory seismicity patterns described in Chapter 
1 involve recognition of abnormal clustering or 
anti-clustering of earthquakes. These include, for example: 
bursts of aftershocks (Keilis-Borok et al., 1980a,b); swarms 
(Caputo et al., 1977; McNalley, 1977; and Evison, 1977a,b):; 
doughnut patterns (Kanamori, 1981); and seismic gaps of the 
second kind (for example, Mogi, 1979). The physical 
mechanism of clustering remains unclear: in particular, it 
is difficult to explain the interconnections between 
earthquakes over such large space and time intervals, which 


are suggested by some of the observed patterns. In the above 
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mentioned patterns no distinction was made between 
mainshocks, foreshocks and aftershocks, all were included. 

Although there is no unambiguous physical definition 
for either a mainshock or an aftershock, the most obvious 
examples of earthquake clustering are mainshock - aftershock 
sequences. This work uses a very simple definition for 
mainshocks and aftershocks which will be discussed in the 
next section. The following questions arise: Do the 
sequences of obvious mainshock - aftershock clusters show a 
higher level of clustering? Are the "mainshocks", defined by 
the simple criteria, themselves clustered on a larger time - 
distance scale than is seen in the mainshock - aftershock 
clustering? If the mainshocks are clustered, are periods of 
abnormal clustering precursory to large earthquakes? Are 
precursory swarms and periods of quiescence observable in 
the sequence of the simple mainshock - aftershock clusters? 

Mhasechapter will discuss some of the statistical 
properties of the sequence of mainshocks in Southern 
California. The question of mainshock clustering is explored 
and, the clusters of mainshocks are described and juxtaposed 
with strong earthquakes. Part 1 uses a simple swarm model of 
clustering. In part 2 several clustering sequences are 
examined in more detail by means of cluster analysis. 

Much of the research on this chapter was carried out at 
the Seismological Laboratory, California Institute of 
Technology while the author was assisting Prof. V. I. 


Keilis-Borok during his tenure there as Fairchild Scholar. 
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Part 1 of this chapter is based in part on the paper "Swarms 
Cte Matnmohockhcmitecouthernucalitonunia by V. 1. 


Keilis-Borok, R. Lamoreaux, C. Johnson, and B. Minster. 


3.2 On the definition of mainshocks and aftershocks 

The question of the definition of an aftershock is a 
fundamental one. A conceptual definition of an aftershock is 
an event which iS a direct result of the redistribution of 
energy resulting during failure of a previous shock. 
Stochastic models of earthquake sequences such as the 
branching model of Kagan and Knopoff (1981) mentioned in 
Chapter 1 present a conceptual problem in dividing 
earthquakes into foreshock - mainshock - aftershock 
sequences if a purely Statistical definition is assumed in 
which mainshocks are defined simply as those events which 
are statistically independent of other events and 
aftershocks and foreshocks are defined as statistically 
dependent events. In such models all or most earthquakes 
will be interdependent. If one wishes to discuss mainshocks, 
aftershocks, and foreshocks some physical rather than 
ae Heherel definition is needed. That is not the purpose of 
this work. The purpose of this chapter is to study 
clustering of earthquake sequences at intermediate time 
scales, time scales of the order of years. A simple 
definition of mainshocks and aftershocks is used to 


eliminate the obvious short term clustering in the 
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Catalogue. The terms mainshock and aftershock used here 
Should be considered in light of this. 
In this work mainshocks and aftershocks were defined 


using the definition of aftershocks given in section 2.3. 


ta-tmsT(Mm); Ram<R(Mm); Ma<Mm 9 


were t is the origin time and M is the magnitude; indices a 
and m indicate aftershock and mainshock respectively; Ram is 
the distance between the epicenters of a mainshock and its 
aftershock; T and R are empirical functions assumed after 
Gardner and Knopoff (1974) and Keilis-Borok et al. (1980b). 
This defintion is simple. What is used is a Simple boxcar 
filter which is dependent on the magnitude of the largest 
event in a Sequence. 

Omori's law that the rate of aftershock generation 
falls off as the inverse of time results in a stong 
Criticism of this approach. If Omori's law is correct for 
all time then the total expected number of aftershocks in 


time T for a given event will be 
T 
NSO aye 10 
=] 


The summation of 1/t does not converge. Essentially this 
implies that each event, if it generates aftershocks, 
generates an infinite tail of aftershocks. The question 
arises, how do we know that a present earthquake is not 
Simply an aftershock of some earthquake which occurred in 


the perhaps distant past? Again at the root of this question 
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lies the exact question of what is the definition of an 
aftershock. After long periods of time the question of 
whether a earthquake is an aftershock of some previous 
event, which occurred at nearly the same location, or is 
perhaps due to changes in local conditions not related to 
the previous event becomes ambiguous. Perhaps the subsequent 
event could be more easily attributed to energy input into 
the system later in time than the previous events through 
some large scale driving force (i.e. plate tectonics). In 
such a case the event should not be termed an aftershock. 
Also, there 1S no guarantee that Omori's law is valid for 
very long time periods. 

If Omori's law is assumed to be correct for all times 
an order of magnitude estimate of the number of aftershocks 
of, for example, all magnitude 6 or greater events can be 
easily made. For the last 45 years there have been on the 
average approximately one earthquake with a magnitude 6 or 
more per two years in the Southern California catalogue. 
Although there are great fluctuations from event to event, 
an estimate of C in Omori's law for a "typical" magnitude 6 
earthquake can be made. Such an estimate gives an 
approximate value of C=2 days for aftershocks greater than 
magnitude 4. Assuming that we have a standard earthquake 
uniformly distributed every T, the expected number of 
"aftershocks" at time T which are generated by these events 


is given by 
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where T, is the average recurrence Deriod and L) is some 
Guuore SUT si eyelesmindicativerotathe lifetime of a 
fault system. For T,=2 years, C=2 days, L=1,000,000 and 
assuming T is small compared to nT, the summation gives 
approximately 10 events per year. This number is Slgnutice mc 
when it is considered that there are on the average 
approximately 70 events/year with magnitudes greater than 4. 
If Omori's law is correct for all time of the order of 10% 
of these events are explanable as "aftershocks" of previous 
large events. These event will, however, appear to be 
Statistically independent due to the different "mainshocks" 
to which they are associated and the long period of time 
from the "mainshock" to "aftershocks" of the same 

emain shocks 

An individual aftershock sequence might be modeled as 

following a time dependent Poisson distribution in which the 
mean number of events per unit time decays as 1/t. If such a 
sequence of events is observed after some period of time t, 
(i.e. the first t. days is not recorded) the series will 
appear random if time periods are observed which are small 
compared to t,. For a Poisson process the variance of the 
distribution of events in any time interval at will be equal 
to the mean number of events in that time interval. viet OF 
example, it 1S assumed that a change in the mean of 25% is 


Significant, a significant change in the variance of the 
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distribution of events following the 1/t time dependence 
will not be noticable unless many time intervals longer than 
t/4 are taken, where t is the time from the initial shock. 
Also, long enough time intervals must be taken to adequately 
estimate the mean number of events if differences are to be 
noticable. Clearly we can not tell if a sequence of shocks 
occurring in a present catalogue are "aftershocks", using 
Omori's law, of some event occurring at a time t, before the 
start of the catalogue if t, iS greater than the present 
length of the catalogue. Such sequences will appear to be 
random. We can not make any definite conclusion about 
properties of the sequences which show variations only at 
bamemscales: Omethe wordemotmor longer than the (length of ithe 
available catalogues. The length of the Southern California 
catalogue is of the order of 50 years. If we assume that the 
seismicity of the last 50 years is representative of the 
preceding 50 years we can estimate the number of magnitude 4 
"aftershocks" in the first year of the catalogue which are 
due to magnitude 6 events occurring in the 50 years which 
preceded the catalogue aS approximately 4. This falls to 
less than 2 after 6 years. This small number of events will 
not significantly effect the Statistics of the overall 
distribution of events. 

The effects of the finite time windows in the 
aftershock definition can also be considered. Consider a 
magnitude 6+ earthquake which occurs at the beginning of the 


catalogue. The number of aftershocks missed, if Omori's law 
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is assumed to be true, due to the finite length of the 
aftershock time window can be estimated for a typical event. 
Assuming the same parameters as before and taking T(M) from 
table 2.1 gives an estimate that 3 to 5 aftershocks with 
magnitudes greater than 4 will be missed in the 50 year 
period of the catalogue. This small number will have little 
effect on the overall statistics of the entire sequence of 
earthquakes. 

Consider another related question. It is well 
documented that large earthquakes on plate boundaries recur 
with Similar magnitudes at nearly the same location. It has 
also been suggested that nearly all the strain release along 
segments of plate boundaries characterized by the occurrence 
of large earthquakes is released by the large earthquakes 
themselves. Approximate return times for the largest events 
in a region can be estimated independently of the seismicity 
patterns uSing constraints from plate tectonics. A review of 
this was given in section 1.4.3. Suppose that all 
earthquakes occurring on these segments are aftershocks of 
the series of largest earthquakes. The question arises, 
after what period of time will the probability of an 
earthquake being an aftershock of the last large event be 
equal to the probability of being an aftershock of all the 
preceding large events. Assuming the large earthquakes occur 
at regular intervals of T, and all generate aftershocks 


initially at the same rate the time is given by 
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T=To/Et 

where the sum is taken over all cycles. If 100000 cycles 
are taken T is canoe aE Th (¢emneMaddle America she 
return times of large earthquakes have been estimated using 
plate tectonic considerations as 35 years (McNalley and 
Minster, i190). This means that 2.5 years after the 
occurrence of a large event we are faced with the ambiguity 
of saying that an event which occurred in the epicentral 
region of the last large event is more likely to be an 
"aftershock" of one the earthquakes occurring in one of the 
distant past cycles of the strain release. The question of 
what is an aftershock of what becomes a unanswerable. 
Clearly, some other physical criteria is needed if the 
SOnce Diao Mon mittihlemtailmoreat temShocksmis mused. 

In this chapter the possibility of long term clustering 
in earthquake sequences is studied. The time length of long 
term clustering which can be studied is restricted to 
periods which are short compared to the length of the 
catalogue. Nothing can be concluded about clustering which 
has time scales which are of the order of the length of the 
eatalogque on longer~- The Southern Californie cacaloguerspans 
approximately 50 °vears,. thisslimitsmenesstudyeorm Vongmtern= 
clustering to periods which are a fraction of this. Here 
ELONGstenm §CluSteringawills referetomciustering on time 


scales of 1 year to possibly 10 years. 
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A valid way to look for longer term clustering is to 
first filter out the short term variations which could 
possibly obscure the picture. This was the approach taken 
here. The purpose of the box car definitions for aftershocks 
was to remove the obvious short term clustering which occurs 
following a shock. In this study "mainshock" then refers to 
the largest event in one of these primary, box car defined 
clusters. "Aftershocks" and "foreshocks" are those events 
which occur respectively after and before the "mainshock" in 
the box car window. Foreshocks can have foreshocks and 
aftershocks can have aftershocks in this definition. The 
definitions here are obviously not unique. Mainshocks 
defined by the above criteria for a particular choice of 
R(M) and T(M) will not necessarily be mainshocks by another 
Gnivreria® 

Piesboxcarsivicerssused here are strong. Even if the 
1/t dependence suggested by Omori's law holds for very long 
periods of time most "aftershocks" occurring in the length 
of the catalogue will fall within the boxcar windows. The 
few aftershocks which are missed will be widely spread in 
time. These "missed aftershocks" will not affect the overall 


Statistics of the distribution of residual events. 
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3.3 Statistics of mainshocks 

The Southern California catalogue 1932 through 1980 
(Hileman et al., 1973; Freidman et al., 1976; and Subsequent 
Preliminary Epicenter Listings) provided the data for the 
study. Two separate mainshock catalogues were created using 
the two versions of the thresholds T(M) and R(M) given in 
fables vsima0Gms ieee nesthresholdsmin Table 2.1 were 
suggested in Gardner and Knopoff, (1974); Table 3.1 gives a 
Simplified version, used in the global test of premonitory 
clustering of earthquakes (Keilis-Borok et al., 1980b). 

The number of mainshocks in different time and 
magnitude ranges is shown in Tables 3.2 and 3.3. A map of 
mainshocks with Mm25.0 (using the threshold of Table 2.1) is 
given in Figure 3.1. It is interesting to notice the general 
quiescence in the occurrence of mainshocks from 1968 to 1978 


CHOU beE So 2)ie 


3.4 Clustering of mainshocks 

Knopoff and Gardner (1974) found that the distribution 
of mainshocks, using the definition of table 2.1),.1n Shone 
time intervals for the whole of Southern California is not 
distinguishable from a Poisson distribution. However, this 
does not eliminate the possibility that the sequences of 
mainshocks in smaller areas deviate from the Poisson 
Gistriburionsand, winweparticular sithatweclusters Of mainsShocks 
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Simplified windows for the identification of 
aftershocks (after Keilis-Borok et al., 1980b). 





M R (km) T (days) 
Shay OR et ha’ 5020 46.0 
450 e- eae 4 505 0 91.0 
Are gee bn4 5050 182.0 
Some Or 4 330), 18 365.0 
> eOn4 5080 1305-0 


Table 3.1 gives a simplified version of R(M) and T(M) used 
in the identification of aftershocks (after Keilis-Borok et 


aie yemoo0bD):. 
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Figure 3.1 shows a map of mainshocks for Southern California 
with magnitudes greater than 5.0. The aftershock thresholds 
of table 2.1 were used. The boundaries of the 14 subregions 


are also given. 
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Figure 3.2 gives a plot. of the cumulative number of 
mainshocks Ns, M24.0, from the beginning of the catalogue. 
5Ns gives the difference between the cumulative number of 
Mainshocks and the number predicted from the average number 
of events per year over the time interval ending 6 years 
before the date given. Note the general quiescence in the 
occurrence of mainshocks from 1968 to 1978. 
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of mainshocks in rectangular boxes of size Ax, Ay, S were 
considered. The parameters Ax, Ay, and s refer to the 
dimensions of the box in latitude, longitude and time 
respectively. 

Let n be the number of earthquakes in one of these 
boxes and let N be the total number of mainshocks in a 
Single rectangle Ax by ay for the entire period of time 
spanned by the catalogue. The two-dimensional histograms 
P(n,N) were tabulated. P(n,N) gives the number of times all 
rectangles Ax by Ay, which contained N events over the 
entire catalogue, had n events in the time intervals of 
duration s. The following values of parameters were assumed: 
AX = Ay = .4°, s= 1, 3 and 6 years; minimal magnitude M, = 
3.5 and 4.0. Two examples of these histograms are given in 
Tables 3.4 and 3.5. __ 

The distributions P(n,N) were compared to the expected 
Poisonian Milseriputs ones with the same average number of 
events per time window, using the Chi-squared criteria 
(e.g., Cramer, 1946). The values of P in each line were 
grouped in such a way that each group had at least 5 
members. The results are given in Table 3.6. Table 3.6 shows 
that the deviations between the observed distributions and 
the corresponding Poisson distributions are significant. In 
all but two cases the hypothesis thatthe distributions 
P(n,N) follow the Poisson law can be rejected at less than 
the .001 level of significance. Additional confirmation of 


the existence of clustering is that for n = 0 the observed 
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values of P are systematically greater than frequencies 


predicted by Poisson distributions. 


3.5 Clustering of mainshocks without foreshocks 

In the preceding section foreshocks were included in 
the analysis. The analysis demonstrated that the sequences 
of mainshocks plus foreshocks as defined by Tables 2.1 and 
3.1 show significant deviations from Poisson distributions. 
The deviations were such as to suggest a Significant degree 
of clustering, however the significance of the contribution 
of the foreshocks was not considered. The following question 
remains: is the observed clustering due simply to the 
presence of foreshocks? 

Gardner and Knopoff (1974) show that if aftershocks and 
foreshocks are eliminated from the catalogue by means of the 
criteria of Table 2.1, the distribution of numbers of 
mainshocks in 10 day intervals over all of Southern 
California is indistinguishable from a Poisson distribution. 

To test for local clustering of mainshocks with 
foreshocks removed, the Chi-square analysis of the preceding 
section was repeated on the catalogue of mainshocks as 
defined by Table 2.1 excluding foreshocks. P(n,N) was 
constructed for various values of Ax, Ay, Ss, and M,. The 
results are given in Table 3.7. 

The sequence of earthquakes (including aftershocks and 


foreshocks) shows a high degree of clustering. The removal 
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Table 3.4 gives an example of P(n,N) using the mainshock 


deLinitlonloretabplerc wand mainsShocksuwitheMes.5. 
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Number of Events in a Cell per 6 Years 
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Table 3.6 gives a Chi-squared test of the 

two dimensional histograms. In the table 

Ss gives the time window used in counting the 

number of events ina .4 by .4 degree square of 
latitude and longitude, n gives the number of degrees 
for the Chi-squared test, x gives the value of x% 
and P is the probability that a random variable 
which follows the Chi-square distribution with n 

n degrees of freedom is less than or equal to x. 
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Table 3.7 gives a Chi-squared comparison of the two 
dimensional histograms P(n,N), excluding foreshocks, to 
Poisson distributions. In the table s gives the time window 
used in counting the number of events, M the minimal 
magnitude, n the number of degrees of freedom, z the value 
of x?, p the probability that a random variable which 
follows the Chi-square distribution with n degrees of 
freedom is less than or equal to z. "I" gives the number of 
ranges which had variances greater than the variances 
predicted by a Poisson process and q gives the probability 
of 1 or more of I independent measurements being too large 
assuming the probability of measuring a variance which is 
too large is the same as the probability of measuring a 
variance that is too small. Clustering is indicated by large 
values of p and small values of q. 


Ch1-saquared test of distribution of mainshocks 


without foresnocks to a Porssonian distribution. 
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Area divided into 1 degree by 1 degree squares 


s M n Zz p I 1 q 
42 days 325 6 3.47 252 als) 2 654 
100 days <}65) Wel PASS VA 9382 27, WS O26 
100 days 4.0 5 3.16 325 17 4 -994 
1 year @).5 8) <\o; tas), SIE 997 24 20 -001 
1 year 4.0 14 een Great © 696 We 7 -834 
2 years Zhe (3) 292 44.00 989 24 22 - 0001 
2 years 4.0 Us AS. Se 9332 14 11 029 
2 years ABSS 5 2.74 260 7 4 =500 
3 years 253 Psy ZY 997 23 20 . 0002 
3 years 4.0 uth 19.24 943 22 9 O73 
3 years 4.5 6 93.46 851 6 4 .344 
4 years 325 10 22.94 -989 ie) iy, - 008 
4 years 4.0 fy She} 5 PAS Bois a ie | 8 113 
4 years 4.5 3) Grol 745 4 3 sae 
5 years 4.0 iO) aht/ 5 SIS) 9334 We 12 - 0002 
5 years 4.5 5 8.66 otyarts 6 6 O16 
6 years 4.0 YU KO) SERENE 10 8 O55 
6 years ar 4 8.86 SEIS) 6 5 109 





s n i p 

10 days €} 515) 4 Seon 720 
15 days 4.0 4 5S) G5Z 
100 days 3.5 1 8.15 300 
100 days 4.0 7 (Sie US) 482 
1 year 4.5 4 2.66 384 

1 year 5} (0) 3 2.13 454 
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of aftershocks and foreshocks in the manner described here 
and by Gardner and Knopoff (1974) can be viewed as a strong 
filtering process which removes the obvious mainshock - 
aftershock and mainshock - foreshock clusters from the 
sequence. This shifts the distribution of remaining events 
(mainshocks) away from a clustered distribution and towards 
a uniform distribution. 

A distribution which shows clustering will have a 
larger variance than a Poisson distribution, while a 
distribution which tends towards a uniform distribution will 
have a smaller variance. In terms of P(n,N) this means that 
clustered distributions should generally have systematically 
.more occurrences of cases in the tails of the distribution. 
In particular P(0,N) has an increased probability of being 
greater than that predicted by a Poisson distribution. For a 
uniform distribution the effect is the opposite. 

As a further test for clustering the following was 
considered: for the various values of Ax, Ay, s and M, the 
resulting histograms P(n,N) were grouped into ranges of N 
such that the ranges contained more than 25 cases (this was 
relaxed to 16 cases for s24 years due to the decrease in the 
number of cases). For each range N the variances of the 
distributions P(n,N) were compared to those predicted by a 
Poisson distribution .ginetTablesst7eules oi vesethemnumberms on 
ranges for a given Ax, Ay, S and M, and "i" gives the number 
of ranges which had variances greater than the variance 


predicted for a Poisson distribution. In Table 3.7 q gives 
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the probability of i or more of I independent measurements 
being too large assuming the probability of measuring a 
variance which is too large is the same as the probability 
of measuring a variance which is too small. A comparison of 
P(0,N) to the corresponding values predicted for a Poisson 
distribution gives almost identical results. Clustering is 
indicated by large values of p and small values of q. 

For the whole Southern California region with s=10 
days, 15 days, and 100 days and 1 year in the chosen 
magnitude ranges the results are in agreement with Gardner 
and Knopoff (1974). The distributions appear to follow a 
Poisson distribution. 

For Ax and say = .4° The degree of clustering is 
Significantly less than the case which includes foreshocks. 
The hypothesis that the distributions follow a Poisson 
distribution can not be rejected at the 95% level (except 
DOCmEMj ea SO wiS=-G8aNdeG years )ieThere is, “however sca 
suggestion of weak clustering for values of s of 3 and 6 
years for all cases. Apparently, a Significant amount of the 
observed clustering in the previous section was due to 
foreshocks. The tabulated cases suggest that differences 
between the distributions and Poisson distributions becomes 
progressively more significant as Ax, Ay, and s increase. 
This is indicated by increasing values of p and decreasing 
values of q aS Ax, Ay, and §s increase. 

These results Lace the hypothesis that earthquake 


sequences are clustered on longer time and distance scales 
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than is indicated by mainshock - Serato clustering, 
although the clustering is weaker than that suggested by the 
previous section. This clustering is evident on time scales 
and distance scales which are an order of magnitude larger 
than the time and distance windows used to eliminate 
aftershocks and foreshocks. 

The failure to detect clustering on distance scales of 
the order of the entire Southern California region is not 
Surprising. The catalogue of mainshocks will be dominated in 
number by events with magnitudes near the minimum magnitude 
cut-off. The statistical properties of the ensemble will be 
determined by this lowest magnitude range, larger events 
will not contribute significantly. Clustering is suggested 
to be a local effect. The time and distance scales of the 
clustering depends on the magnitude range considered. For 
the mimimum magnitudes considered here this range is an 
order of magnitude less than the size of the entire Southern 
California region. By using the entire Southern California 
region a superposition of several clustered distributions is 
obtained. The superposition of many clustered distributions 
may appear random. The clustering of mainshocks is only 
evident in time windows of the order of years. Due to the 
short duration of the catalogue (49 years) it is not 
possible to test for long-period, non-local clustering in 
the entire catalogue. 

The results presented here suggest that as the minimial 


magnitude is increased clustering becomes evident at larger 
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distance - time scales. Unfortunately the duration of the 
catalogue is insufficient to test this hypothesis. For M25.0 
time windows of 10 to 15 years might be needed. This would 
allow only 3 or 4 time intervals using the available 
catalogue. 

The hypothesis that Anere are related events in the 
catalogue over long periods of time and large distances, 
such as is suggested by the stochastic branching models of 
earnenquake soccuyrence, ersenoteiniconflict ewitth ithe 
observations presented here. If, for example, the Kagan and 
Knopoff branching model, in which all shocks can be thought 
of as "children" of those shocks occurring before them in 
time, iS correct no length of box car would be sufficient to 
eliminate all dependent events. The Kagan and Knopoff model 
was discussed in Section 1.11. In such a model clustering 
would always be observed at time scales longer than the 
length of the boxcar. The short length of available 
catalogues procludes the possibility of testing for the 


presence of such long term clustering. 


3.6 Part 1: Quiescence, swarms of mainshocks and strong 


earthquakes 
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3.6.1 Hypothesis 

The relevance of the clustering of earthquakes to the 
process of preparation of a strong earthquake has been 
reported in many studies. Some long-term premonitory 
seismicity patterns are directly formed by different types 
of clusters. The review of these patterns can be found in 
Chapter 1 as well as in Keilis-Borok et al. (1980b) and 
Kanamori (1981). Some long-term patterns (i.e. bursts of 
seismicity) have been shown to give good results in 
retrospective prediction tests. However, they do not 
indicate the place of the coming earthquake within the 
region. They indicate only a time-window (3-4 years for the 
bursts of seismicity) as do most other long-term patterns. 

Some other clustering-related premonitory patterns, 
which have been suggested, imply the following sequence: 
quiescence > activation and/or abnormally large clusters > 
strong earthquake. The several suggested patterns of this 
type use different definitions of quiescence, activation and 
clusters. A review of these patterns is given in Chapter 1. 
These patterns are designed to indicate not only the time, 
bee also the place of the coming large earthquake. However, 
so far they do not show "acceptable scores" even in 
retrospective prediction due to an excessive number of false 
alarms. 

Many abnormally large clusters, encountered in the 
study of the these premonitory patterns, are formed by a 


mainshock and its aftershocks (the aftershocks of only the 
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strongest earthquakes are eliminated from consideration in 

the definition of some patterns). In the following sections 
the hypothesis that the mainshocks themselves form a similar 
sequence (quiescence + activation and/or clustering) before 


a strong earthquake will be tested. 


3.6.2 Data analysis -- general scheme 

The above mentioned hypothesis was considered in the 
following way. 

(i) Southern California was divided in 14 areas, shown 
in Figure 3.1. These areas represent the zones of major 
faults; the location of these faults was taken from Hileman 
et al. (1973). Obviously the division, shown on Figure 3.1 
is by no means unique, especially where the faults are dense 
and not separated by zones of relatively low seismicity. 

(ii) The lists of mainshocks in each area were checked 
to see if the sequence {quiescence > abnormal clustering = 
Strong earthquake} occurred. For this purpose the following 
functions were tabulated: 

N(t,M,) -- the number of mainshocks in a sliding time window 
from oe) to t, with magnitude Mm 2 M,; 

R(t,M) -- maximal number of mainshocks (in the same 
time-window and magnitude range) with epicenters that can be 
placed within a rectangular cell Ax by Ay; 

N(t,M,.) -- average value of N(t,M,) for the time interval 
from the beginning of the catalogue to t. 


Swarms of mainshocks are defined by the conditions 
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RUeyemaxi1/2N(t) cl: N(t)2N(t) ie 


This same definition was used in Sapuroye etal -s (11977 )¥ton 
earthquakes in Italy, however the Italian Study sar cenot 
Separate mainshocks and aftershocks. 

Curie Hiunce vonseN( f)PEN(t) NUL) mandeswarmse were 
juxtaposed with the time of occurrence of StoONnGe\ Mme —6 4) 
mainshocks in each area. The hypothesis implies that strong 
earthquakes are preceded by a swarm (which is one of the 
types of abnormal clusters), and/or by a maximum of N(t) or 
N(t)/N(t), which are the measures of activation: either one 
has to be preceded by a mimimum of N(t) or N(t)/N(t), which 
would correspond to quiescence. Quiescence and activation 


ahceiOleLlOLrma lt yeaetined here. 


3.6.3 Data analysis -- results 

The catalogue of mainshocks, using the thresholds R(M) 
and T(MJ@ofL Table 3.1, was analyzed. To secure the 
completeness of the catalogue the minimum magnitude 
threshold M, was chosen as 4.0 except for regions of low 
activity. This more or less predetermined the following 
choice of parameters: sax = ay = .4°, s = & VeauSieG=o mel test 
the end of the catalogue the average number of mainshocks N 
in the 6-year windows was less than 6; the threshold M, was 
LoWwereGuouiiitCienuly sOmObDtaInvaneN eOLei rom 6mtore a These 
Parameters are not claimed to be optimal. The results of the 


analysis are represented on figures 3.3 to 3.10. The strong 
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(Mm 2 6.4) mainshocks are listed for convenience in Table 
Syl oa 

Vertical lines indicate the moment of strong (M 2 6.4) 
mainshocks. The numbers on the upper time scale indicates 
the year of occurrence and magnitude of the strong 
mainshocks outside the area represented on the plot. This 
makes it easier to check whether the occurrence of the 
pattern in one area is connected to the strong mainshocks in 
another area. 

Values of N(t), N(t)/N(t) and R(T), according to the 
GCeminiltiensormrnese tunctions,srerer to the end of the 
sliding 6-year time window. This window was moved in one 
year steps. The starting time for each area was shifted in 
such a way that the last strong mainshock falls at the 
beginning of a time window (the situation just before the 
strong mainshock is given in the window preceding the strong 
event). Accordingly, to avoid misunderstanding, the vertical 
line indicating the moment of the strong event was displaced 
by ~1 mm to the right. The actual moment is 
indistinguishable from the preceding time mark. The values 
of R are indicated only for swarms, each double circle 
Signifies the occurrence of a swarm. It should be noted, 
however, that consecutive double circles often correspond to 
one and the same swarm which 1S concentrated on a shorter 
time interval than 6-years. With the choice of 6-year 
intervals we could estimate N(t), N(t) and R(t) only after 


1937 as the catalogue begins in 1932. The estimates of N(t) 
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Se ee Se ewe +5 SRLee. ey 


Year 


Month 


Region 


ee NE cA ea A LENNON abate) 


1934 
1934 
1934 
1940 
1942 
1948 
Lo52 
1956 
1968 
1973 
£979 
1980 


Day Latitude Longitude Magnitude 
1 30 88233 =113550 6.6 
LZ 30 32.29 alas) S18) 6.5 
aL 32.00 114.75 These 
= the) BWA TUS =e) se ]0) sel 
10 oe B27 -116.00 ne) 
eZ 4 SSSI IS) Ubsliss 52 ¥e! 55 
21 35.00 12 Oz If 3° 
SRA ILS) rei l 5.92 6.8 
Bislaily) ~116.13 6.4 
34.41 -118.40 6.4 
10 5 32501 SAS) 6 Se 6.6 
5 22 Bi56 -118.79 (8) 


Owens Valley 
Laguna Salada 
Laguna Salada 
Imperial Valley 
Elsinore 
Southern Mojave 
Kern County 
Ensenada 

San Jacinto 
Santa Barbara 
Imperial Valley 
Owens Valley 


a 


Table 3.8 gives a listing of large events, M26.4, in 
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become stable much later. Negative focal depth is indicated 
in the catalogue for some earthquakes (most of them have 
epicenters in the Owens Valley and Imperial Valley areas). 
They were elliminated from the calculations represented in 
Figures 3.3 to 3.10. The results obtained by including these 
earthquakes do not differ significantly from those 


presented. 


3.6.4 Summary of the results 

The areas considered can be divided into two groups. 
The first group consists of 7 areas for which a lower 
threshold magnitude of M, = 4 could be maintained (Figures 
353 tO)>- 6). The) data for@thissgroup seem to confirm the 
expected pattern. The three last strong earthquakes -- Owens 
Valley, 1980; Imperial Valley, 1979; and San Fernando, 1971 
-- are preceded first by quiescence (represented by clear 
minimums of N(t) and N(t)/N(t) and then by the swarms. The 
San Miguel earthquake, 1956, is preceded by a swarm, but not 
by a pronounced minimum in N(t)/N(t). N(t) does show a 
minimum from 1941 to 1952, however N(t) had not achieved a 
stable value. For the 1948 Desert Hot Springs earthquake the 
results are negative. It was preceded by activation (an 
increase in N(t)) but not by a swarm or quiescence. Swarms 
were not diagnosed in the S. Mojave area at any time. For 
the Parkfield area the results are positive in the sense 
that neither the quiescence > swarm pattern or strong 


earthquakes occurred there (there were no false alarms). The 
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i r tivation 
i 3.3 gives the test of the Quiescence + ac 
ri A het rae the Owens Valley and Imperial Valley regions. 
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Figure 3.4 gives the test of the quiescence + activation 
hypothesis for the Santa Barbara and Ensenada regions. 
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Figure 3.5 gives the test of the quiescence + activation 
hypothesis for the Southern Mojave and Laguna Salada 
regions. 
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area is presently in a state of quiescence. For the Laguna 
Salada area a "false alarm" occurred in 1957. This "false 
alarm" may be connected with a general quiescence - 
activation associated with the neighboring Ensenada region 
before the San Miguel earthquake of 1956. 

The second group includes the 7 remaining areas where 
mainshocks with magnitudes less than 4 were considered 
(bigures = 3./@to 3.10)MsThe expected pattern was not observed 
in these areas. Stronger mainshocks, Borrego Mountain, 1968, 
and Kern County, 1952, were preceded by swarms, but not by 
clear quiescence. The same was true for the strongest 
mainshock in the Southern Sierra region, M = 6.3. In the 
Riverside area we have a "false alarm" and in the Elsinore 
area — two "false alarms”. In spite of the genéral poor 
score for these regions, the recent swarm in the San Jacinto 
area deserves attention, since it was preceded DYgamc leat 
minimum of N(t) and N(t)/N(t). This swarm includes a 
mainshock in 1981. The adjacent Southern Mojave area shows a 
Similar pattern. The two areas, San Jacinto and the Southern 
Mojave, are on opposite sides of the present seismic gap 
Hoendir ved invthatcherset al.) (1975) eeThisaqapsrs tocatedman 
the San Jacinto fault in the northwest corner of the San 
Jacinto region (see figure 3.1). 

In both groups the score could probably be improved by 
the change of parameters, first of all by the use of the 
thresholds from Table ™2.1¢and by narrowing Ax and Ay in ‘the 


definition of a swarm. However, aS a posteriori improvement, 
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Figure 3.7 gives the test of the quiescence + activation 
hypothesis for the San Jacinto and Kern County regions. 
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Figure 3.8 gives the test of the quiescence > activation 
hypothesis for the Southern Sierra and Los Angeles Basin 


regions. 
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Figure 3.9 gives the test of the quiescence eC —aveation 
hypothesis for thesRiversivde andseousinore sreqions. 
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Figure 3.10 gives the test of the quiescence + activation 
hypothesis for the Northern Mojave region. 
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it would hardly change the conclusions. Figure 3.11 gives a 


rough qualitative summary of the results for all areas. 


JO, OeCONGCLUSTONS-) pant | 

Mainshocks of Southern California, while forming a 
Poisson sequence in this region as a whole, are sometimes 
clustered in the space-time domain. This clustering will be 
more thoroughly examined in the following sections of this 
chapter were a more general definition of clusters is used. 

There are qualitative confirmations of the hypothesis 
that’ mainshocks of magnitude Mm = 4 form a "quiescence > 
Swarm" pattern before a mainshocks of Mm = 6.4 in the same 
fault zone. For a more qualitative test of the hypothesis 
the follwing are required: quantitative definitions of 
quiescence; more objective separation of the fault zones, or 
a method to identify quiescence and clustering without such 
separation; and a statistical model of the sequence of 


mainshocks. 


3.7 Part 2: cluster analysis of mainshocks in Southern 
California 

The first part of this chapter demonstrated that the 
sequence of mainshocks in Southern California, as defined by 
the simple criteria of Table 2.1, shows a significant degree 
of local clustering. Periods of low clustering (quiescence) 


followed by high clustering (Swarms) were suggested to be 


















. 


wiiviteove: Gay opnecs ¢ tht etjRagee ae 
1s Oey ye iets © ha. azilevp asin 


er roraui snes ele 


PL my «6 sri Qt aD aa \ 


tds nj] svreageg none fot 
veta-¢oa0aG off Ai Begesgayee 
gtingae yleégubwods saan 

vh e's LOeTSang 


s18ay 
‘salen (3082 


“kik ‘Jeqq “iseawe 


SiG? tga 

- @ ([HoOnenee tag 

TT ticoehl af Bodgens 
Oo: 228s 


etvodeniant 





Qgeviiet wi a lanise Do af@yfons +etenls 4% eee 


aiaaliles - 


cr) me 
v 


i or sAinhe> ip Qano Aika. 3a 1386 ee f 


Sel, 
_ ' 





ron an a, Neer 





110 


60 10 80 







6.5 ape Mt. | 


opoodgo 





63 ERSENADA) 


1qdoo 
6 


l dq qokbooo ig 
iba eidies oar sere 
| a) RAAT Het 











LO eet | r LAGUNA SALADA 
-| 7.1 Delta be est pe eS . i 
| 
[det [| Parry aa 
se [bs Liens 
| 


40 


Va 


Figure 3.11 gives a rough summary of the results of of part | 
1. Periods of quiescence is indicated by dark circles and 
Swarms are indicated by open circles. The minimal magnitude 
for the region and large earthquakes occurring in the region 
are also indicated. 
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premonitory to a future strong earthquake near the quiescent 
> swarm region. The swarm model of clustering was very 
Simplistic and the characteristics of individual swarms were 
not considered. This section uses cluster analysis to 
investigate the clustering of the mainshocks on a more 
formal level. The sequences of clusters in several areas are 
examined. 

Cluster analysis is a means of classifying a group of 
individuals (objects, cases, specimens etc.) into sub-groups 
such that the "Similarity" between individuals within each 
Sub-group is high while the "Similarity" between individuals 
of different sub-groups is low. Although the application of 
cluster analysis is most familiar in the social sciences, 
the need for classification schemes is common to all fields 
of study. With an adequate classification scheme, the amount 
of data necessary to describe a population can be greatly 
reduced (i.e. the population can be characterized by the 
characteristics of the sub-groups rather than the 
characteristics of all members of the population). Remote 
sensing is an example of an area in the geophysical sciences 
which nae made extensive use of cluster analysis (for 
example: Holmes and MacDonald, 1969; Landgrebe et al., 1969; 
and Landgrebe and Phillips, 1967). The classification 
problems encountered in pattern recognition and those in 
cluster analysis are very similar. 

The amount of literature on cluster analysis is large. 


The nature and technique of applications of cluster analysis 
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vary Significantly from problem to problem. Only the 
concepts of cluster analysis which are used in this chapter 
will be discussed. For a comprehensive review, the 
interested reader is refered to to the following sources: 
Wishart (1978); Jardine and Sibson (1971): Duran and Odell 


(1974); and Spath, (1980). 


3.7.1 The clustering problem and mainshocks 

Pet Om-micqnmads, Gs, +... denote the set of N 
mainshocks in Southern California. This set forms the 
population for the study; the individual mainshocks are the 
elements. Each earthquake is characterized by a set of 
observables X which are used to define the Similarity or 
dissimilarity between elements. In theory these observables 
could include any characterisitics of the mainshocks such as 
location (latitude, longitude), focal depth, time of 
occurrence, focal mechanism, magnitude etc. Only latitude, 
longitude, and time of occurrence will be used in this 
study. Let J be an integer less than N. The problem posed by 
cluster analysis is then to determine J subsets (clusters) 
of Q such that each element is assigned to one and only one 
cluster and such that individuals assigned to the same group 
are "Similar" while individuals in different groups are 
iditferent™, the measure of ’similarity and«dissimilarity 
being determined by the set of observables X . One group may 
be a residual set which contains the individuals that are 


not Similar to other groups of individuals. 
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The terms Similarity and dissimilarity need to be 
quantitatively defined. Needed is a measure of the 
Similarity or dissimilarity between two elements, an element 
and a cluster, and between two clusters. Many such measures 
have been proposed (see for example Wishart, 1978). The most 
familiar are measures of dissimilarity based on the ordinary 


Euclidian distance which are used in this study. 


3.7.2 On Uniqueness 

Any classification scheme devised is inherently 
non-unique. The most obvious source of non-uniqueness enters 
in J the number of clusters determined. Obviously, a good 
solution.depends on a reasonable choice of J and even for a 
given J there is no guarantee of uniqueness. Some criterion 
is needed to measure the "goodness" of the solution. If such 
a criteron can be properly constructed, the partitioning can 
be formulated as an optimization problem based on the 
desired criterion. 

Even when such a problem can be formulated, the problem 
of non-uniqueness still exists and different "goodness" 
ena could give different results. The problem of the 
relative scaling of the various observables is another 
source of non-uniqueness. For example, when time and spatial 
coordinates are included, the relative weights of the two 
must be defined in the dissimilarity criteria (some unit of 
time must be equated to some unit of distance). In general, 


the clusters will vary under different scales. Also, some 
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computational algorithm is needed to solve the clustering 
problem. There is no guarantee that a computed solution will 
give an absolute optimum and not just a local optimum. 

A direct way to solve the clustering problem for a 
given optimization criterion is to evaluate the optimization 
criteria for each choice of clustering alternative. The 
unique solution would then be the choice yielding the 
optimal value. However, this method, termed clustering by 
complete enumeration (Duran and Odell, 1974) is impractical 
for any reasonable problem. The number of possible 
pareiulonssOmeatset sor nm objectsmintow subsets none of 
which is empty is given by S(n,j), where S(n,j) is the 
Stirling number of the second kind (Abramowitz and Stegun, 
1968). 

For a Six year interval there are on the average 200 
mainshocks with magnitudes greater than 3.5 in Southern 
California. If ten clusters are desired, complete 
enumeration would require the evaluation of approximately 
10'*?* cases. Some efficient strategy is needed to limit the 
amount of computational time required to some feasible 
values 

Many efficient clustering algorithms have been 
suggested in the literature. Three different strategies are 
used in this study, density analysis, hierarchical fusion, 
and iterative relocation. These algorithms are taken from 
CLUSTAN, a cluster analysis Pa peva published by Edinburgh 


University (Wishart, 1978). 
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JerinvgeronsSmandenotationsertrom wrshart 


total number of individuals. 
number of variables per observation. 


total number of clusters, 
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individual. 


eluster. 


Nie Cem imbe wrote nd v1 GUallcumin the ath 


oer NemameanavaluesoLeathes1th variable ror the 


Kenmeelusterm. 


the ith 


Gen cous 


cluster. 


cluster. 


Z1 .... a vector of length m giving the center of 


Dij .... the ordinary Euclidian distance between the 


OmeLnemttoneandajthecluster. 
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WicheC UStercmimandej@are joined (ives 1(i,j) = Brr(ity) — 
petal ye te jae ay) 

The error sum of squares, refered to as the 
within-group sum of squares by Duran and Odell (1974), is 
the sum of the distances from the center of a cluster to 
each individual contained in the eicecar? As such it is a 
measure of the scatter around the cluster center. Two of the 
cluster techiques used here (hierarchical fusion and 
iterative relocation) use this functional as the 
optimization criterion. Solutions are sought which minimize 


RNeéstovalwerrorn, sum of squares for J clusters. 


3.7.4 Algorithms used 


3.7.4.1 Hierarchical fusion 

Methods of hierarchical fusion start with N clusters 
(each individual is a starting cluster). In each of N-1 
cycles the two "closest" clusters are fused forming a single 
cluster. The final step consists of a single cluster of size 
N. The method can be thought of as giving the sequence of 
overlapping classifications which define the hierarchical 
manner in which each individual and/or group of individuals 
joins the entire population. The method of hierarchical 
fusion used here was first suggested by Ward (1963). In 
Ward's method each cycle combines the two clusters whose 
fusion yields the least increase in the total error sum of 


squares. Ward's method finds minimum variance spherical 
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clusters (Wishart, 1978). 


3.7.4.2 Iterative relocation 

Iterative relocation finds a local optimum for J 
clusters, starting with an initial classification of the 
individuals into J initial clusters. Each individual is 
removed from its parent cluster. The similarity between the 
individual and each cluster is compared and the individual 
is placed (relocated) in the cluster which results in the 
least total error sum of squares. The process is repeated 
until there are no relocations during one complete cycle (a 
scan of each individual). In general the process can be 
shown to converge for reasonable data. However, convergence 
will not in general find a global optimum in the error sum 
of squares for the J clusters. Rather, convergence will be 
towards a local optimum in the neighborhood of the initial 
elustern classification (Wishart, 1978). A useful grouping by 
iterative relocation depends on a good starting 
GlassiiiGationuwhich iS not toomdistant. -cromecherfanal 
desired solution. In this study iterative relocation will be 
used to refine initial classifications generated by other 


methods. 


3.7.4.3 Number of natural clusters: density analysis 
Suppose @(x,y,t) is some continuous function which 

gives a theoretical model of the density of epicenters in 

Space and time. Wishart (1969) suggests that the number of 


clusters for the lowest level of clustering at which the 
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Classifications are "natural" or taxonomically significant 
1S Given by the number of local maxima in Such a density 
function. The lowest level of classification corresponds to 
the hierarchical level with the maximum number Ofaclusterns: 
This will be taken as the number of clusters of interest in 
this study. 

To determine the number of maxima it is necessSary to 
construct an estimate for the continuous density function 
using the discrete set of epicentral locations given by the 
data. Some averaging process is necesSary. One method is to 
divide the space into small cells in Tati tudes longitudes and 
time. An estimate of the density at a location could then be 
given in terms of the number of events Pen cell isa sano 
a good method since local maxima which are close could be 
fused in the averaging process. Wishart (1978) Suggests that 
the density estimate at a given location be based on the 
average distances to some number of nearest events. 

A pseudo density G@(x,y,t) is defined here as follows: 
let UDite #Di 2 ese. Dik Hebe tthe tdistances from the ith event 


to its k nearest neighbors. Then d(x,y,t) is given by 


nO yes D 16 
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The number of nearest neighbors k used in the above 
estimate determines the amount of averaging used in the 


calculation of qd. An increase in k is analogous to 
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increasing the volume of a cell over which average densities 
are estimated. Two local maxima which are close will remain 
distinct if each is a result of more than k/2 events. 
Physically @ gives an estimate of the density in the 
neighborhood of an individual by using the estimated radius 
which encloses k/2 neighbors, this radius being estimated as 
the average distance to the k nearest neighbors. Wishart 
(1978) points out that this procedure will give smoothed 
density estimates. The densities at the maxima will be 
underestimated and the densities at saddle points will be 
overestimated. Maxima due to less than k/2 individuals will 


not be detected. 


3.7.5 Initial clustering constraints 

In order to search for clusters of mainshocks some 
preconception of what constitutes a cluster is necessary. 
Thesresults of part 4 of this chapter will be used) tor this 
purpose. In the first part of this chapter Swarms were 
defined as groups of mainshocks occurring within 6 years or 
less of each other, and occurring within a rectangle of 14° 
by (4°, The results of Section 3.4 can be used suggest more 
appropriate time and distance scale for the clusters. 

Section 3.4 suggests that "mainshocks" (as defined 
here) with magnitudes greater than 4.0 show a significant 
degree of clustering for time periods of the order of 3 to 6 
years and distances of the order of 100 km. Somewhat shorter 


time and distance scales were suggested for the case with 
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minimal magnitude equal to 3.5. The DOSStDINIIty of ever 
longer term clustering in the catalogue could not be tested. 
The question asked here is, what sequences of earthquakes 
give rise to clustering characteristics for these time and 
distance scales? Although this observation is not well 
constrained and does not give a optimal estimate of possible 
Spatial and temporal sizes for the clusters, it can be used 
to give an order of magnitude estimate of these Sizes. In 
Particular; ethes results, of, Section 2.4 give an order of 
magnitude estimate for the scaling of time relative to 
distance. This suggested normalization, to look for 
clustering as was indicated in Section Shit eel S eS OMO RV eaGS 


equalmtom100kme 


3.7.5.1 Dissimilarity criteria 
In terms of latitude (y), longitude (x), and time Dij 


.1S given as 


: ; heh 
exp) VEY) itiety) | 2 7 
' Sl X5 YG t5 J 


where xo, y. and t, are normalization factors. The 
constants xo, yo, and t, are chosen so that at the center of 
the region (35° latitude) a specific unit of distance, r,, 
is equal to a specific unit of time. Examples used here 
include: 50 km. = 2 years; 50 km. = 3 years; 75 km. = 2 
years; and 30 km. = 2 years. These result in normalizations 
of 100 km. = 6.67, 6.0, 4.0 and 2.57 years. There is no 


physical basis for these choses other than the observations 
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QSaSeation 3.4. 

AS an example xo=.548393°, y,=.449507°, and t,=2 years 
give the following properties to Dij: 

wane SNescentenrorethesregionm(35° latitude) two tyears 
of time is equated to 50 kilometers of spatial distance. 

~~ Two events occurring at the same place and separated 
DYMO yeansmwuiwtihegivesa Dig of &:; 

-- At the center of the region two events occurring at 


the same time and separated by 150km will give a Dij of 1. 


3.7.5.2 Residual events 

The results of Section 3.4 suggest that clustering of 
mainshocks, although present, may be weak. Some earthquakes 
llayebeslatevemoved=shromeareas Of clustering. local clusters 
are sought. These local clusters need not include all 
events. For this reason a residual set is created which 
includes those events not found to be sufficiently similar 
to the clusters found. 

For practical purposes clusters with very few 
individuals are not of interest. Only clusters with more 
than five individuals were considered. This minimum number 
was chosen arbitrarily. 

To create the residual set the following constraints 
were applied to the events: 
ljeliea clusterawvasefound toscontainsless than five 
individuals the individuals were placed in the residual set 


if the individuals were not sufficiently close to another 
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Pomerere 

2. Let k be the number of the cluster to which the ith event 
1s assigned. If Dik > T the event is moved to the residual. 

Dik is the dissimilarity between the ith event and centroid 

of the kth cluster with the ith event removed and T is some 

threshold value. 

Events which are in the residual set do not necessarily 
remain in the residual set. If the centroid of the clusters 
changes so that the minimum dissimilarity between an event 
and the closest cluster is less than T , the event is 
assigned to that cluster. At the end of the clustering 
process the residual set is searched for the existence of 
missed clusters. The value of T limits the time and Spatial 
distance that a cluster can span. For example;manchoicevot aT 
= 1.00 limits’the time span of a cluster to less than 3t, 
years and the spatial span of a cluster to less than 3r, at 


the center of the region. T will equal 1.00 in this study. 


3.7.5.3 General clustering procedure 

The procedure used to find clusters is as follows. A 
density function @ was constructed as indicated in section 
3.6.4.3 to estimate the maximum number of clusters for the 
analysis. Several values of k were used to construct @. The 
results presented here are based on k=7 for M24.0 and k=11 
for M23.5. The number of local maxima ( J ) in the density 
function were determined. This number was equated to the 


maximum number of clusters searched for in the study. In 
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general the number of clusters that resulted was less than 
Js several or the local maxima were found in areas of low 
density and did not survive the threshold criteria. 

Two separate procedures were used to generate initial 
classification of all objects into J initial clusters. 
First, events were grouped into J clusters by assigning each 
event to a cluster that had an assumed centroid 
corresponding to the location of the local maximum (i.e. 
each event was assigned to the nearest local maximum). 
Second, Ward's method of hierarchical fusion was used to 
generate J+I cluster categories. Ward's method could result 
in the splitting of a cluster due to just one local maximum 
into several clusters. "I" was chosen by looking at all 
hierarchical classifications with J or more clusters and 
picking as the desired initial classification the 
hierarchical level with the minimum number of clusters which 
had at least one cluster for each dense point. In both 
methods, all events are classified into clusters (no 
threshold criterion was used). 

Iterative relocation was then used to refine the two 
initial classifications by creating a residual set and the 
final results were compared. The results of the separate 
procedures differed very little. The results using 
hierarchical fusion to generate initial classifications are 
presented here. In general they resulted in final 
classifications with lower total error sum of squares than 


did the other procedure. The entire process was repeated on 
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the residual set to search for missed clusters. 


3.7.6 Results 

The output of cluster analysis consists of a listing of 
the J groupings of the earthquakes, each grouping 
constitutes a single cluster. For Southern California such 
listings are rather large, consisting of 100 to 300 clusters 
(depending on the clustering parameters) of the 
approximately 2000 mainshocks with M23.5. A complete listing 
of cluster groupings will not be given here. Rather, several 
illustrative examples will be discussed. A listing of 
several complete cluster analysis groupings in computer 
readable form is provided for the interested reader on tape 
number 007917, volume label LMRX10 stored at the University 
of Alberta Computer Centre. A copy of the tape is available 
upon request from Dr. Edo Nyland, Department of Physis, 
University of Alberta. The tape consists of 9 separate 
files. File 1 contains a complete description of the 
contents and formats of all other files. Listings and 
explanations of several Fortran programs needed for the 
interpretation of cluster analysis results are also 
contained on the tape along with the following cluster 
analysis examples: 
1. Cluster analysis for the northern section of Southern 
California (see Figure 3.12 for region boundaries) with 
M24.0, k=7, and a normalization of 50km=3years. 


2. Analysis for the southern region of Southern California 
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With M24.0, k=7, and a normalization of 30km=2years. 
3. Analysis of the northern section of Southern California 
with M23.5 , k=11, and a normalization of SOkm=2years. 
4. Analysis of the southern section of Southern California 
with M23.5, k=11, and a normalization of 30km=2yrs. 
The interested reader is invited to analyze these examples 
Or to generate other classifications using the procedure 
outlined on the tape. 
Two examples of ences clusters are given in Table 
Seo eaDOec a Onandmrigurels. |2eDuelto computing cost 
considerations it was necessary to divide the Southern 
California region into two segments for the cluster 
analysis. The boundaries of these segments are given on 
figure 3.12. Table 3.9 gives a listing of the centroids of 
clusters found in the Northern section with normalization 
factors of x,=.548393, y,=.449507, and t,.=3 years (i.e. 50 
km. equal to 3 years). Also given are the standard 
deviations in latitude, longitude and time for each cluster 
and the total Sigma of each cluster. Sigma is defined as 
iS cme “ 
kei 
Sigma is roughly proportional to a summation of the 
fracture areas of the earthquake sources (Keilis-Borok and 
Malinovskaya, 1964). 
jablems SURO ivesudesimmlanslistingmromeciusters = nethe 
SOUCNELMecect lOnmwithex,=.3 29036), eyn= 26007040 andat, =2 


years (i.e. 30 km. equal to 2 years). Figure 3.12 shows the 
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centroid for al) events 

PESJ@) ES. GH) VeE43'257 F=1187 1 0.8024 1954.35 13.8063 
centroid for the residual events 

52 soo 1.5295 -118.90 1.1945 1951.99 15.0746 


centroids of the 33 clusters 
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8B 34.47 OWE) = WexXe).S 0.1669 1937.78 1.8129 Ue} o SIS}<} 
SG ZS O20 2am liga O) 0.4279 1939.13 0.8896 Tf atsheyy 
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7 34.90 ORS Soa SEOS ORO SS oaieos 1.3806 22.061 
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GS SoS Qo BANA SN 7 BIE ORNS SS a O72e97 1.9467 8.586 
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Table 3.9, qives a tabulation, of the centroids 
found in ae northern ‘regionsof Southern) Galifiorpnia for the 


case M,=4.0 and 50 km: equal” to™3"years (ive.Nty=3) years, 
Tee oto393se and yp=.449507.)" 
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centroid for al) events 


Ses} lel. 9793 1.2433 -116.14 1.0142 TSS Siti, 14.2710 
centroid for the residual events 
yey hed ae) 1.6645 -116.53 1.4260 1955.49 14.7124 
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Figure 3.12 gives the locations of the cluster centroids 
listed in tables 3.9 and 3.10. The boundary between the 
northern and southern regions is indicated by the dashed 
ne. 





Meats} 


locations of the above cluster centers on a map of Southern 
California. 

From Figure 3.12 it can be seen that the cluster 
centers themselves form groups in space. The groups of 
cluster centers were instumental in the choices of 
boundaries for the divisions of the Southern California 
Catalogue into the smaller regions analyized in Part 1 of 
this chapter (Figure 3.1). The remainder of this chapter 
will give a discussion of the sequence of clusters in 


several of these subregions. 


3.7.6.1 Clusters in the Owens Valley region 

The» Owens Valley region (Figure 3.1) is the most 
isolated of the regions. Figure 3.13 gives a time magnitude 
plot of mainshocks, M24, in the Owens Valley region defined 
by figure 3.1. Earthquakes not found to be members of 
clusters are shown with a square symbol. The clusters were 
determined with a normalization 50 km.= 3 years, and K=7. 
The time spans of the various clusters are indicated by 
horizontal bars above the plot. Figures 3.14a and 3.14b a 
latitude, longitude plot for each cluster in the region. The 
quiescence and activation before the 1980 Owens Valley 
earthquake is clearly evident in the plots. The quiescence 
and activation is indicated by number of events, energy 
release, and compactness of the clusters. The cluster 
occurring between 1971 and 1976 has the fewest members and 


lowest Sigma of any cluster in the region. The cluster also 
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Figure 3.13 gives a time magnitude plot of earthquakes in 
the Owens Valley region (figure 3.1). Earthquakes not found 
to be cluster members are indicated bu a square symbol. The 
times spanned by the clusters are indicated by horizontal 
bars. 
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Figure 3.14a gives latitude longitude plot of the clusters 
found in the Owens Valley region. 
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Figure 3.14b gives a continuation of Figure 3.14a. 
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Shows a large spread in distance and time compared to the 
other clusters. The cluster which includes the Owens Valley 
event has just the opposite characteristics, most members, 
largest energy release, and it is one of the most compact 
clusters in time and space. There are no other similar 
quiescent > activation periods for the Owens Valley 


sequence. 


3.7.6.2 Clusters in the Southern Sierra region 

The Southern Sierra region (Figure 3.1) is an area of 
much lower seismicity rate than is Owens Valleyemraguresse5> 
gives a latitude - longitude plotyoft clusters in the region. 
The clustering parameters are the same as they were for 
Owens Valley. The clusters found are very weak containing 
few members (5 or 6) and are not very localized in Space. 
There is one exception. The cluster containing the 1946 
magnitude 6.3 earthquake contains 13 members. Two 
earthquakes in the cluster occur after the 6.3 event, both 
in 1948. The cluster is much more concentrated in Space than 
are other clusters in the region. The main activity in the 
cluster occurs on the lineament north of the epicenter of 
the magnitude 6.3 event (see Figure 3.1). It is this group 
of epicenters which generates the swarm found in Part 1 for 


this region (see Figure 3.8). 
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Figure 3.15 gives a latitude longitude plot of the clusters 
found in the the Southern Sierra region. 
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3.7.6.3 Clusters in the San Jacinto region 

Figure 3.7 shows a clear quiescence followed by a Swarm 
in the San Jacinto region during the last time windows of 
the catalogue. This pattern is of particular interest due to 
the suggested current high seismic potential of the region 
(for example Thatcher et al., 1975). 

BUGULeS = Au ba ,ers | 6D sand 93 s1OCmOsvecmaG latitude — 
longitude plot of clusters for M23.5 in the region. A 
clustering normalization of 30 km. equal to 2 years was 
used. Figure 3.17 gives a time magnitude plot of the cluster 
members. The clusters display several interesting features: 
ie nemlargesteevent inethe reqiony, the 1942 magnitude: 6.5 
event, was preceded by anomalous activity, although the 
suggested quiescence > activation pattern is not clear. 
Pemelicecqulescences= activation patternmoccunring satethe end 
of the catalogue is evident in figure 3.16. 

Activity before 1937 is low and the clusters for those 
periods are not spatially concentrated. During the period 
1941 to 1946 two very tight clusters were found. One of 
these clusters contains the 1942 magnitude 6.5 event. 
Fourteen of the 22 events in these clusters occur before the 
1942 event and are concentrated on the same fault as the 
magnitude 6.5 event. The only other cluster which is so 
Spatially concentrated is the last cluster. 

The 1968 magnitude 6.4 event is not precéded by a 
Similar pattern. The epicentral region of the 1968 event is 


quiet and no activation is evident before its occurrence. 
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Figure 3.16a gives a latitude longitude plot of clusters 
found in the San Jacinto region. 
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Figure 3.16c gives a continuation of 3.16a and 3.16b. 
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Figure 3.17 gives a time magnitude plot for members of the 
clusters in the San Jacinto region. The times spanned by the 
clusters are indicated by horizontal bars. 
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Interestingly the only events in the clusters which occur 
near the future epicenter during the period 1954 to 1968 
form a "doughnut pattern" around the epicenter. The doughnut 
pattern is formed by events occurring from 1961 to mid 1963, 
seven years prior to the major event. The entire fault 
System is very quiet after 1968 until 1975 ness the 
Spatially concentrated activity begins on the north section 
of fault. Events in the cluster occurring from 1972 to 1976 
are mostly on the fault system to the southwest. The spatial 
concentration of events in the last cluster looks 
Superficially very similar to the activity preceding the 


1942 major event. 


3.7.6.4 Other regions 

In the preceding three examples the quiescence > 
activation patterns observed in part 1 are also visible in 
ENesclusrers thus was true for other areas also. Clusters 
in the Imperial Valley region shows the quiescence > 
activation pattern before the 1979 magnitude 6.6 event. 
Similar results are seen in the Santa Barbara region before 
the magnitude 6.4 San Fernando event and in the Ensenada 
region before the 1956 magnitude 6.8 San Miguel event. 

Regions in which the quieScent + activation pattern 
were not observed show similar cluster characteristics. 
Clusters in the Kern County and Desert Hot Springs regions 
show no obvious patterns before the magnitude 7.2 Kern 


County event and the magnitude 6.5 Desert Hot Springs event. 
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The low seismicity regions in the south, North Mojave, 
Riverside, Elsinore, and Los Angeles Basin have few cluster 


centers. 


3.7.7 Density representation of clustering 

Other means of characterizing the clustering of 
earthquakes are possible such as the pseudo density function 
Ceceninedminmcecimonmseo 04.30. Fiqureses alse 3.19, 3205 sand 
3.21 show examples of @ for mainshocks with M23.5. Here @ is 


defined in the following way: 





Daz [En=2v*, Onn ¥ir i 
hs ¥5 


where Dni is the distance from the nth event to its ith 


nearest neighbor and is given by 


dixy)=10(18/ 5p, ) 
I 
3} 


20 


x refers to latitude and y refers to longitude. The 
catalogue of mainshocks was divided into 6 year intervals 
and the function @ was computed for these intervals. 
Specifically d 1s the density functuonmorysectioneseonees 
without time, calculated for 6 year time intervals with 
k=ige Thejdistance Din 1S normalized sso ithatel00mkme sequals 
one unit of distance in the equation. 

Figure 3.18 gives an average of @ for Southern 
California using sliding 6 year time windows stepped one 
year. The average was constructed as follows: @ was 


calculated for the six year intervals from 1932 to 1981. The 
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Figure 3.18 gives a contour Plot of the average @ for the 
Southern California region uSing sliding 6 year intervals 
Stepped one year. Anomalous values of @ have been eliminated 
from the average. ; 
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Figure 3.20 gives a 
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Figure 3.21 gives a map of @ for the years 1975 to 1981, 





148 


values of @ at the epicentral locations were used to 
estimate values of @ at regular grid points spaced .25°. 
These estimates were bases on weighted averages of the 
epicentral points using a weighting factor which decreased 
ecm, (0 lseance stromatheagnidepoint)eAnmoverallwaverage for 
each grid point was constructed. The average values were 
then recomputed after eliminating all values found to be 
anomalous. Anomalous values were defined as values at the 
grid points which differed by two standard deviations from 
the grid point means. The Surface II graphics package, as 
implemented on the University of Alberta computer, was used © 
to generate the contour mapsS and grid estimates (Sampson, 
797.5.) -2 

Figure 3.19 and 3.20 give maps of the @ for the periods 
PAS eo lands! 950 etom.95cr. aChangesmin «0 imethe wicinuiry 
of the magnitude 6.8 1956 San Miguel earthquake can be seen. 
Duaincgmenempervods1|945ito: 1951 MF fgourem3 no) sO Gin iehe 
Vac minty aot ethemiuturemepicenter, kl 5iglatittiudems: = li6 a0 
longitude) is low, approximately 10 compared to a long term 
average of 40, with no pronounced peak in the region. The 
wanita 4950 eto: 19568. (Bigures3 220)s showsmand istincuspeakmor 
65 in the same region. The large peak in the center of the 
region in Figure 3.20 is due to activation following the 
1952 Kern County earthquake. It should be noted that the 
catalogue is not complete for M23.5 for early years in the 
vicinity of the San Miguel earthquake. This makes comparison 


between the early values of @ and the overall average 
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unreliable. 

Figure 3.21 gives a map of @ for the last time period 
of the catalogue (1975 to 1981). The activations preceding 
and following the Owens Valley and Imperial Valley 
earthquakes are evident. Recent activation in the San 
Jacinto region, although not as pronounced, can also be 
seen. The central region of Southern California appears to 
be "quieter" than usual. Other representation using 


functionals similar to @ are possible. 


Sv. oOo Conclusions part 2 

The results of cluster analysis support the hypothesis 
that mainshocks (as defined here) form a "Quiescence > 
clustered" pattern before many of the largest earthquakes in 
Southern California. Although the number of cases presented 
is too few for definite conclusion, there is a suggestion 
that the activation stage of the quiescence + activation 
pattern is characterized by the strongest clusters in the 
region (ie, the most concentrated clusters in space and 
time). Some statistical model of regional clustering is 
needed to test this. A model that provides a quantitative 
definition of what constitutes anomalous aC GiViniy 
(activation and quiescence) for a region is needed. Cluster 
analysis will be useful in developing such a model. 

Although cluster analysis is computationally expensive 
and the results are tedious to interpret, it appears to be a 


useful tool in the analysis of earthquake Sequences. One of 
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the obvious advantages of cluster analysis over the analysis 
presented in Part lof this chapter is»that, cluster 

analysis does not require an a priori regionalization as did 
the swarms pattern (i.e. the regionalization given in Figure 


8 auler: 


3.8 Concluding remarks 

If aftershocks and foreshocks are removed, using the 
Simple definition of Table 2.1, the catalogue of residual 
events ("mainshocks") in Southern California with M24.0 
still shows a significant degree of clustering on time 
scales of 3 to 6 years and distances of the order of 100 km. 
The possibility of clustering on larger time and distance 
scales could not be tested due to the short length of the 
catalogue. The hypothesis that earthquakes are 
interconnected over large time and distance scales is 
Supported by the calculations presented here. 

Anomalous groups of the mainshocks, Swarms or 
concentrated clusters, preceded by quiescence were found to 
precede, most but not all, large events in the Southern 
California region. This quiescent > activation pattern in 
mainshocks is similar to some patterns observed for all 
earthquakes by other investigators. The time scale of the 
quiescence > activation pattern for mainshocks is much 


longer than the usual proposed patterns. 
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The separation of mainshocks from aftershocks as 
described here can be viewed as a strong filtering process 
which removes short period variations in the catalogue. 
Short term anomalies may obsure longer term variations. The 
removal of aftershocks, as defined here, from the catalogue 
made the identification of the quiescence + activation 
clearer. The distribution of events in the residual 
catalogue is closer to a Poisson distribution than is the 
distribution of events in the highly clustered catalogue 
which includes aftershocks. Deviations from Poisson behavior 
in the residual catalogue are, thus, more statistically 
Significant than in the "complete" catalogue. This technique 
may be useful elsewhere. 

To my knowledge, this work is the first attempt at 
applying the mathematical technique of cluster analysis to 
the study of earthquake sequences. The analysis presented 
here obviously leaves many gaps. It does, however, show the 
usefulness of the technique in the study of anomalous 
seismicity patterns. 

Cluster analysis is particularly adept at picking out 
regional periods of anomalous activity and quiescence. It 
appears that the technique would be very useful in the 
location and study of seismic gaps in, for example, Middle 
America. Another possible use of cluster analysis would be 
in the development of a statistical model of the clustering 
of earthquakes in space and time. Such a model requires some 


physical basis. 
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ae improvements can be made in the cluster 
analysis techique presented here. More physical metrics 
which also incorporate the earthquake magnitudes could be 
constructed. Also, a more physical means of determining the 


normalization of time with respect to distance in needed. 
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4. Bursts of seismicity near the Cocos - North American 


Caribbean triple junction 


4.1 Introduction 

Three related seismicity patterns, precursory to the 
largest earthquakes in a region, have been proposed by 
Keilis-Borok et al. (1980a,b,c). These three patterns, 
collectively termed "bursts of seismicity", consist of the 
abnormal clustering of earthquakes in time, energy and space 
before a major earthquake and were briefly described in 
Chapter 1. This chapter will apply the patterns; pattern B, 
pattern Swarms, and pattern Sigmaniro retrospective 
prediction of large earthquakes which have occurred in the 
vicinity of the Cocos - North American - Caribbean triple 
pin cdt Ol. 

Pattern B, "bursts of aftershocks", consists of a 
medium magnitude mainshock which has an anomalous number of 
aftershocks concentrated at the beginning of the aftershock 
sequence. Pattern Swarms consists of a "Swarm" of mainshocks 
where’a "Swarm" is defined aS a group of moderate events 
concentrated in space and time and occurring during a time 
interval when the overall seismicity 1S not below average. 
Pattern Sigma consists roughly of an increase in the 
cummulative seismic energy released, to the 2/3 power, ina 
sliding time window. Sigma is roughly proportional to a 


Summation of the fracture areas of the earthquake sources 
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(Keilis-Borok and Malinovskaya, 1964). Pattern Sigma is 
identified as a peak in this summation. 

In algorithmic form the patterns have been applied with 
some success to many regions, a review of which, is given in 
Section 1.6.3. Bursts of seismicity patterns are regional in 
nature. As such they do not indicate the location of the 
future earthquake precisely. In most cases the patterns take 
place at some distance (hundreds of kilometers) from the 
mainshocks which they precede. Accordingly, the occurrence 
of the patterns may be interpreted as indicating a 
Significant increase in the probability of strong earthquake 
somewhere in the region within the next several years 
(Kerlas-Borok et aly, 1980b): 

There is no adequate physical theory to explain the 
occurrence of these patterns. Various speculations suggest 
that non-linear features of "friction" on fault zones 
(Barenblatt et al. 1981) may play a role. Other suggestions 
include the obvious one that an increased density of small 
asperities (Kanamori, 1981) appears on the fault zone. The 
large separation between the patterns and the subsequent 
large event may indicate that both patterns and large event 
are symptoms of the same, as yet not understood, underlying 


Cause. 
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4.2 The data set 

For the Middle America region the only earthquake 
catalogue available in machine readable form was the NOAA 
PDE catalogue covering the years 1898 through 1979. The PDE 
listings for January, 1980 through July, 1980 were added to 
this. The nature of the NOAA data set presents two 
non-trivial problems in homogeneity, completeness of 
recorded events in time and consistency of magnitude 


estimates in time. 


4.2.1 Magnitude consistency 

Since the NOAA catalogue is compiled from many data 
sources there is the possibility of inconsistencies of 
magnitude estimates. Different sources have used different 
magnitude scales. The following scales are present in the 
NOAA catalogue: Ml (the Richter 1935 magnitude scale), Ms 
(surface wave magnitude), Mb (body wave magnitude) and Mu 
(unspecified magnitude). The most commonly used magnitude 
for statistical studies is Ms. For this reason Ms is used in 
this study. 

Unfortunately, the majority of events in the NOAA 
catalogue do not give an estimate of Ms. It is therefore 
necessary to attempt conversions of other magnitude 
estimates to the Ms scale. Bloom and Erdmann (1979) show 
that for statistical purposes Ml and Mu are roughly 
equivalent to Ms. After 1963 the most frequently recorded 


magnitude in the catalogue is Mb with no other magnitude 
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indicator given. It was necessary to make estimates of Ms 
from Mb for most earthquakes from 1963. This is complicated 
by saturation of the Mb scale relative to the Ms scale for 
larger magnitude events (Kanamori and Anderson, 1975, and 
Kanamori, 1977). Large errors in estimated Ms values are 
possible for individual events. This constitutes the single 
largest difficulty in the use of pattern Sigma which 
Critically depends on accurate magnitude estimates. 

Many authors have approached the problem of conversion 
of Mb to Ms (i.e. Geller, 1976; Bloom and Erdmann, 1979; 
Kagan and Knopoff, 1980). Here the approach taken by Kagan 
and Knopoff (1980b) is used. For the conversion of Mb to Ms 
an empirical curve of the average values of Ms for a given 
Mb was constructed using all earthquakes in the Middle 
America region which had an estimate of both Mb and Ms. The 
resulting relation is similar to those given by Kagan and 
Knopoff (1980) and Bloom and Erdmann (1979). It consists of 


two segments: 


Ms MOmiOreMDiES 26:20 
Msi =a oMb. a3, Otome Mbe ao. 0 
For the analysis the following rules were followed for 
the selection of a magnitude estimate M for an individual 
event. 
M = Ms, if Ms is present 
GieMSenOtsOLeSeniL. 


M = Ml or Mu, if Ml is not present 


If only Mb is present: 
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M is given by the above empirical relation. 


4.2.2 Completeness in Time 

The completeness of the catalogue limits the time 
window of the study. Studies of correlations between large 
events and smaller events above some magnitude M, are 
limited to the study of only those events which occur after 
the catalogue becomes relatively complete for events with 
magnitudes M and above. Table 4.1 gives the number of events 
in various magnitude ranges for 5 year intervals starting in 
1900. The table includes events in the entire Middle America 
region. Before 1910 only events with reported magnitudes 
greater than 7.5 are present. Numbers of magnitude 6.5 to 7 
events start to become significant after 1910, numbers of 
Some Omo mec Vell ouarter e930; "and numbers of 5. tows. onevents 
after 1955. The number of events with magnitudes less than 5 
increases continually. 

The patterns will be tested on the available data 
starting in 1930 with the understanding that the catalogue 
PomnGt Vet yeCOMple GuLrOLr early years a hGBLeOUMuGm LO tmeC cea, 


years could be biased due to this. 


4.3 The region 
Figure 3.1 shows the region of interest along with the 
locations of large ( M26.0 ) earthquakes listed for the area 


in the NOAA catalogue. The area is a region of complex 
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Magnitude interval 


Year COMET UME OMOR OO Oem ig Da7 20+ 
SHC) SNCS 0 0 0 0 0 0 8 
(Sholay Ene 0 0 0 0 0 0 3 
ToT OeS tS 0 0 0 0 1 é 4 
eo Sia 20 BW 0 0 0 0 Z 6 
Pela S225) 54 0 0 0 1 2 3 
(Tore TES 6 125 0 Oa 1 i 8 
ISIN OS SRC rS) 225 0 0 2 34 i s) 
183551940 249 0 0 10 8 s) 4 
1940-1945 BS) 0 oe 8 iZ 13 6 
fo 4551950 49 0 0 5 6 11 il 
fo0si So5 318 2 9 exe 34 1g, 7 
Vishay 1 eN0) 424 us 43 26m 30 18 Zz 
PS0S 965 Sch Manes) 50 qe 11 S| 3 
sis SiAe Sif Fas) | 83 oh ial S 2 
TOROS S 20.0 mars One) 30 35 ues 4 3 
Hoes VEG 33 164 OCmemaG Sie! 12 4 4 


DoDlems imei esmtNemnumberS sOteevents sli stedeing ice NOAA 
catalogue in various magnitude ranges of M for five year 
pee nvalsestarting in) 19007 Themtablceincdudes@event sm nauncs 
entire Middle America region. 
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Figure 4.1 gives an outine of the area used-in the bursts of 
seismicity study along with large earthquakes-( M > Si) 
listed in the NOAA catalogue for the area. 
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tectonics which contains segments of the boundaries between 
the Cocos and North American plates, the Cocos and Caribbean 
plates and the Caribbean and North American plates. The left 
lateral Motagua-Polochic fault zones define the Caribbean - 
North American plate boundary. The Cocos plate is subducting 
under the North American and Caribbean plates at a rate of 
8.3 and 6.7 cm/year respectively at the Middle American 
trench near the triple junction (McNally and Minster, 1981). 
The nature and location of the Cocos - North American - 
Caribbean triple junction is poorly understood. 

The plate boundary segments near the unstable triple 
junction have been the sites of numerous large shallow 
earthquakes. A list of large events ( M 2 7.1 ) which have 
occurred in the region since 1920 is given in Table 4.2. 

The destructive February 4, 1976 Ms=7.5 Motagua 
earthquake occurred on the North American - Caribbean plate 
boundary. The Middle American trench region in the vicinity 
of the proposed triple junction (near the junction of the 
Middle America trench and the western extension of the trend 
of the Motagua - Polochic fault zones) has been the site of 
several M27.0 earthquakes since the turn of the century, the 
March 29, 1973 Ms=7.3 being the most recent. However, the 
segment of the trench northeast of the triple junction has 
not been the site of a historically recorded large (Ms27.0) 
earthquake. This, the southern end of the North American - 
Cocos plate boundary, is the site of the subduction of the 


Tehuantepec Ridge. This region may not have the potential 
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YEAR MONTH DAY CAE LONG. MAG. DEPTH 
1926 2 TSO 007.00 (ies 0 
is isis Ine RL Aho “Shiai (ees 0 
1942 8 6 14.00 see U G23 0 
1944 Ge Zhe e152 00 Bd oO 0 0 
1946 6 (Ono U -94.00 4 100 
1550 1 0iieeS etl 4550 moe 1.3 0 
shes! (lah) Gee Shepste, sie)y ete 7.4 0 
1970 4 29 14.52 soem OU je 33 
1976 2 DEEL Ons 2 waehs in TEL te = 


Table 4.2 gives a list of large ( M2 7.1 ) events occurring 
in the study region since 1920. 
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for large shallow earthquakes (Kelleher and McCann, 1976). 

The segment of the Cocos - Caribbean plate boundary 
south-east of the triple junction has been the site of the 
largest earthquake in the region, the August 6, 1942 M=7.9 
event. The segment of the trench just north-west of the site 
of the 1942 earthquake was the site of the Nov. 17, 1953 
M=7.3 earthquake. Neither of these segments of the trench 
have ruptured since. 

On the basis of historical seismicity several authors 
have identified the site of the 1942 earthquake and the 
segment of the trench south-east of the 1942 earthquake as a 
seismic gap of the first kind (Kelleher et al., 1973; McCann 
et al., 1979). The lack of large events in the region since 
1942 has been clearly shown by Kelleher et al. 1973. McNally 
and Minster (1981) point out that estimated cumulative 
seismic slip in the same area has been low since the 1942 
earthquake. 

Earthquake recurrence intervals in the Middle America 
region have been found to be highly regular by several 
authors. Rikitake (1976) calculates the average return 
Baca for large shallow earthquakes in the Middle American 
region as 34.5 +3.6 years. McNally and Minster (1981) 
estimate the average return period as 32.5 +8.4 years on the 
basiuseot recorded data 1698-1979 eandeasSeso. i224. 0Syeans on 
the basis of historical data 1542-1979 in the Oaxaca region. 
On the basis of estimated average return times and the low 


amount of cumulative seismic slip McNally and Minster (1981) 
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suggest that the Guatemala coast deserves special attention 
as a zone of presently high seismic risk. 

Figure 4.1 shows that the region of highest seismicity 
in the area (at least for M26.0) is in the vicinity of the 
proposed triple junction. It is suggested here that the 
nature of the seismicity this complex region is indicative 
of the state of stress on the adjacent plate segments and 
that periods of high energy release and earthquake 
clustering in the area of the triple junction occur when one 


Or more of the boundaries is near major failure. 


4.4 Bursts of seismicity 

It is first necessary to define some terms. Strong 
events are the set of largest events in a region, those with 
magnitudes above some threshold M,. Earthquakes with 
magnitudes less than or equal to M, but greater than some 
minimum threshold m are termed the background seismicity. 
Relatively large shocks in the background seismicity, those 
with magnitudes between M, (a lower limit) and M, (an upper 
limit) are called main shocks. Aftershocks of major events 
are excluded from this classification. 

Some rules are needed to analyze the success of a given 
pattern. If a pattern is valid then the occurrences of the 
pattern should show a marked tendency to fall within some 
time window preceding major events. The length of this time 


window may depend on several parameters such as the 
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magnitude of the strong events and the region under 
consideration. The length of this time window will be termed 
t.. A success is then defined as the occurrence of a pattern 
Or a sequence of patterns which is followed by a major event 
within time t,.. A failure will be defined as a major event 
which was not preceded by a pattern within time t.. A false 
alarm is defined as an occurrence of a pattern which is not 


followed by a major event within time t,. 


4.4.1 Swarms 

A swarm is defined as a group of earthquakes 
concentrated in time and in space and occurring at at time 
when the overall seismicity in the region in not below 
average. Some additional definitions (Keilis-Borok et al., 
1980b) are needed. 

--An aftershock is defined by the following: Consider 
two earthquakes with time sequence numbers i and j, i>j. The 
second earthquake is an aftershock of the first if the 
following conditions are satisfied: The distance between 
their epicenters is less than R(Mi); the time difference 
ti-t3 < T(Mi); and Mj < Mi, where T(M) and R(M) are 
empirical functions. These thresholds must be set a priori. 

R(M)is defined as: R(M)N= 75ekini et (M)m1 Sutakenmerom 
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--N(t) = the total number of of earthquakes with M 
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occurring in the entire region within the time interval from 
GacmtOmt. 

--N(t) = a running average of N(t) from the beginning of 
the catalogue to time t. 

--n(t) is the same as N(t) but excludes aftershocks of 


strong earthquakes. 


IV 
= 


--r(t) = the maximum number of earthquakes with M 
occurring within the time interval extending from t - s tot 
and occurring within an arbitrarily located rectangle R of 
dimension Ax by Ay in latitude and longitude within the 
region (i.e. r{t) is a measure of the clustering of the 
mainshocks). The parameters s, Ax, and Ay are free 
Parameters of the algorithm. 

A group of earthquakes within a rectangle Ax by Ay is 
Galledeaiswarmiif nt) 2 oN(t) and nit) = en(t). 

The following choices were made for the free 
DabametenS:Mjje=)/..5;)'S = =2 yearns: ae=5)05) 8 ==. 5s xa Ay 
= 1°; and t, = 3 years. The values used for s, Ax and Ay are 
greater than those used in California (Keilis-Borok et al., 
1980b). These values were increased due to the larger 
maximum magnitude of earthquakes and poorer epicentral 
locations of events in Middle America relative to those in 
California. 

Since the minimum magnitude threshold of reported 
events in the catalogue changes significantly between 1950 
and 1963 two separate time segments have been considered. 


The results are given in Figure 4.2. A swarm is said to have 
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Figure 4.2 shows plots of n(t)/N(t) and r(t)/n(t) for the 
periods 1932 to 1959 (a) and 1952 to 1980 (b). Large 
earthquakes miMece7. 2a) ace indica Ledeby, vertical lines. 
Dashed lines indicate slightly smaller events of events 
occurring just outside the study region. 
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eceunrred whem bothen(h)/N(t)2170° ander(t) /n(t)> i502 The 
occurrences of swarms are marked in the figures. For the 
period of time 1932 to 1950 there are not many events in the 
catalogue. However, using the available data several swarms 
are observed. The Dec. 14, 1935 M=7.3 event occurred too 
G€arly to test the pattern. The first clear swarm occurred 
during the time window ending on January 1, 1941. This 
precedes the M=7.9 1942 event by approximately 1.5 years. 
After the 1943 large event the average number of events in 
the catalogue dropped and there were no swarms until the 
strong swarm during the time interval ending on January 1, 
1951. The swarm continued during the next 2 years. The time 
between the beginning of the swarm and the occurrence of the 
M=7.3 event in 1953 is just less than 3 years. Activity 
either remained high or the character of the catalogue 
changes after the 1953 event. Figure 4.2a indicates a false 
alarm starting in mid 1955. This may be due to activation 
before the Oct. 24, 1956 M=7.3 event which occurred just 
south of the region considered here. 

Figure 4.2b begins with the end of the swarm preceding 
the 1953 event. The false alarm in 1955 disappears as 
MUeNeAN tees Ss LessSmthan st elheresarcenoscwarins mt ron lo oomtco 
1970. The pattern fails to detect a swarm before the April 
29, 1970 M=7.3 event. Activation due to this event produces 
a false alarm (it should be noted that much of this 
activation is due to a foreshock sequence that preceded the 


strong event by a day; if the preceding time window had 
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included this foreshock activity a swarm would have been 
detected). The February 4, 1976 M=7.5 Motagua earthquake is 
preceded by a Swarm occurring during the time period ending 
January 1, 1975. Interestingly, a Swarm occurs at the end of 
the data set in the time window ending January 1, 1980. All 
Swarms occur within a degree of the intersection of the 
trend of the Motagua-Polochic fault zones and the Middle 


America Trench. 


4.4.2 Bursts of aftershocks 

According to Keilis-Borok et al. (1980a, 1980b) the 
occurrence of a "burst of aftershocks" (pattern B) indicates 
a Significant increase in the probability of a strong (M 2 
M,) earthquake in the same region within the next t, years. 
"Pattern B" consists of a mainshock in the medium magnitude 
range with an anomalous number of aftershocks concentrated 
at the beginning of the aftershock sequence. Specifically 
the pattern iS a mainshock with magnitude Mi such that M, < 
Mi < M, and bi(e) = B , where bi(e) is the number of 
aftershocks in the first e days following the mainshock and 
B is a threshold above which bi(e) is considered to be 
anomalous. 

The number of aftershocks was counted with e=2 days. 
Table 4.3 gives a listing of all mainshocks with 4 or more 
aftershocks in the 2 day period following the mainshock. 
Also listed are the large earthquakes ( M 2 7.2 ) occurring 


in the region after 1950. Unfortunately, the NOAA catalogue 
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1s not sufficiently complete to identify anomalous levels of 
aftershock activity with any reliability. No estimate of 
magnitude is provided for most aftershocks contained in the 
Catal OGuemmithewimistingein Table 46 3™includes all 
aftershocks: no lower magnitude threshold was used. 

Even with the small numbers of aftershocks the pattern 
SCOLes wel whitest. = so YearS, aS was) the case for the 
SOULNGMmccmOLntamstudyrmDe=—soeattersnocks, anc Me =—7 3 
the pattern scores 2 successes, one false alarm, and one 
failure to predict. The false alarm is generated by the 
magnitude 7.1 event which occurred on Oct. 23, 1950, three 
years and 24 days before the magnitude 7.3 event in 1953. If 
t. is relaxed slightly to include an extra month the pattern 
would score 3 successes, no failures and no false alarms. If 
B is lowered to 4 aftershocks one more alarm occurs at the 


end of the catalogue. 


4.4.3 Sigma 


Pattern Sigma consists of a peak of the function 


t 
SoC! Oeamer 
t-s 
where the summation is taken over a sliding time window of 
duration t-s to t and over the magnitude range from M, to 
M,. This pattern was first introduced by Keilis-Borok and 
Malinovskaya (1964). Specifically pattern Sigma is said to 
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which defines anomalous activity. Since § is an exponential 
Function of magnitude its successful application depends on 
consistent and accurate relative magnitude estimates over 
the entire length of time used in the study. For the NOAA 
catalogue this Condition weliemete atuali, probably! only 
occurs after 1963. Regardless of this, the pattern is tested 
here on the data from 1932. 

The values of the parameters were selected as above: 
S=c2myCaGS eM (=i svmeM, =. 0 aeMo=/ 0mm, OF 2D=. 9 leand: h=4.5. 
Figure 4.3a and 4.3b give plots of § from 1932 to 1960 and 
form 1960 to July 1980. The success of pattern § depends 
Critically on the choice of §. 

One feature apparent from a study of § is that values 
of § before 1963 (the time of significant improvement in the 
catalogue) are, in general, much larger than values after 
1963 even though the catalogue is more complete after 1963. 
This could be due to higher seismic energy release before 
1963 or to over-estimates of earthquake magnitudes before 
1963 relative to those after 1963. The latter could result 
if Ms was systematically under-estimated by the empirical 
ofan of Section 4.2.1. Even with the low quality data, 
pattern Sigma works wells fon secanlyeyeansS mle mls Ome Vent 
occurs too early to test the pattern, there is a peak before 
the 1943 event, a possible false alarm in 1946 which can be 
eliminated with an appropriate choice of §, and the largest 
peak precedes the 1953 event. After the 1953 event seismic 


energy release remains relatively high, producing one or two 
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Figure 4.3 shows DUctsmo tase lOtsL ne Bem 0d Smo 3 2et Om 960 (a) 
and 1957 to 1980 (b). Large events are indicated as on 


Figure 4.2. 
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false alarms between 1954 and 1960 depending on the choice 
i Wy 

The results after 1963 are complicated by the possible 
inconsistencies between early and later magnitude estimates. 
If the luxury of choosing separate § before and after 1963 
is allowed, the the following conclusions can be drawn:. The 
1970 event is not preceded by a peak (again activation 
indicated in the time windows following the event are in 
part due to the foreshock sequence mentioned above) and a 
weak peak precedes the 1976 event. The catalogue ends with a 


strong peak (even by pre-1963 standards). 


4.5 Concluding remarks 

Given the quality of the data set, the patterns seem to 
score well. Without including the 1936 event pattern Swarms 
scores 3 successes (1943, 1953, 1976) one failure (1970) and 
one possible false alarm (1956, excluding the false alarm 
generated by the foreshocks and aftershocks of the 1970 
event). With the most optimistic choice of § on two separate 
time intervals pattern Sigma scores the same successes and 
failures as does pattern swarms with a longer false alarm 
time in 1956. In the most optimistic case pattern B scores 3 
successes, no failures and no false alarms. All occurrences 
of the patterns are close to the proposed location of the 
unstable triple junction. Success of the patterns may be due 


to the complex nature of the tectonics in the triple 
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junction area. 

The last available data indicate the recent occurrence 
of both pattern Swarms, pattern Sigma, and possibly pattern 
B in the region due to events occurring during 1979 and 
early 1980. Based on the past history of the region and the 
lack of large events on some segments of the trench for more 
than 30 years, this region warrants attention as a possible 
Site for the occurrence of a large ( M27.3 ) shallow 
earthquake within the next 3 years. 

It should be noted, however, that this suggestion is 
based on data which is not of as high a quality as that used 
for Similar studies in California. The NOAA PDE file is 
incomplete for low magnitudes in its early years and 
magnitude estimates are, to say the least, inconsistent. 
Even the identification of the largest earthquakes in the 
region for the early years iS Subject to question. 

For example the 1953 earthquake which was identified as 
having a magnitude of 7.3 in NOAA is listed as magnitude 7.1 
by McNally and Minster (1981). This would then not qualify 
aS a strong earthquake by the criteria used here. If the 
threshold for the definition of a large event is lowered to 
M,=7.1 two other earthquakes would be added to the list of 
major events. Inclusion of these two events would improve 
the score for Sigma but not for swarms and bursts of 
aftershocks. 

If the patterns are considered only for the relatively 


reliable data after 1963 the score for both Sigma and swarms 
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1s one success, one failure and the current alarm while 
pattern B scores 2 successes. The comments presented here 
should be considered in the light of this and the fact that 
the patterns depend heavily on the inconsistent magnitude 


estimates. 
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5. Concluding Remarks 


5.1 Summary 

In this thesis several aspects of the clustering of 
earthquakes related to precursory seismicity patterns have 
been explored. Although there doesn't exist an unambiguous 
definition for aftershocks, the most obvious indications of 
earthquake clustering are mainshock - aftershock sequences. 
Using a Simple definition for an aftershock, Chapter 2 
proposes a self-similarity in this clustering, a scaling law 
for the occurrence of aftershocks. The scaling law suggests 
that number of aftershocks does not depend on the magnitude 
of the mainshock if the aftershocks are counted ina fixed 
magnitude range below that of the mainshock. This scaling 
law reduces the number of observations necessary for 
Statistical studies of aftershock sequences and imposes a 
constraint on acceptable physical models of aftershocks. It 
is also suggested that foreshocks, on the average, generate 
larger numbers of aftershocks than do mainshocks. The 
scaling law supports the hypothesis that the distribution of 
irregularities in material strengths on a fault zone shows 
no characteristic scale. 

In Chapter 3 the clustering of earthquakes on longer 
time and larger distance scales is explored. One way to look 
FOrmlOngeoaLermEClUSLelING m1 Ge Lomo bmi tc OULEEnemSHOLL 


term variations which could possibly obscure the picture. 
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This was the approach taken here. To study the possibility 
of long-term clustering "foreshocks" and "aftershocks" were 
removed from the catalogue using a simple boxcar filter 
definition for foreshocks and aftershocks. The purpose of 
the box car definitions was to remove the obvious short term 
clustering which sometimes precedes and followings a shock. 
In this study "mainshock" then refers to the largest event 
in one of these primary, box car defined clusters. 
"Aftershocks" and "foreshocks" are those events which occur 
respectively after and before the "mainshock" in the box car 
window. In this definition a “foreshock" - "mainshock" - 
"aftershock" sequence can be thought of as first order 
clustering which is perhaps imbedded in a hierarchy of 
higher order, longer term clustering. 

It was demonstrated that even if foreshocks and 
aftershocks are elimiminated, through such a strong 
filtering process, the catalogue of residual events 
("mainshocks") still shows weak but significant clustering 
characteristics. The clustering is evident on time scales of 
the order of 3 to 6 years and distance scales of the order 
of 100 kilometers for mainshocks with magnitudes greater 
than or equal to 4. The possibility of longer term 
clustering could not be studied due to the short length of 
the catalogue. 

The hypothesis that there are related events in the 
catalogue over long periods of time and large distances, 


such as is suggested by the stochastic branching models of 
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earthquake occurrence, iS not in conflict with the 
observations presented here. If, for example, the Kagan and 
Knopoff branching model, in which all shocks can be thought 
of as "children" of those shocks occurring before them in 
time, iS correct no length of box car would be sufficient to 
eliminate all dependent events. The Kagan and Knopoff model 
was discussed in Section 1.11. In such a model clustering 
would always be observed at time scales longer than the 
length of the boxcar. Unfortunately, the short length of 
available catalogues procludes the possibility of testing 
for the presence of such long term clustering. 

Anomalous periods of clustering and anti-clustering in 
the sequence of residual events were juxtaposed with the 
occurrences of large earthquakes in the same region. 
Anomalous groups of residual events, Swarms or concentrated 
clusters of "mainshocks", were found to precede, most but 
not all, large earthquakes in the Southern California 
region. The quiescence > activation pattern observed here 
for mainshocks is Similar to precursory patterns observed by 
other investigators, although the time scale is in general 
longer. 

Short term anomalies in the catalogue caused by 
aftershock sequences may obsure longer term variations. The 
removal of aftershocks, as defined here, from the catalogue 
made the identification of the quiescence + activation 
clearer. The distribution of events in the residual 


catalogue is closer to a Poisson distribution than is the 
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distribution of events in the highly clustered catalogue 
which includes aftershocks. Deviations from Poisson behavior 
in the residual catalogue are, thus, more statistically 
Significant than in the "complete" catalogue. The technique 
used here may be useful elsewhere. 

Chapter 3 also demonstrates that the mathematical 
technique of cluster analysis is useful in the study of 
anomalous seismicity patterns. It is suggested that the 
activation stage of the quiescence > activation pattern is 
characterized by the strongest clusters in a region (i.e., 
the most concentrated clusters in space and time). It is 
possible in this technique to quantify the strength of a 
cluster and eliminate the need for apriori regionalization 
in space and time. 

The characteristics of the quiescence > activation 
patterns, at least for the last three large earthquakes in 
Southern California, are superficially very similar with 
periods of quiescence Starting approximately 10 years before 
the large events and periods of activation starting 2 to 3 
years prior to the large event. The period of quiescence are 
characterized by the lowest activity in the regions. Last 
available data indicates the recent occurrence of a similar 
patternein the Sangdacinto sedi On sltwerscmsugges  cOm chong: nis 
area deserves attention as a region of presently high 
seismic risk. 

Proposed possible physical explanations for observed 


clustering (as opposed to computational algorithms) are not 
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numerous. Examples include inhomogeneities in material 
Strengths along a fault zone ("asperities") and non-linear 
dynamic friction on fault zones. These were discussed in 
section 1.11. The asperity model gives a qualitative 
explanation of the quiescence + activation pattern observed 
in Chapter. 3. Precursory patterns such as quiescence and 
Swarms may reflect the physical properties of fault zones 
rather than the physics of the failure process. 

Chapter 4 applies three proposed precursory seismic 
patterns termed "bursts of seismicity" to retrospective 
prediction of large earthquakes which have occurred in the 
VieGinityrol thescocos = North American® "Caribbean triple 
junction. The three par eence Uburstsmor®aitershocks”, 
Swarms, and sigma, are based on abnormal clustering of 
earthquakes in time, space and energy. The patterns give 
respectable results although the available data set 1S not 
Cit hemhrghestmoualiuy.. Allsoccurrences of the patterns are 
close to the proposed location of the unstable Cocos - North 
American -— Caribbean triple junction. Success of the 
patterns may be due to the complex nature of the tectonics 
in the triple junction area. There is no adequate physical 


explanation foregchesubursts Ole selsmicity — pattericr. 
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5.2 Further extensions of this work 
There are many gaps in this study. The most obvious gap 

is a lack of any physical theory which explains the obverved 
clustering. Any extension of this work should be towards a 
Viable physical model model which offers an explanation of 
the observations presented here: the self-similarity in the 
occurrence of aftershocks; the observed clustering of 
mainshocks; the quiescence > activation pattern observed 
before many large earthquakes; and the "bursts of 
seismicity" patterns. A physical model for foreshocks, 
aftershocks, and mainshocks should be sought. Occurrences of 
precursory patterns have been independently studied by many 
investigators, and are well documented. Any model of 
sequences of earthquakes which fails to explain the well 
documented patterns such as swarms and quiescence suffers a 
severe deficiency. A viable model should also have the 
characteristics of self-similarity. 

mee buSterwanalysus appears tombesamuselulitooletor the 
study of seismicity patterns. The location and study of 
seismic gaps could be approached using the technique. 
Bro cher possible use of cluster analysis would be in the 
development of a statistical model of clustering in space 
and time. A more physical means of determining the 
normalization of time with respect to distance is needed as 
well as a more physical metric for the determination of 
dissimilarity which incorporates the earthquake magnitudes. 


However, these improvements require some input from a 
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theoretical model. 


In conclusion it may be worthwhile to speculate on 


useful directions for study of the physics of the earthquake 


generating mechanism. The earthquake generating mechanism 


must in the final analysis be a thermodynamic one. Several 


constraints must apply to any model of this process. 


ee 


Most of the time an area prone to earthquakes is either 
in thermodynamic equilibrium or undergoing processes 
which are thermodynamically (but not mechanically) 
reversible. 

The thermodynamic system iS non-linear and as a 
consequence the potentials that describe it (whatever 
they are) can have more than one stationary point. 

The system can be characterized by internal and external 
Daramerers: 

Slow changes in the external parameters could change 
stable equilibrium points to unstable ones and 
transitions from unstable states to stable states can be 
earthquakes. 

In such a non-linear system cluster activity could 
ee 

If this notion is correct models that isolate elements 
of the process are unlikely to succeed. Perhaps an 


earthquake is a collective phenomenon. 


Perhaps a future research project could address the 


construction of a macroscopic thermodynamic model, with 


probabilistic properties, of the earthquake generating 






ehom tanbterceds 
alitvetiew evan oi osterioiehae i 
nt .¢4. Lp Pheee no snotzoeTi® fykeag 


a 


” ‘tno 2a fr) i=é pri IwIotey 
*. x a? aA ze hes Bil j i Sui7 iz sas os 
Viane ‘sua aintes7enes 7 


eni? wld Se gach eset —_ 
_ 


vocapads- sad & 


en iE b 
7 
Sia iis ity feel its lise sigue | : 
* 
7109 2> 2rarrt tara ; ’ 


jsopos2a8 | a. 
Pv) PDLV 1 a 45 ss ¢ve@ deeail~-naue s oowe nt s€ a 

23 S08 2 /- 

: Spayests: = sesva: vert tisbe® Joeuteo al. colada uM Lica it~ il 
7 ity apelin’ .beenwe co: viedilow or san agigy wad) 7 








“process. a 
+ 





Abt PeQu 
€ v? 
- % 
ral 
/ 
‘4 
8 » 
a. » 
f ’ “4 7 Se 
o ‘ £ ; + ws + €<s! ‘ 


poe. "Ge Soo ax « [cua Divi. &e 
“50iuce ~ Vine 1¢)es lori: Jaethe 
Gece £( oy asl}, Seta oe 
Ea ~ 


: - 








184 


References 


Abramowitz, M., and I. A. Segun (1968). Handbook of 
Mathmat ical Functions, Dover Publications, New York. 


Acharya, H. K. (1979). A method to determine the duration of 


quiessence in a seismic gap; Geophys. Res. Lett. 6, 
681-684. 


Aki, K. (1979). Characterization of barriers on an 
earthquake fault; JU. Geophys. Res. 84, 6141-6148. 


Aelen a Creme GlO75)7mGeologicalscriteria stor evaluating 
Se1SsmiCuly moll GeO il woOG ee AINer 66, 104411057. 


Alen ae CemmeR eee oes CMAMONC  sGarr urRilchter, and J. .M. 
Nordquist (1965). Relationship between seismicity 
and geological structure in the southern California 
GEequon weet | /eeSe lS] MSOC AMNeIT e855, 8/53a79 7% 


Anderson, J. G. (1979). Estimating the seismicity from 
geological structure for seismic-risk studies; Bul]. 
Sei smoln wooGarAniens 69 74135-1587 


Ando, M. (1975). Source mechanism and tectonic significance 
of historical earthquakes along the Nankai trough, 
Japan; Tectonophysics 27, 119-140. 


BarenDlatt mec el ekVie IKentSSBOLOK  yanceM eam wava Sit k 
(1981). A model of the clustering of earthquakes; 
Geophys. J. Roy. Astron. Soc. (in press). 


Bloons hieD eeeandun. 8s BE dina nim el S70) Mee uence 
magnitude - time relations in the NGSDC earthquake 
data, file- Bul J. Selsmola .SoC. Amer we 69 2 0S5 2090. 


Borovik, N., L. Misharina, and A. Treskov (1971). On the 
possibility of strong earthquakes in Pribakalia in 
thee hutuce: micU«eiGa «6 5G laUoonerhVS: @oOlMG leant |, 
isa Gr 


sor0ne1e 244 






et Chew , .1 608 . i « CSE@RRRTER, 
diaf « Ph -ii ‘avatl | Ser Foes! esr See 


jee a yA ,(2Te.) © A (eee 
aa | ‘“éscee a Al gonae@eige 
ia-T oa 


AS Sea's 20 (every, shea fee 
fLival AAKUPHIIES i 7 7 


oftye bee (ere! t.. ,49e7sBnA 
in “ROT fhitaicvs Llesrggss 


i! 
1 te | at 


i 


, - ~~ ty 2 r & i 5 aM " IBIGK 
tea fanbeoseid de | - 
5 fi + r rs } f Gye iy 


a so 2 r ari e An et os cA [ {ea * uv | a 4 +teldasseH _ 
tdavadoisr24 %- mitateuvic ets. .9o tetom &)) ftSers 7 


enh). soe .foarek OR ok .eleoee 7 


‘ z ina } eG. oa 
staat conaiae, Tina eet tat 
ne ved 
. = 





185 


Brady, Be T.)(1974)eeTheory of earthquakes = 1.’ A scale 
independent theory of rock failure; Pure and Appi. 
Geophys. 112, 701-725. 


Brady, B. T. (1975). Theory of earthquakes - II. Inclusion 
theory of crustal earthquakes; Pure and Appl. 
Geophys. 113, 149-168. 


Brady, Beedle 97 Gavretheony ofeearthquakes - 711 In 9Inelusion 
theory of deep earthquakes; Pure and Appi. Geophys. 
Las 2a 


Brady, B. T. (1976b). Theory of earthquakes - IV. General 
implications for earthquake prediction; Pure and 
Appl. Geophys. 114, 1031-1073. 


Brady; Bem mlovoc)emDynamicsaoteraultmqrowth = sAmphysical 
basis for aftershock sequences; Pure and Appi. 
Geophys. 114, 727-739. 


Brune, J. N. (1968). Seismic moment, seismicity and rate of 
slip along major fault zones; U. Geophys. Res. 73, 
117-784 


Bute CeaLGEPERED Wea cshyaands Rem OssBunrOrd «i S77)™ 
Steady-state seismic slip: a precise recurrence 
model; Geophys. Res. Lett. 4, 91-94. 


Bute we HGe wommeH ee hl Uuke mand Resi WeSsSonmil 7 4)) 
Premonitory vertical migration of micro-earthquakes 
in Central California - evidence of dilatancy 
biasing; Geophys. Res. Lett. 1, 221-224. 


Byerlee, J. (1978). A review of rock mechanics studies in 
the United States pertinent to earthquake 
prediction; Pure and Appl. Geophys. 116, 586-602. 


Caputo), Me, PB. Gasperiniy )V. skeilis-Borok) @U meMarceliiy sand 
I. Rotvain (1977). Earthquake swarms as foreunners 
of strong earthquakes in Italy; Annali di Geofisica, 
MX XKONS or eec OS oO. 


Chinnery, M. A. (1969). Barthquake magnitude and source 
parameters; Bul/. Seismol. Soc. Amer. 59, 1969-1982. 


€2! - ; 






wisse A .D+ get coupe eae 10 eae (ACE: ee 
=i 2a isnivijgal é20009o eer cosen 
| greror 5). eee 


(ca! ati ~amee . (2eeT) Lal pe 
ar Lt cyselt sieateld v2 oe vee ) 7 
7 


Gos “oT euiede6e Geeh tn ~toes ) 
rot OT ; 


A: Gia 
Noo 
aj 7 
siee , ide » 4sau rT ; 
y om f ; 7 2 - 
b u 
, 
"a é =) - 4 > , alia 7 
& J2 a 
; on ; 
- 
- 


6 BALIN A ie eee eee sive 
. r 3e3 tac0mess 
(ied Lattdeda : 


WW) .@ynese toateaia , 25 


a7 Nd) BG.Wy co © ci teytus oO: tp wetvet A .(BV8!) vb eee 

_ TUivaSA Tks SA sali seg ketasa Ses lay edz _ 

| Mats: F*  ckosd Aqeh tie et poiesiieag! 0 
7 r. oY _ 


a a, . 
v4 ; 7 ; 


: 7 Pa = i ‘a | A ba 


Pes my * iota - 





’ 
; : 


186 


Cramer, H. (1946). Mathematical Methods of Statistics, 
Princeton University Press, Princeton. 


Das, S. meanders Akiw(1977a)e0A numenical study “of 
two-dimensional spontaneous rupture propagtion; 
Geophys. JU. Roy. Astron. Soc. 50, 643-668. 


Das, S. , mandek eAkimal977Zb)earaulteplane with ebarriers: a 
versatile earthquake model; JU. Geophys. Res. 82, 
sietanel = Syey7/ i}, 


Davies, G. F., and J. N. Brune (1971). Regional and.global 
fault slip rates from seismicity; Nature 229, 
Ose Cees 


Delsemme, J., and A. T. Smith (1979). Spectral analysis of 
earthquake migration in South America; Pure and 
Appl. Geophys. 117, 1258-1285. 


DODGOVO Sy ae emeD ees eles BOKOV anda eels Milachkins (1979) 7 
Estimation of the size of earthquake preparation 
zones; Pure and Appl. Geophys. 117, 1025-1044. 


Duda, S. (1965). Secular seismic release in the 
circum-Pacific belt; Tectonophysics 2, 409-452. 


Duran, B. S., and P. L. Odell (1974). Cluster Analysis A 
Survey, Springer-Verlag, New York. 


Hz rewonsky Aw alGuA «1G. PrOZOLOVm 961 ee SpacCeauime 
clustering of earthquakes; (in press). 


Ellsworth sew. Loy sAstGealbindh,. Weg ePLeSCOCi, sande aac. 
Herd (1981). The 1906 San Fransisco earthquake and 
the seismic cycle; (in press). 


Engdahl, E. R., and C. Kesslinger (1977). Seismological 
precursors to a magnitude 5 earthquake in the 
Central Aleutian islands: *Joun) Phys) ‘Earth 25; 
S2ase5200K 


Evernden, J. F. (1970). Study of regional seismicity and 
associated problems; Bul]. Seismol. Soc. Amer. 60, 


‘yi Tel ate 
qi at 
a+ 
=) 
» g : ee : ale By = 
Se. Se cigterses oor 









atest esi ecpartsun: . (pedis “a 
i sears eo rene) at ness aaa 


4 | ayn f = A = ey ' ind Al Soe ve : leeal 
y) openness el he sim 1 -veeeet ) 
Sy Get,” ee ee .).2yigos 
or a . 1c Dat on bas ,a@ ,aed ¢ 
y La} i e.,7t8.9¥ 
tce~ gga 
; . fod » ot 2 Gaetend a 
; * sige? - 
4 »7 ¢rf ; ; 
2 .? A bee .20 . Smiteaiad : 
S1gim otavpti ftes 
> Tie ee : \OGr 
: 
oa sarin « * , 
n i 2,4.) pydeievendes a 
witguigead | 
‘v1. 4 Baas 
" 
~a2 . (dee?) .& (abe 
i | ve & ew = er! > 
: bia 2 Bae 
“ i {3 ; a", us =“ 
- 
; 7: 
: . A ,yaenove. ae ¢. 
, iu mri 
~~ 
Ta 
7 diet © 24.) <0 <8 dpe : 
esi cad SUP eat .( ge eae 


»  ‘hH0%e 


ola @paes 


ao ae 
ae 


pat} 


187 
393-446. 


Evernden, J. F. (1971). Location capability of various 
seismic networks; Bul]. Seismol. Soc. Amer. 61, 
24i=2565 


Evison, F. F. (1977a). Fluctuations of seismicity before 
major earthquakes; Nature 266, 710-712. 


Evison, F. F. (1977b). The precursory earthquake swarm; 
Phystebanth eel anet. Tnteriorsy 15, =pl9=p23% 


Evison, F. F. (1977c). Precursory seismic sequences in New 
Zealand; N. Z. Joun. of Geol. and Geophys. 20, 
P29 —s4 1 


Evison, Ff. 7F.9(1977c). Long-term seismic ‘precursor to the 
1968 Inangahua earthquake, New Zealand; N. Z. Joun. 
of Geol. and Geophys. 21, 531-534. 


Fedotov, S. A. (1965). Regularities of the distribution of 
strong earthquakes in Kamchatka, the Kurile Islands 
and northeastern Japan: Tran. Acad. of Sci. USSR, 
DOStTRaPhVS Mab ar byes 6)7966-93% 


HeCOLOV OCIA NE EAM eDOLD Tl kandem Vi gN@ ae MOLOZOV san sl. 
MyachkingsV.88.) Prebrazensky,; sand) G.. A. sobovey 
(1970). Investigation on earthquake prediction in 
Kamchatka; Tectonophysics 9, 249-258. 


ErLediian Mesias WhitCOMmbD, mC. shea ely anGmurmeAr 
Hileman (1976). Seismicity of the southern 
California region, 1 January 1972 to 31 December 
1974, Contribution No. 2734, Division of Geological 
and Planetary Sciences, California Institute of 
Technology. 


Huis iG eC manGeAtmE Gl Lindhi@!07.o) PAmchalGcutimerat tie hone 
orientation between foreshocks and aftershocks of 
the Galway Lake earthquake, M=5.2, 1975, Mojave 
Desert» Calitornias(abstract) 3 hectonopnys! Cs 52, 
GOu=602% 


Gardnens JaiK.) andebsknopoff (1974)%) Is the sequence ‘of 


are~Zee 





o =, 4 F Ye } | a nos tvea ; 


, vy : =" noe iva 


vt" 
7 al iys 
i fi * _ aa 
, Neale 
Le 7 
: i o - : +" 
Ser) a 428 \vogobel | 
' tore hie 
5 vi : : Siri 
1 o } i | obeg 
i» J vie 
j er) 
* iIe2h Deas 





Rite Bip 6 aBl A 92,0 moocidli » .% 4, BM pheneiaee 
32 3. ~ yYsinieagied . (ate) nage | 

Ae veeunes |. nOlos sins tee 
Re mele 7lG. Cals .64, moisadiasao3., at 
2550 gicvaliied »accesing: y1e8 ae 














Lie 


ae 





4s 
— 









188 


earthquakes in southern California, with aftershocks 
removed Poissonian?; Bull]. Seismo]. Soc. Amer. 64, 
hei eysl— ile 


GaSpe minis h aM Caputo ABV. = leeKel laus~Borok mG: Marcelli, 
and I. M. Rotvain (1978). A swarms of weak tremors 
as precursors of strong earthquakes in Italy, in 
Vichislitelnaya seismologiya, (Computational 
seismology), vol. 11, Nauka Press, Moscow, 3-13. 


Geltand, =. eM. mOl.mL eGubermanyaVieslemkelliSs=sBorok, cand E. 
Ja. Ranzman (1972). Criteria of high seismicity, 
determined by pattern recognition; Tectonophysics 
1S. CHGane | 


Gelfand ale aes rele GUDEr Man eM -meomaKalezkajya,iy. 1. 
Ker lus -BOLoOkeh Ua smranzmaneel.  M. Rotvain- and) M: 
P. Zhidkov (1974). Recognition of places where 
strong earthquakes may occur, II: Four regions of 
Asia Minor and S-E Europa, in Vichislitelnaya 
seismologiya, (Computational seismology), vol. 7, 
Nauka Press, Moscow, 54-71. 


Geltvand) =i aMe Oh tmAGmGUDErMan, OV cml KealissBonok, gb. 
Knopoff, F. Press, E. Ya. Ranzman, I. M. Rotwain, 
and A. M. Sadovsky (1976). Pattern Recognition 
applied to earthquake epicenters in California; 
Phys. Earth Planet. Interiors 11, 227-283. 


Geller hind MMO CNR Ee SCalingmrelaclonmroumear enguakes Source 
parameters and magnitudes; Bul]. Seismol. Soc. Amer. 
66n ae DUNE iio 237. 


GorshnkovieAcmi Mim CaDULO, BVem lt mhetli SSBOLOKmana N. 
Ofitserova, E Ya. Rantsman, and I. M. Rotvain 
(1979). Pattern recognition applied to sites of 
possible occurrence of strong earthquakes. IX Italy, 
M > 6.0, in Vichislitelnaya seismologiya, 
(Computational seismology), vol. 12, Nauka Press, 
Moscow, 3-1/7. 


GutenbercwmB mondeGe HemRichLer smO54) es Seismicity or the 
Earth and Associated Phenomena, Princeton University 
Press, New Jersey. 


Hasegawa, A., T. Hasegawa, and S. Hori (1975). Premonitory 


sas # 7 


nS tides oc he Sia ante iba ini Set AL. ced abodes 
= a ea 1 5 ( i@eqmegicl Ong 
Vat * feat 





ht oe 9 bed tegie 
ron «5 { bv 
ee be ty af : 
| inig@aly | ’ 


t 148 \{ 37eu 8 ie 4 : : 
M onSjLeq 7 
Parks at 7 
“eT ciTe rep 
ars am | ; 
7 


— 
—_ 
= 
—_ 
rh 
ra 
a 


ao: 


P 7 a & 
“a “tf 
wah i 7 
‘ i FF ’ 
in te? bee) 
’ 1 foruhan 
, eo iel-at tf has 
hi » el Pa) r ly 
ia i] 


r aa 5 Sele 334 ‘ : ese! 4., io vosudaved 
biysusesM 1 & ta > oes Y 4. ates easis0 ; 
at a4 2 ff J r pr? ‘ i? , | 
‘Ohad’ ni opsoupea- Ot. ves ae 20h 3 22 OS a 42 eacg 7 
A rts Bay Valea vol «Os Rr 
~~ aici “St ity .) yor'cittaa teil dy duciee4 
Li % - a. 









eels) 


Variation in seismic velocity related to the 
Southeastern Akita earthquake of 1970; Joun. Phys. 
Earcthm@2orm@ 1092203" 


HilemanymulsA tA Ce ReeALLen manded.M, Nordquist (1974) % 
Seismicity of the southern California region, 1 
January 1932 to 31 December 1972, Contribution No. 
2385, Division of Geological and Planetary Sciences, 
California Institute of Technology. 


Holmes, R. A., and R. B. MacDonald (1969). The physical 
basis of systems design for remote sensing in 
agriculture; Proceedings IEEE 57, 629-639. 


Imamura, A. (1928). On the seismic activity of central 
Japan; Jap. JU. Astron. Geophys. 6, 119-137. 


Inouye, W. (1965). On the seismicity in the epicentral 
region and its neighborhood before the Niigata 
earthquake; Kenshin Jiho 29, 31-36. 


Ishida, M., and H., Kanamori (1977). The spatio-temporal 
variation of seismicity before the 1971 San Fernando 
earthquake, California; Geophys. Res. Lett. 4, 
345-346. 


Ishida, Ma,sandeH) Kanamori@(1978). \Thertoreshockvactivity 
of the 1971 San Fernando earthquake, California; 
BUTT Sel STO MN SOCHeAmMes. 68> 126551279. 


Ishida, M., and H. Kanamori (1980). Temporal variation of 
seismicity and spectrum of small earthquakes 
preceding the 1952 Kern County, California 
earthquake; Bul]. Seismol. Soc. Amer. 70, 509-527. 


Ito, K. (1980). Periodicity and chaos in great earthquake 
occurrence? U. Geophys. Res. 85, 1399-1408. 


Jardine, «Ne; and Re Sibson y(197/ 1). Mathematical taxonomy, 
John Wiley & Sons Ltd., London. 


Jones, L. M., and P. Molnar (1976). Frequency of foreshocks; 
Nature 262, 677-679. 





it BS inh fev Geer iee fi 
s iHyoans tae avg - La Ale 


£Os-e5t 


no 12 gaa? 
‘antiguel 
52 ARs 


Petit’ .n ook 
a re “ye tisleete2 
ra ' Luts yo ; 4 ol | i a 7 ' @ “Ine 


Livia .aees 
t 7+ 10 ed Ae 
* "5 . oy 
= } : \ e° att 
eS 
a 

tL) oe 

if 3 
i ! ee OL 
% B O35 
5 > Wf 

a 
; 3 

or 

é - 
, 4 
or 9 Ph ‘\ 2 40 
. * t ‘ 
+4 as, % i ‘i 
i, (onc? ré 4 ( i. | 
7 ; ; : Y2s-5 , y os 
7 ey i4 iz 7 75 ons 12 Dimer $e 
4 a 7 a ny ‘ 7 te ie | 4 at Aft ‘ f | er | 10 
iy fj ‘J p if 4 ‘ 

a-toer OX i y ie | Se oe ,; é@lavodises 


-{0Ge 0 
aWEyMA) ks 1 STGP LIVIS 
we 


PIWipioi=s inert RG y th 20S tbo 0799 
Mine beet As ean 











J. ,ranelig 


BTUMemi 


,4adistel 


shi nei: 


a sat 








190 


JOnecwmua Maj, ance. Molnar9@1979) ‘Sometcharactenistics of 
foreshocks and their possible relationship to 
earthquake prediction and premonitory slip on 
faults; JU. Geophys. Res. 84, 3596-3608. 


Kagan, Y., and L. Knopoff (1976). Statistical search for 
non-random features of the seismicity of strong 
earthquakes; Phys. Earth Planet. Interiors 12, 
Pe) ew Mh ote 


Kaganpesjeand mo eknopoft 60197) saBbauthquake risk prediction 
as a stochastic process; Phys. Earth Planet. 
interions 14475 .97-108.. 


Kagan, ve neand .beorknopolimulo ie iasStatusticaliistudy of the 
occurrence of shallow earthquakes; Geophys. uJ. Roy. 
AS CRON OCHIED > RO a0 0- 


Kagan eo nvoem mhnOpO mie | 9d0a jnmopatraledistribution of 
earthquakes: the two-point correlation function; 
GEODAVS Si MIROY SAS ERON: 25001862; £30353 20% 


Kagan, Y., and-&. Knopoff (1980b). Dependence of seismicity 
onmdeptn: sculls, Sersmol: Socm.Amer 78/0), Vein a1822.. 


Kanamori, H. (1977a). Seismic and aseismic slip along 
subduction zones and their tectonic implications, in 
Island Arcs, Deep Sea Trenches and Back-Arc Basins, 
AGU, Washington, 163-174. 


Kanamori, H. (1977b). The energy release in great 
earthquakes; U. Geophys. Res. 82, 2981-2987. 


Kanamori, H. (1981). The nature of seismicity patterns 
before large earthquakes; Maurice Ewing Series, 4, 
AGU (in press). 


Kanamori, H., and K. Abe (1979). Reevaluation of the 
turn-of-the-century seismicity peak; JU. Geophys. 
Res. 64 Gisi=6 139. 


Kanamori, H., and D. L. Anderson (1975). Theoretical basis 
of some empirical relations in seismology; Bull. 
Sai snolteesocw Amen. 65, m1 0us—1095:. 






a *> yaes eadde sent Te) ) gaqtos A cine os yd 
f vty nae : . reo £266 Sae S35 pert wie” 

i Vet hae st id MIS Sshest Clave 
Aes Seek 1 sor et Re ‘uy tesion? 


is f=i3a3 7 ‘Sent «1 Diag , aaa 
4 >i > eT Pile: WaT roe _ i 
; “ty. 19 7) AS er AP. iné= Tes 


‘ ¢°¢ 4 T 


- - 
vi ce pes ; 
7 a> @ a 
Pihat: kant t : 
i i f .. ~aapan 
a wi) be be | { 
oy v\ Te 
i 7 ; See aa 
4 "14% 
bee 4 i @ "Ex! 
fi 7 7 o-% a 
= 4) 
' M4 tvomenga 
% niue " 
_ | 
} + 1% ~ ; L © bk 
f raéy 4 iA - 
(Srenk a TTC) «H | AT 
Tis ay .2CCE? eo Geos bupAs IAS 


‘ee seikg ig 33 wey ’i 14 wtugon of? ,/ °8EI) oh  fsehenen : 
De eit a's oe exisupdsias aprai #34 sa 
7 | ivy age 


cer 





1a 


Kanamori, H., and G. S. Stewart (1978). Seismological 
aspects of the Guatemala earthquake of February 21, 
1976; JU. Geophys. Res. 83, 3427-3434. 


Kant orovivch alia vieetGleM aeMOUchan, eve) lakers zBorok, and 
Eee Vee Viilwkovich (1970 iawAestatisticalemodel fof 
Seismicity and an an eStimate of the baSic seismic 
ehLeC ts mie lam Ca ae oC UU SoOneEnVS 50) J0 seanrtiy.5, 
852510 2 


Kasahara, K. (1973). Earthquake fault studies in Japan; 
Haitian (Wetelits tks See, (Meigen (0h, (tear 


Kellis-Bororwe Viesliepiv «9M. ePodgayetskaya, and A. G. Prosorov 
G19 7) Onetherlocalestatistics of catalogs of 
earthquakes, in Vichislitelnaya seismologiya, 
(Computational seismology), vol. 5, Nauka Press, 
MOSCOWR MD D= 7.9 


Kevla sasOuOk vi melee eK OpOlt mange. 9M. ROtVvain (198025) 
Bursts of aftershocks, long-term precursors of 
strong earthquakes; Nature 283, 259-263. 


Ke isls-BOnOk mv else KN ODOL te meleemM, eROCVa1N and s1).8M. 
Sidorenko (1980b). Bursts of seismicity as long-term 
precursors of strong earthquakes; JU. Geophys. Res. 
SiS enous ao les 


Kewl s-SOuOM = aVEElE Ll KNODOL france Cemen, MALLenaM So OUCHE 
Long-term premonitory seismicity patterns in Tibet 
and the Himalayas; JU. Geophys. Res. 85, 813-820. 


Keilis-Borok, V. I., and L. N. Malinovskaya (1964). One 
regularity in the occurrence of strong earthquakes; 
JU. Geophys. Res. 69, 3019-3024. 


Keilis-Borok, V. I., V. M. Podgayetskoya, and A. G. 
Prozoroff (1971). On the local statistics of 
catalogs of earthquakes; Computational Seismology 5, 
Sil TAS) 


Kevlis-sOLok av mie cndmlameM.e RotvVai ge 978)h5.0n 
soon-to-break seismic gaps, in Proc. Conf. VI: 
Methodology for-Identifying Seismic Gaps and 
Soon-to-break Gaps, U. S. Geol. Surv., Open-File 







inde: 2 AiG call ie 


“twiedties <fegeene mee aye we 


Px Pa ‘ LY a) iD ie Sr.) Ns ay ee b 


Reis) 4 o eel’ ot (ApiveIesten 
phat t seek it 
iS Mie IGA YIsSRRee 
a c): sessetie 


, oat “88 


z 


ijod@ . (eter) sh jevedeeen 
. , sy } (ie 


a? Gir ASLOE “alten 
‘- 
9 ‘Sup tee © | 
. i a 
: 4 MW 


i-@ l,i ison 


& a 
KS [ 
: vie} E 
“)) 3”"Sid ta 
4 » a os oe r ‘9° 
{ . x 5" " - 
~ é S . = { 
4 ti a 


7 ‘ ‘ 'g : “Of ‘ _s ' v i. 6-27; 11268 ~ 


7 Ap ‘Si Ss) eee iL VV ,AaSveB=ai isan 


ag ni 2iseloees 
—ss ~ sr’. —., \ 
3? s488F ,cyisoe 4. 


‘eo .otatacoyephgt..M .v ..1. .¥ <aeeaee 
tea4 iagzol e424 oO 435004 Ho voaeas 
ach aeons : 





12 
REpOL Ee Om74578s57-594" 


Kelis BOGOK aan le emrand lem M-mROLVaInsi41079), .TWwo 
long-range precursors of strong earthquakes, in 
Vichislitelnaya seismologiya, (Computational 
seismology), vol. 12, Nauka Press, Moscow, 18-27. 


Kelleher, J. A. (1970). Space-time seismicity of the 
Alaska-Aleutian seismic zone; JU. Geophys. Res. 75, 
br 4 5e5756" 


Kelleher, J. A. (1972). Rupture zones of large South 
American earthquakes and some predictions; Ju. 
Geophys. Res. 77, 2087-2103. 


Kelleher, J. A., and W. McCann (1976). Bouyant zones, great 
earthquakes, and unstable boundaries of subduction; 
J. Geophys. Res. 81, 4885-4896. 


Kelleher, J. A., and W. McCann (1977). Bathymetric highs and 
the development of convergent plate boundaries, in 
Island Arcs, Deep Sea Trenches, and Back-Arc Basins, 
AGU, Washington, 99-114. 


Kellehermpeud sen leroy Kespeanded. sOliver (19/3). Possible 
criteria for predicting earthquake locations and 
their application to major plate boundaries of the 
Pacific and Caribbean; U. Geophys. Res. 78, 

2557 22505), 


Kelleherpard A yerandaeeSavi Om 1.97.5) emma S ia OliLo hao L 
seismicity before large strike slip and thrust-type 
earthquakes; JU. Geophys. Res. 80, 260-271. 


Kellehets, ds Aw, 0.) Savino, 1. ROWLGELe mands WemMc Gain 
(1974). Why and where great thrust earthquakes occur 
along island arcs; JU. Geophys. Res. 79, 4889-4899. 


Kelleher lene nov Kes weandaulnmoliverm(ho/3)emrossible 
criteria for predicting earthquake locations and 
their application to major plate boundaries of the 
Pacific and Caribbean; JU. Geophys. Res. 78, 

25 ag a2 oon 


ei 


. idan 


OCT THE | ANe-00 2aeqee) 


(efor 1PSey raas 
: 


‘yas grore® 3¢ 
wy) .e i 
7 - + i > 
Pa wv 
) 
; s [ 
e,% 
‘ 
J 
by ~/. 
f | > 7 
a 
ror? rT i? 
cat wit — hiv = 


hd own ie: »wor 


tee fol in aa fee 


a) 


1. Peat yal a¥ derodmak fell 


| S20RITVGeTs srrats ack 


| ,A 40 4 tee tier 


Nai tn ever tat! Gal iN 
ST, ue 45 OE 
weise gaoituels-spsuaie 


eevle-dbVe 


i aes ee 


a He ee PPS 18s 
7% _ i‘ YEXSa) 


ofa 76s 


> nm : 

‘a . 1 ia] oo) (4 
“zie | 
rel ' us 

p ' 2 7 ‘7 be 
? 7a 
j 
oR 
1200 £ inas 
rloe 


7 i ved siti “ws sl ‘ Lia 
tosei S203 7 24nigelse 
ayrnacdes., aw 8-6¢ sold 720 


oH Le 18) ; vag sh 4s A _ , isdn 


o3ede Sho. ee We 


+o )2220 Gualel gagie 


a 
= 


. vacdol ied 


 rodetise 


“sielieg 










1938 


Khattri, K., and M. Wyss (1978). Precursory variation of 
seismicity rate in the Assam area, India; Geology 6, 
685-688. 


Knopofifsel, m97 1) eA ustochasticemode! tforathesoecurrence of 
Main-Sequence earthquakes; Rev. Geophys. Space Phys. 
Cee i 75a 108% 


Knopoftt, ih. wandsYs Khagane (1977 )-mAnalysuswotsthe theory tot 
extremes as applied to earthquake problems; Ju. 
Geophys. Res. 82, 5647-5657. 


ROSobOkOVeEVEN anol seROtValnm(a 97 7) -mePattern MecOgniLion of 
earthquake prone areas, VI. magnitude M 2 7, in 
Vichislitelnaya seismologiya, (Computat ional 
seismology), vol. 10, Nauka Press, Moscow, 88-105. 


Kristy, M. J., and D. W. Simpson (1980). Seismicity changes 
preceding two recent central Asian earthquakes; Ju. 
Geophys. Res. 85, 4829-4837. 


Landgrebe, D. A., and LARS Staff (1969). LARSYAA, a 
processing system for airborne earth resource data, 
LARS Information Note 091968, Purdue University, 
Layfayette, Indiana. 


landgrebe, DamAmandelem.>. Phiniips (1967). A multichannel 
image data handling system for agriculture remote 
sensing, in Proc. Seminar on Computerized Image 
Handling Techniques, Washington, D. C., XIT 1-10. 


Lay, T., and H. Kanamori (1980). Earthquake doublets in the 
Solomon Islands; Phys. Earth Planet. Interiors 21, 
26 3=504, 


Lukk, A., and V. Ponomarev (1972). Trends in the variation 
of the seismic background in the course of time; 
Tous AcaesSCiPeUSSREPAVS  SOlIdSEAaCE 5S7m3y 2 12270, 


McEvilly, T., W. Bakun, and K. Casady (1967). The Parkfield, 
California earthquake of 1966 Cocos plate; Bull. 
Seismol, Soc: sAmer.*57, 51221. 


McCann, W. R., S. P. Nishenko, L. R. Sykes, and J. Kraus 














. Me tasydg7t ,f STC Gaye Mm baw oe 
ary OD soot dev (mies G00 28 Goat. Viste 


- S0as9ed- 


1 Oe FEtSSI4 i VOT? of 4 tigeore 
2 ““ 94 gS one 2c hapa aia 


v3, «89! = : a 


e nett ’ png wah gos tcaguaa . uy . 


ka 85 : l 9 IAS 
' ‘ a .~ 2) : 
ne «+ Vo ehodgeéor 
J POT: % 7He 
7 i ¥ 
. se .Vreraa 
Y =e iS 
* est) 
a2 ’ et4 yet - 
= » ‘. 
é > » Ey a pat: Pave 
LJ af 4 Ss | 
z wa. a i= 20hneg 
’ t 1% Pa Le 54 a 
on 
i. C , Ur +4 
» be : 
; jacana Mt Oi se 
| Ot Sa) ie ye avrag 2" :abonlaei> son@iaa ; 


+ 0¢98ee 





fe ni sAnenT ~LSTGti ve ioanonod, . Vv Sen 40h ene 
ae of. ot. Sagospdoed oimed oo ate Seaou4 
met A583 Dhjoe Joy ear a we 
Bhs ech BAM. tates” A i 
- : c>. “ao ah — 























a 









-; 






194 


(1979). Seismic gaps and plate tectonic: Seismic 
potential for major boundaries; Pure and Appl. 
Geophys. 117, 1082-1147. 


McNally, K. C. (1977). Patterns of earthquake clustering 
preceding moderate earthquakes, central and southern 
GaliiicuiiameabSeraGte -mrUSe DO; ingo. 


McNally, K. C. (1981). Plate subduction and prediction of 
earthquakes along the Middle America trench; (in 
press). 


McNally ,@keCeanded, BY Minster (1981). Non-uniform 
seismic slip rates along the Middle America Trench; 
JU. Geophys. Res. 86, 4949-4959. 


Minster ed eB e wanders Hemoordanms e976). Present Dayaplate 
motions: U. Geophys. Res. 83, 5331-5354. 


Mizone, M., M. Nakamura, Y. Ishiketa, and N. Seta (1978). 
Earthquake prediction from micro-earthquake 
observation in the vicinity of Wakayama city, 
Northwestern part of the Kii peninsula, central 
Japan) coun. Phys. Earth 26), 397-416. 


Mogi, K. (1968a). Some features of recent seismic activity 
invandenear wWapan, 91s Sui) es banthquake Res. Insc. 
WOKVOSL6 ye L225 ai2 36% 


Mogi, K. (1968b). Sequential occurrences of recent great 
earthquakes; Joun. Phys. Earth 16, 30-36. 


Mogie K.n(196ec je Migrationpof Selsmic acuiviby- bUie 
Earthquake Res. Inst. Tokyo 46, 53-74. 


Mogi, K. (1969a). Some features of recent seismic activity 
jneand neapedapan,e2, PACtIVity sbelore sandwatter 
great earthquakes; Bul]. Earthquake Res. Inst. Tokyo 
AT eo Daa leis. 


Mogi, K. (1969b). Relationship between the occurrence of 
great earthquakes and tectonic structure; Bull. 
Earthquake Res. Inst. Tokyo 47, 429-451. 


$¢! 














| ‘ 
.. Sa ee" ; ee oe ee 
rei a 
i 
f Pn f 
1° DO .D > eal 5 We ee 
% 
(A ostantiaqiiss 
tee ia 
589 oad, A Vie 
; 0 \ijoe ise 
ae 
28024 708 
s7yhe “ ,eobgee 
\, 4 i aad y Op at 
. :. , Fo 
‘ wn ] 4 
i t 7 < ‘ lay ay 
: . 2 > ae 
' Lae o> hy at °A it a 
, \ a a nf oh 
<s 
supse .dGeah), ail eee 
: 2 b j iovntsyes 7 ? 


(i est : Stren } weiter . Gheeea (ton 
; if os , Oe  Swaecs) as Sy “+t Zywignes ‘a 










pores a 


: wiva? Sime ioe Inesat ic: egoydeet 
Sant 


=. =e ee © ee oe ‘ 
4 pn a 1 






f 7 


195 


Mogi, K. (1979). Two kinds of Gaps; Pure and App]. Geophys. 
Went e2o1486% 


Molnar, P. (1979). Earthquake recurrence intervals and plate 
tectonics; Bul]. Seismol. Soc. Amer. 69, 115-133. 


Nur, A. (1978). Nonuniform friction as a physical basis for 
earthquake mechanics; Pure and App]. Geophys. 116, 
964-991. 


Ohtake, M. T. (1976). Search for precursors of the 1974 
Izu-Hanto-Oki earthquake, Japan; Pure and Appl. 
Geophys. 114, 1083-1093. 


Ohtake, M., T. Matumoto and G. V. Lathan (1977). Seismicity 
gap near Oaxaca, Southern Mexico as a probable 
precursor to a large earthquake; Pure and Appl. 
Geophys. 115, 375-385. 


Papazachos; B. (1974). On certain aftershock and foreshock 
parameters in the area of Greece; Ann. Geofis. 27, 
297 —5 115% 


Papazachos, B. (1975). Foreshocks and earthquake prediction; 
PECTONODRVSICS 28, 2132226. 


Papazachos, B., M. «Delibasis, N- Liapis, G. Moumoulidis, and 
G. Purcaru (1967). Aftershock sequences of some 
large earthquakes in the region of Greece; Ann. 
GeOkiste20,; 1593. 


Powelle, JevAv, and)S:. Ja Dudag( 19/75a)eeAgsteatiseicall study 
of earthquake occurrence; Pure and Appl. Geophys. 
113, 447-460. 


Powell, Je A.) andeS.ad. Duda (1975D)" sAngearthquake 
migration pattern of high statistical Significance; 
Pure and Appl. Geophys. 113, 461-466. 


Prozorov, A. G., and A. Dziewonski (1981). A method of 
studying variations in the clustering properties of 
earthquakes: application to the analysis of global 
seismicity; (in press). 


ve 
CC Gas ae 
‘ 
i] 
a | 
? nh 
a | | - = & = © 
| 






















~ 2. .a3ed to aiged ovt'){ Re 
| sett ~ST POT ee 


ot eteusegzeR . (272A 


wre ona 4S. Thi teataedaa” aa 
ps; 
ah wteSinunot «( BSt) cae 
MAF paypinEi sed es sui steS : ¥ 
/?=b oe 
" 
‘te, : M ,oo KITS 
’ one 1, eo 


“Ty .M shaesee 
Tt 
158% 
- 7 GT : iy ey 
; « ; ,endentagadg 
) 7S a4 x 1 
yi 


‘ae .€. ,@0i06eeq8% 
7 r = a4 e 
ve ar / it he ? ya? 


he Ets BS * » tise eee 


jrekee i ae (eC) } otan wwe 


[Vives teas Spas 


| , Of ,er tee | ae 
: 7 
J 5 . = f18 , A ols 5 Liew, 
os S.fsi ines sf#acpdize@ Fo 7s 


elSb-Tae (TER? 








196 


Quittmeyer, R. C. (1979). Seismicity variations in the 
Makran region of Pakistan and Iran: relation to 
great earthquakes; Pure and Appl. Geophys. 117, 
ee ae Or 


Rikatake, T. (1974). Probability of earthquake occurrence as 
estimated from crustal strain; Tectonophysics 23, 
Paste ke Shh Se 


Riki takcmelimGl OT oaye Stati sticssoreuleimate, strain of the 
earths crust and probability of earthquake 
occurrence; Tectonophysics 26, 1-21. 


Rikitake, T. (1975b). Earthquake precursors; Bull]. Seismol. 
Sows Tsien es) Wee) ee 


Rikitake, T. (1976a). Recurrence of great earthquakes at 
subduction zones: Tectonophysics 35, 335-362. 


Rikitake, T. (1976b). Earthquake Prediction, Elsevier, 
Amsterdam. 


Rikirakerl0lo/9)e)Classificatvonuvor earthquake 
precursors; Tectonophysics 54, 293-309. 


Sadovsky, M., I. Nersesov, S. Nigmatullaev, L. Latynina, A. 
Lukk, A. Semenov, I. Sembireva, and V. Ulomov 
(1972). The process preceding strong earthquakes in 
some regions of middle Asia; Tectonophysics 14, 
295233, 0a 


Sampson, R. J. (1975). Surface II Graphics System, Kansas 
Geological Survey, Lawrence, Kansas. 


Sauber, J., and P. Talwani (1980). Application of 
Keilis-Borok and McNally algorithms to earthquakes 
in the Lake Jocassee area, South Carolina; Phys. 
Earth Planet. Interiors 21, 267-281. 


Savage, J. C., and W. H. Prescott (1978). Asthenosphere 
readjustment and the earthquake cycle; JU. Geophys. 
Res. Soi cOl aos ols 











I ae | 


. | ) ae 
ct. ban bem SRETeI0* D> Ae 

oft Tel? Lawt; ie Roles? assiai on a4 : 
y "4 +ES2aRD ON I 168 - Vest. a ne of) 
EST -CRET We) | 


’ »! . 7 


ea 
lc es ae 


m2? es ears 9e 


wee Gs 4 1 


7 


,8iS7 ie tte 


i S29 Ors erp) 





MELO taal gr . Ge!) | inawelg? 4 E 


a a aol 





-. : 
i ae eee 
|. ip 


19Oe 


Sehouz Cm asad?) weAgphySsical interpretation fof the 
Haicheng earthquake prediction; Nature 267, 121-124. 


SGholamaC BH Pull BReESYkesmiandey y Pe Aggarwal ¥(1973) % 
Earthquake prediction: a physical basis; Science 
[Sie OS=5)h0¢ 


Sekiya, H. (1977). Anomalous seismic activity and earthquake 
Prediceron UCU AEehVSmeLasti. 25), 89094573. 


Shimazakiven. pander. Nakatas(1980) 2 )Time-predictable 
recurrence model for large earthquakes; Geophys. 
RES CCE Gems alo 252. 


Sieh, K. (1978). Prehistoric large earthquakes produced by 
slip on the San Andreas fault at Pallett Creek, 
California; JU. Geophys. Res. 83, 3907-3948. 


GongheS #eheneLeeAStiz and) J. Havskov (1981). Seismic gaps 
and recurrence periods of large earthquakes along 
the Mexican subduction zone: a reexamination; Bull. 
Seismol. Soc. Amer. 71, 827-843. 


Spath, H. (1980). Cluster Analysis Algorithms for Data 
Reduction and Classification of Objects, Ellis 
Horwood Publishers, Chichester, England. 


Suyehiro, S. (1966). Difference between aftershocks and 
foreshocks in the relationship of magnitude to 
frequency of occurrence for the great Chilean 
earthquake of 1960; Bull. Seismo]. Soc. Amer. 56, 
TS 522005 


Suyehiro, St vand He Sekiya (1972) .~Roneshocks and 
earthquake prediction; Tectonophysics 14, 219-225. 


Sykes, L. R. (1971). Afershock zones of great earthquakes, 
seismicity gaps, and earthquake prediction for 
Alaska and the Aleutians; JU. Geophys. Res. 76, 
SOZ1=3081" 


Sykes, LewRepeand RaiGssQuittmeyer (1981). Repeat times of 
great earthquakes along simple plate boundaries; (in 
press). 


a mo 








Ay ~oigniemyreies (aptabig a Otte)! pe 
Pa tant , Oo) US. we one ocetieee 





‘ - 
> ‘a 
later hanya pM onl! a ew . fossa: 4 r* 


; a 16n i.clceriqg’ & tho PSotires @seggccres 
of6-a208 taf 













. mt “ee | hich Ogee est CTR) HW Jevised 
<a, 2 wis (noise fbeia or a 


&@ ..8. tdesenaite 


~ rr + OC ee? of ,faika 


i a 


( ‘ ‘ ‘ | f 
$ 8 ; a je ry 
r > 
; = 


het ao Wer, SA) ote 
t e cs Th) Poe 
sar Oe 
Hg wo Sap Sop bd ‘Se cl) «2 <OFINRSSVE 
i yn? anc 6 3 ei 12e707 
= 5 9 ‘oy 0 iW YST 2 
: 
ar » \ ie yt FE WIPLITEs 
OOS cer ~ 


sie e* Steapig. ; “are oyLie¢ in Dri th 2k oti deyue 
SI PEW eh Tes wetethyiy steupis tee) si kee 
. 7 | oh ao. 1) 
Saas ¥ a pala fe 7 } 


i 
) 
w 
¥ 
> 
ah 
y 
r 
cy 
me 
7 
a 
_ 





a) 


a 





198 


Talwani, p. (1979). An empirical earthquake prediction 
model; Phys. Earth Planet. Interiors 18, 288-302. 


Thatcher mew mA i Leman, ande?, CG. Hanks (1977). Seismic 
Slip distribution along the San Jacinto fault zone, 
Southern Calrrornia,, and its: implications; Bul]. 
Geol esoc = Amen= 867 1140—1146.. 


Toksoz Mee Nee wate Ee eeohakal sand Ast. Michael (197.9)'. 
Space-time migration of earthquakes along the north 
Anatolian fault zone and seismic gaps; Pure and 
Appin Geopoys. 117, 125821270- 


Utsu, T.) (1972)-" Aftershocks and earthquake statistics; JU. 
Fac. Sci. Hokkaido Univ. Ser. VII, Geophysics N3, 
379-444", 


Utsu, T. (1972). Aftershocks and earthquake statistics; Ju. 
Fac. Sci. Hokkaido Univ. Ser. VII, Geophysics N4, 
1-42. 


Utsu, T. (1974). Space-time pattern of large earthquakes 
OccurcuImgqmom. ene Pacriic coast of ther Japanese 
molands mUOuUne es FOVS) Earth 22, 320-342. 


Vere-Jones, D. (1975). Stocastic models for earthquake 
sequences; Geophys. JU. Roy. Astron. Soc. 42, 
Sai oZo 


Vere: loneseeD."\ 1076). A branching, model™torecrack 
propagation; Pure and Appi. Geophys. 114, 711-725. 


Vere-Jones, D. (1978). Barthquake prediction: A 
Statisticiamisavicws UOUn-—haVS bance. Oo mm2 alnG. 


Weimu-G ce MSI Xl Urge eee chia Ope Sec naOn mance nm OU 
(1978). On seismic gaps previous to certain great 
earthquakes occurring in North China; Acta 
GEODhYS | Capo lWiGa 2lye 21 Sache 


Wesson Re ML ma OmWe LEE SWwOr theo 73)-ceismicity 
preceding moderate earthquakes in California; J. 
Geophys. Res. 78, 8527-8543. 


eS 
7 






Vee = Lavoe #3a4 ees ‘A eo TETE Be 
df newline | ttn 2 ga Giees i ie 


ie 7% es inmew <> T Ane. .wémelia iA * : Wf vasebnedn 
| ; ' sf, Ge ets on ~270 . 2 resin is fe | 8 
; ; 4 taut 4 z vn? 2 . 7. eae B44 
ee ee \tym 


im oe | fied ofl 1h ,. 9 Veneta 
; : . ‘ Ey ate oar) "3 28eS : 
ri ; i ad ; (otena 


| re Wg =e - } "} ha, 


aoe .( OVOP) £7 sea 
mo: osnol-e4ay 
LA aie) A s(?') 8 Sento waar 


Me 27 4 Tole eee ic a4 


“— We gol ts. » 50s taal coe a enol Ree | 1 
ease -ESh..( 45; TIS LS Ss. rarv 4° 2B! ‘salsa | 


oe As. 2. Gide a 
. sinatse no vt re ona 


ak ad 29 





122 


WittcCCoMmb auLmH me. Garmanyyeand D.)i.. Anderson (1973). 
Earthquake prediction: variation of seismic 
velocities before the San Fernando earthquake; 
Nalturem 81632-6035. 


WishantyeD selooOoMeNUMerIcal classification method for 
deriving natural classes; Nature 221, 97-98. 


Wishart, D. (1978). CLUSTAN User Manual (third edition), 
Program Library Unit, Edinburgh University, 
EGinbunghymococtland. 


Wood, M. D., and S. S. Allen (1973). Recurrence of seismic 
migration along the central California segment of 
the San Andreas fault system; Nature 244, 213-215. 


Wurm eer me Memey UCT arem WU ek. CaoOTe HemsGhen ye W. Hang, KK. Tain, 
and S. Lu (1976). Certain characteristics of the 
Haicheng earthquake (M = 7.3) sequence; Acta 
Geonnvs soinica 197 109-117. 


Wyss, M. (1973). Towards a physical understanding of the 
earthquake frequency distribution; Geophys. Uu. Roy. 
ASCCONMm— SOG. 3 l2,) 3402560. 


Wyss, M., and R. E. Habermann (1979). Seismic quiescence 
precursory to a past and a future Kurile Island 
earthquake; Pure and App]. Geophys. 117, 1195-1211. 


Wyss, M., and W. H. K. Lee (1973). Time variations of the 
average earthquake magnitude in Central California; 
Proc. Conf. Tecton. Prob. San Andreas Fauit System, 
Sieelabreleia! (lebuc, qeWieli 5 Aish, ace ier 


Yanashina, Ke, and) Yo Inoue (1979). A doughnut—shaped 
pattern of seismic activity preceding the Shimane 
earthquake of 1978; Nature 278, 48-50. 


YamasShina, K. (1981). Some empirical rules on foreshocks and 
earthquake prediction; (in press). 


Vequlalp a l)seM amo Cm sme lime Wl Oms(a0)/ 45) eens Lats tila! 
prediction of the occurrence of maximum magnitude 
earthquakes; Bull]. Seismol. Soc. Amer. 64, 393-414. 


pl) 


ees. Ru? tH 
ie i 


uo : 


cel A Th PP 3 mea * ane. i wt ao higeaie 











a. (oven? «Go a Se ee 
¥ iv tsar + hapapale . 
i its > G A apt “Epled ae 22 phir eri oF. 


a 
ae) Cta St 


eu ' , \SSSRi ee a2 ste iW 
6b 22 od u i. viveé ; a 


i f ‘ 
a CCl) .@ > Paedeae 
7 “2 ’ ; s vr 7 Ge) 74 : 
’ ip32vea sae - 
j { és j i 4 en Jor 
@ i 
td ' é *— 


A i, a “ies 


f ‘ ny re on 2274 


7 ra ba mare) = - , 2096 ' 
é cms. Ssaaieve 
fa ’ \ oor 
vi £4, res Ties J é. 











Bie: | (Lilie oepiee fo adegeee 


ne <8 Fa Ne vatep te depict: F 



















a aa 


vee 


ue ee” 








peas. Poca 2a einai 
ise seaereat Ce 8 oi 
met Ht bile i oo 


nem 


Aidan pied 
Being oka 


dheyeititi 
an 
et 
bre: 


fae ie 


er 
Satire 

ind oApaiin ded 
samen 


a 
ented #463 

» - ae 
eee aa 


Poe 
Peneiceteates 8 ee 


Fred 
dina 
é 


hy cat barrels a) 
ae 


jin0 
one ve 


ehedie4 


wa Tate 


roi 
APNE Y iat aad 
y 


1 hae 


al Habitat aah 4 ii 


hela dt 
Wipe ah we 
iN oe mena ir Hi 


g/L ae) 
Hoy epee 


Seep eden! 
Hecti 


PEE i 
cone Ap 


hii 
tte SSL ep ay 
Uh By aH Dns g Bey 
i Vy, Wilaaa ia 
SP RAN Ht EN ee 


a 
tipi Mire 
Aneel ry rand 


atti if 
th ‘ih Hye Ape 
peal vnita eae n wie 


ae 
Pua soll 
Pe) ihe id AES OBR WAU ed ati ee 
NY 


Vivhalba ante hedvh ie ealeai fl 
+ My 


polio? = ceed at 
peace 


ight aed wild 
Sa ligelAeto. peat nai 


A eteaieane pate ae 
ie 


Sei 


Sirap ascent 
Daag olde 
SEAT tee clint has 
pike 


ees yepsicusss 3 
Las NOREEN TEA 


“| 
erie 
Sheyea cabal 


ahaha Neg eNe A ALS 


a te sede noah 


Hie al 

me Peete) 

btbobe de dese seigrvnevantts pealens 
Had wa bbb) 9 than Hedise 


nh ii 

HOY stad 
itahenene heey 
ful S 


if 
jhe 


FE 
ahilon sa WERE nied 

Ling aaa as eit, 

ots 


tinnkdyay obi) 


eh Wells aide 6 dio 


7 


anche ay 


reais Bin 
inf WE 8) Ry pad “pheensoe) Frito 


W Blsreo ui sds here 
ise ¥EN sip eo tt caren Nh 
Prien : 
Pate ad aaieh 
eatin 
spe diez igh eo fea 


ff Webs poeta 
maT 


Erased yeahs 


trp PH 
tmues mearnsii ure Ti re 

Ah NPA Ta 

Jrastiney eo ey 

EY shigetetraintee 


Sauce 
ATOR 
Sie er ane AY 


smear? Mee: 
SENS 


Sen anb rer a 
Sure” Phin fy raty mt 
Paap pleas setae “ 


esyegreste BSE Hr 
aia RU pwr ubtes 
Cehcacae leita 
sus vor 
Seance taba a a akties 
min ys se eas Tee oat ft 
ae ny las fadienel 
peebteett ob ape LN Hed nL ae 


ah 
austerity: 
) 


eas mae 
ae rs 


or 


mye teacahy 
Sarid ees 


ACen ini Chena 
bad rasyiaissla ye iat 
pica iT : 
sara pyee 
ies 


+H 
etal 
ae 


Hoy 
a rah aa 


) Tiere mines 
ae mad ry 
Rakobad 


FFU rae aweMs Wave 
Vaid sweety sits 
SAU tae hy 


Mey 


Sy ‘ 
HEBER! 
avirdey 


M ra 
iar sy 
yh 


i 
wensen aie 
Bins 

ie 


Burhep ade alert 
t agai 4) Tatas 


eitasie 


spy 


4 
DARE VA aid Neon te 
f ce ra 


Leet ih 

fed ANY 

EAU Maat) Ri Red Well d) nh aelcnah 
a 


eesacpte UR net 


val 


er to Path ER ae heh 


‘ 
ae if With yy? ipa 


Gepey Su Mellnyy 


pessrsyave nest 


Ns 
j 


Rye) 
naa 


Heber y nell beth Babin? 


Sasa pS tis 
oh WHA Wr, 


% 

Wainbisa wae 1 
me es J Tie \ 
ENR NST a hate 
YEW Kea res 


Pandey abe es 


bawaes) Bis 


DyUadets RAMUS 
Anh ee) 


4 
VNR Uhatah (Nt cea 
Nt Hy a 
TANG 


4 

i} 

a Pea flenikye hive peg as 
y Ay UN 

we visit 


WE VSN 


We als tit y 
Heh) ty ined ak oh t Wied 3 em Ayes \ 


Ht 


desis 


i; 
Ants a Wty 


asta 
ii As ae he 


buteisesbpund 
HW NEVO Maths 


bed ot 
at Si wy Ud 
penn 


bebe Ad 
PAM EN wl 
avis 
Peak tata ity 
boasted 
I 


CEE 
iN A Bela hey 
Vado, howe h 
ie sonia iy 
uta Ni 
wR ahd oh Huan 
YONG DeWine 
Naina tah 
iia yl iE 


LN ont tte 


Diciia neva lie) it! 





